Random Matrices and Chaos in Nuclear Physics 
Part 1: Nuclear Structure 
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The authors review the evidence for the applicability of random-matrix theory to nuclear spectra. 
In analogy to systems with few degrees of freedom, one speaks of chaos (more accurately: quantum 
chaos) in nuclei whenever random-matrix predictions are fulfilled. An introduction into the basic 
concepts of random— matrix theory is followed by a survey over the extant experimental informa- 
tion on spectral fluctuations, including a discussion of the violation of a symmetry or invariance 
property. Chaos in nuclear models is discussed for the spherical shell model, for the deformed 
shell model, and for the interacting boson model. Evidence for chaos also comes from random- 
matrix ensembles patterned after the shell model such as the embedded two— body ensemble, the 
two— body random ensemble, and the constrained ensembles. All this evidence points to the fact 
that chaos is a generic property of nuclear spectra, except for the ground-state regions of strongly 
deformed nuclei. 



PACS numbers: 24.60.l_z,25.70.Gh,21.60.Cs,21.60Ev 

I. INTRODUCTION 

In the 1970s chaos became a household word of physi- 
cists. Chaos had been known to mathematicians, astro- 
physicists and other specialists since the early years of 
the century. It took the advent of the personal com- 
puter to make chaos generally known: The exponential 
divergence of the trajectories of a chaotic system with 
Hamiltonian dynamics could easily be simulated. Many 
Hamiltonian systems, especially those with few degrees 
of freedom, have been analyzed since. The subject has 
been reviewed, for instance, by Gutzwiller (1990). 

In the wake of this development, physicists in several 
fields became interested in quantum manifestations of 
classical chaos ("quantum chaos"). Again, this caused 
a flurry of activity, especially in the study of systems 
with few degrees of freedom, see McDonald and Kaufman 
(1979), Casati et al. (1980), and Berry (1981). The work 
culminated in the Bohigas-Giannoni-Schmit (1984) con- 
jecture. The conjecture connects the spectral fluctuation 
properties of quantum systems which are chaotic in the 
classical limit with predictions of random-matrix theory 
(RMT). A summary of these developments is given by 
Haake (2001). 

Independently of these developments and preceding 
them, RMT had been developed in the framework of nu- 
clear physics by Wigner and Dyson - see papers in the 
compilation by Porter (1965). Data accumulated in the 
1960s and analyzed in the 1980s provided evidence that 
nuclear spectra follow RMT predictions. The wide in- 
terest enjoyed by RMT in the 1980s among practitioners 
of quantum chaos had repercussions on nuclear physics: 
Nuclei are many-body systems, and chaos manifests it- 
self here in ways different from those of few-degrees-of- 
freedom systems. At the same time, nuclei are paradig- 
matic for the applicability of RMT to fermionic many- 



body systems and for the occurrence of chaos in such 
systems. Thus much work was done to establish RMT 
and to analyze and interpret quantum chaos in nuclei. 

The present review is the first part of a two-part se- 
ries dealing with RMT and chaos in nuclear physics. It 
focuses on spectral properties of nuclei, while the sec- 
ond part will deal with RMT and chaos in nuclear reac- 
tions. Topics such as compound-nucleus scattering, Er- 
icson fluctuations, isobaric analog resonances, and parity 
violation in epithermal neutron scattering are not dealt 
with here. The paper is intended as an introductory re- 
view to the field. It is aimed mainly at two groups of 
physicists, those who work on quantum chaos in fields 
different from nuclear physics, and nuclear physicists who 
wish to learn about RMT and chaos in nuclei. We accord- 
ingly assume no prior knowledge of RMT and chaos nor 
do we assume any detailed knowledge of nuclear physics. 
We aim at giving a comprehensive survey which is fo- 
cused on concepts and illustrative examples while deriva- 
tions and formulas are kept to a minimum. The last 
comprehensive review of the field has been given in Re- 
views of Modern Physics over 25 years ago by Brody et 
al. (1981); a short review was later published by Bohigas 
and Weidenmuller (1988). Wherever possible we have 
avoided giving a large number of references in favor of 
citing review articles: Readability of the article was our 
primary concern, followed by completeness. 

In Section [II] we motivate RMT and introduce those 
concepts of RMT which are frequently used in nuclear 
physics. We pay particular attention to the Gaussian or- 
thogonal ensemble of random matrices (GOE). We estab- 
lish the connection between RMT and quantum chaos. In 
Section irm we describe those applications of RMT to nu- 
clear spectra which do not make use of specific nuclear- 
structure concepts. The comparison between RMT pre- 
dictions and spectroscopic data is used to establish the 



domain of applicability of RMT to nuclear spectra. Spe- 
cial attention is devoted to the breaking of isospin sym- 
metry, and to a test of time-reversal invariance. The 
role of RMT and chaos in nuclear models is described in 
Section [TV] We focus attention on the two most impor- 
tant nuclear-structure models, the spherical shell model 
and the collective model in two of its versions, but also 
mention a number of other applications of RMT. Nuclei 
(and other fermionic many-body systems) are governed 
by the mean field (in nuclei: the shell model) and a resid- 
ual interaction dominated by two-body forces. In such 
systems, the structure of the Hamiltonian is very differ- 
ent from that of a typical GOE matrix. That difference 
has given rise to a number of random-matrix ensembles 
which are closer in structure to the mean-field approach 
than is the GOE. These are treated in Section [V] Sec- 
tion |VI] contains a summary and conclusions. 



II. RANDOM MATRICES 
A. Why Random Matrices? 

Random Matrices were introduced into Nuclear 
Physics in the 1960s by Wigner (See Wigner's papers 
in Porter (1965)). That step was preceded (and, in our 
view, probably motivated) by Bohr's insight that nuclei 
are systems of great complexity. It is useful to recall the 
arguments that led Bohr (1936) to this insight. 

Experiments in the 1930s, especially by Fermi and his 
group in Rome on neutron scattering by light nuclei, had 
revealed the existence of numerous narrow resonances 
(Fermi et al, 1934, 1935). We show in Fig. □] data of 
a similar type taken in the 1950s by Rainwater and his 
group at Columbia University (Garg et al., 1964). That 
group used time-of-flight spectroscopy of slow neutrons 
to measure the total neutron cross section on a number 
of heavy even-even nuclei (nuclei with even numbers of 
protons and neutrons). The cross section versus neutron 
energy E n shown in the panels of Fig. [1] for the target 
nucleus 232 Th displays narrow resonances with widths < 
1 eV and spacings of about 20 eV. The target nucleus 
232 Th has spin zero and positive parity; the incident 
slow neutrons carry zero angular momentum and have 
spin 1/2. Therefore the resonances all have spin/parity 
1/2+. These resonances correspond to excited states of 
the "compound nucleus" 233 Th with an excitation en- 
ergy slightly above the neutron separation energy of 4.786 
MeV (the "neutron threshold"). The number of reso- 
nances observed in each compound nucleus was limited 
by the resolution of the spectrometer and was never much 
larger than 200. Similar data on proton resonances at 
the Coulomb barrier in lighter nuclei were later taken by 
the Triangle Universities group (Wilson et al, 1975). To- 
gether these data form what has been called the "nuclear 
data ensemble" by Haq et al. (1982) and Bohigas et al. 
(1983). 

Bohr argued that the existence of numerous narrowly 




FIG. 1 The total neutron cross s ectio n on Th versus neu- 
tron energy E n in eV. From Ref. (jl46h as reproduced in Bohr 
and Mottelson (1969), Vol. 1, p. 178. 



spaced and narrow compound-nucleus resonances was in- 
compatible with independent-particle motion and was 
due to strong nucleon-nucleon interactions. Indeed, as- 
suming an independent-particle model with a nuclear ra- 
dius of about 5 fm and a potential well depth of several 
ten MeV, one finds that the single-particle states have 
a typical spacing of several hundred keV and widths of 
the order of ten keV or larger, in complete disagreement 
with the data. To account qualitatively for the data, 
Bohr proposed his "compound-nucleus model" (Fig. [2]): 
The incident nucleon carries kinetic energy (as indicated 
by the billiard cue), collides with the nucleons in the tar- 
get and shares its energy with many nucleons. In units 
of the time for passage of the nucleon through the nu- 
clear interior, it takes the system a long time until one of 
its constituent nucleons acquires sufficient energy to be 
re-emitted from the system. 

Bohr's idea that the nucleus is a complex, strongly 
interacting system was adopted by the community and 
held sway until the discovery of the nuclear shell model 
in 1949. Bohr's idea almost certainly motivated Wigner 
to introduce random matrices. To explain the spirit of 
the approach, we focus attention on nuclear levels with 
the same quantum numbers (total spin J, parity II, and, 
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FIG. 2 Bohr's wooden toy model of the compound nucleus. 
From Nature (1936). 



at least in light nuclei, total isospin T) and ask: Can 
we identify generic spectral properties of a system with 
strong interactions? Fig. [3] shows six spectra, all having 
the same total number of levels, and spanning the same 
total energy interval, and therefore having the same av- 
erage level spacing. The spectra differ only in the way 
the spacings between neighboring levels are distributed. 
For the one-dimensional harmonic oscillator (the right- 
most spectrum), all spacings are identical. The spacing 
distributions differ more and more from a delta function 
as we go ever more to the left. The random-matrix ap- 
proach characterizes spectra by their fluctuation proper- 
ties: The distribution of spacings of nearest neighbors is 
the first and obvious measure for spectral fluctuations. 
It is referred to as the nearest-neighbor spacing (NNS) 
distribution. There are other measures such as the cor- 
relation between nearest spacings, between next-nearest 
spacings, etc. Some of these are introduced below. 

To implement this approach, we need to develop a sta- 
tistical theory of spectra. Random matrices provide the 
tool to do so. Instead of considering the actual nuclear 
Hamiltonian (which was not known in the 1950s) we con- 
sider an ensemble of Hamiltonians (each given in matrix 
form). The ensemble is defined in terms of some prob- 
ability distribution for the matrix elements, hence the 
name random matrices. The ensemble is chosen in such 
a way that the member Hamiltonians incorporate generic 
features. The spectral distribution functions are calcu- 
lated as averages over the ensemble and are compared 
with the actual fluctuation properties of nuclear spectra. 

Canonical random-matrix theory (RMT) as developed 
by Wigner and Dyson - see the compilation by Porter 
(1965) - classifies systems by their symmetry properties. 
Nuclei are invariant under time reversal. The matrix rep- 
resentation of the nuclear Hamiltonian can accordingly 
be chosen real and symmetric. The random-matrix en- 
semble which is considered almost exclusively in the se- 
quel is therefore an ensemble of real and symmetric ma- 
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FIG. 3 Six spectra with 50 levels each and the same mean 
level spacing. From right to left: The one-dimensional har- 
monic oscillator, a sequence of zeros of the Riemann Zeta- 
function, a sequence of eigenvalues of the Sinai billiard (see 
Section IH.F|I , a sequence of resonances seen in neutron scat- 
tering on 166 Er, a sequence of prime numbers, and a set of 
eigenvalues obeying Poisson statistics (see section IlLFfr . From 
Bohigas and Giannoni (1984). 



trices. 

The random-matrix approach does not aim at calculat- 
ing individual spectra and at comparing them with data. 
Rather, one determines the joint probability distribution 
of the eigenvalues and from here calculates certain spec- 
tral fluctuation measures such as the NNS distribution 
as averages over the ensemble. RMT contains one (or, in 
the general case, a number of) input parameter(s). In the 
case of spectral fluctuations, that input parameter is the 
average nuclear level spacing. The fluctuation measures 
predicted by RMT are scaled by the average level spac- 
ing and, thus, parameter-free. If the observed spectral 
fluctuation properties agree with RMT predictions, and 
if no further information on the system is available, one 
concludes that the system is generic. This implies that no 
information beyond the average nuclear level spacing can 
be deduced from the available spectral information. If, 
on the other hand, the data do not agree with RMT pre- 
dictions, this indicates that the spectrum is not generic 
and that the available spectral information may be used 
to deduce further properties of the system. The harmonic 
oscillator in one dimension is a case in point. 

The random-matrix approach to spectral fluctuations 
(and to other properties of complex systems) has some 
similarity to classical thermodynamics. There one is also 
interested in a generic description of systems in terms of 



a few parameters. These parameters (specific heat, mag- 
netic susceptibility, etc.) are system-specific, but within 
the framework of classical thermodynamics, need not 
be determined from the system's Hamiltonian. In that 
sense, classical thermodynamics and random-matrix the- 
ory are phenomenological theories that do not refer to an 
underlying system-specific Hamiltonian. The random- 
matrix approach differs fundamentally from the dynam- 
ical approach used in most fields of physics where one 
integrates the equations of motion and fits a few param- 
eters of an (otherwise known) Hamiltonian to the data. 
Similar to classical thermodynamics, the random-matrix 
approach has been applied to many systems beyond nu- 
clear physics (Guhr et ai, f998). 



B. The Gaussian Orthogonal Ensemble 



In Section Hl. Bl and in Section Hi. CI we define the Gaus- 
sian orthogonal ensemble (GOE) and collect and inter- 
pret a number of results for this ensemble. We also in- 
troduce the Gaussian unitary ensemble (GUE). Proofs 
may be found in the books by Porter (1965) and Mehta 
(2004). For the GOE we consider real and symmetric 
Hamiltonian matrices H in a Hilbert space of dimen- 
sion N. With /x, v = 1, . . . , N the matrix elements obey 
H^v = H vfl = H*. For realistic systems Hilbert space 
is infinite- dimensional, so we consider the limit N — ► oo 
in the sequel. The ensemble is defined in terms of an in- 
tegration over matrix elements. The volume element in 
matrix space 



d[H] = n <!#„„ 



(i) 



^A<L> 



is the product of the differentials dH^ w of the indepen- 
dent matrix elements (i.e., of the matrix elements not 
connected by symmetry). The ensemble is defined by 
the probability density V{H) of the matrices H, 

V{H) d[H] = M) exp | - A T race(tf 2 ) j d[H] . (2) 

Here A/"o is a normalization factor and A a parameter. 
This parameter defines the average level density (see 
Eq. (|I4p below). In applications of the GOE to data, 
A is determined by the empirical average level density. 
The spectral fluctuation properties of the GOE are then 
predicted in a parameter-free fashion. 

The Gaussian weight factor is a cutoff that ensures 
convergence of the ensemble averages of observables for 
large values of the integration variables. We use the sym- 
metry of the matrices to write the trace in the exponent 
as E P <„ 2H lu + E P H lv Tnen the probability density 

V(H)d[H] = Mol[exp{-^LH^}dH^ 



N 



J\cM-7^H%}dH P ° (3) 



is a product of terms each of which depends only on a sin- 
gle matrix element. Therefore the GOE has the following 
properties: The independent matrix elements are uncor- 
rected Gaussian-distributed random variables with zero 
mean value and a second moment given by 
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(4) 



Here the overbar denotes the ensemble average. Defining 
the GOE by these properties is equivalent to the defini- 
tion ©. 

While the form of the probability measure in Eq. © 
is fixed by symmetry requirements, the Gaussian cutoff 
in that equation seems completely arbitrary. However, 
using plausible assumptions one can actually derive that 
factor. Rosenzweig and Porter (1960) have shown that 
the distribution |(5J) is obtained when one assumes that 
the ensemble is orthogonally invariant, and that matrix 
elements not connected by symmetry are statistically in- 
dependent. And Balian (1968) has derived the distribu- 
tion ((2|) from a maximum entropy principle. 

In the GOE every state in Hilbert space is connected 
to itself and to every other state by a matrix element 
of H . Since all non-diagonal matrix elements have the 
same first and second moments, every state is coupled to 
all other states with equal average strength. This results 
in level repulsion between any pair of levels, and in a 
complete mixing of states in Hilbert space. The impor- 
tance of such coupling is seen when we consider a more 
general ensemble with probability density 

V a (H)d[H] = ^oIlM-^}^ 



x 



N 



L[ ex P{- T^2 H 2 p JdH pa (5) 



p<o 
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p<cr 



where the positive parameter a ranges from zero to one. 
For a = all non-diagonal elements vanish, and the 
ensemble (J5J) consists of diagonal matrices with inde- 
pendent, Gaussian-distributed diagonal elements. The 
shape of the average spectrum is Gaussian, there is no 
level repulsion, and the spectral fluctuations are Poisso- 
nian (see Section lll.Fp . For a = 1, the ensemble coincides 
with the GOE. For values of a between these two limits, 
the shape of the spectrum and the spectral fluctuations 
interpolate between those two limiting cases. Significant 
mixing between levels occurs when the mean-square mix- 
ing matrix element H 2 with ji ^ v is roughly equal to 
the square of the mean level spacing. Taking for the latter 
the GOE value d = tt\/N at the center of the semicircle 
(see Eq. (fl4l below), we find that significant mixing oc- 
curs when a is of order 1/yN. We see that for N — > oo, 
mixing sets in as soon as a differs from zero. These ob- 
servations are used in Sections IIII.D.ll and IIII.D.4I 

Reality and symmetry of the matrices H pL , are pre- 
served under orthogonal transformations of the basis. 
The ensemble @ is accordingly chosen in such a way 



that it is invariant under such transformations: With 
each matrix H belonging to the ensemble, all matrices 
obtained from H by orthogonal transformations also be- 
long to the ensemble. As a consequence, there does not 
exist a preferred direction in Hilbert space, and the en- 
semble is generic. Because of that invariance and the 
Gaussian cutoff, the ensemble is referred to as the Gaus- 
sian orthogonal ensemble of random matrices. 

Instead of the N(N + l)/2 integration variables H^ 
with /j, < v used in Eq. |(5J), we may use the N eigenvalues 
E„ of the matrices H and the N(N — l)/2 generators of 
the orthogonal transformation O which diagonalizes H . 
Then the volume element AH takes the form 



dH = dO]\\E IJt -E v \\\dE p 



(G) 



p<v 



The factor dO stands for the Haar measure of the or- 
thogonal group in N dimensions. (The Haar measure is 
the unique invariant measure that can be assigned to ev- 
ery compact group and that is used to define integrals 
over that group (Conway, 1990)). The probability den- 
sity V{H) takes the form 



T(H)d[H] = A^odOexp -_££* 



N 



Jl\E p -E a \]JdE, 



(7) 
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The right-hand side of Eq. (O is the product of two fac- 
tors. One factor depends only on the eigenvalues and 
the other, only on the diagonalizing matrices. It follows 
that the eigenvalues and the eigenvectors of the matri- 
ces H are uncorrelated random variables. In the limit 
N — > oo, the projections of the eigenvectors onto an ar- 
bitrary direction in Hilbert space have a Gaussian dis- 
tribution. Some properties of the eigenvalue distribu- 
tion can be read off directly from Eq. ([7]). The factor 

n 



p<v 



\Ep — E v \ stems from the volume element in ma- 
trix space and reflects the orthogonal invariance of the 
ensemble. It causes the probability density for the eigen- 
values to go to zero as two eigenvalues approach each 
other. This is a manifestation of level repulsion, a basic 
feature of quantum mechanics. 

Another property of the GOE is displayed when we 
write the probability density in the form 

T{H)d[H]=JV dOY[dE t t 
p- 

xexpj-^^i^ + EH^-^l} -(8) 

^ V p<(T ' 

For the interpretation of the eigenvalue distribution in 
Eq. ||5J) we use an analogy to classical statistical mechan- 
ics. We consider the eigenvalues E p as position coordi- 
nates of N particles in one dimension. The probability 
density of the eigenvalues in Eq. ([5]) then has the form of 



the canonical partition function (integrated over the mo- 
mentum variables) of a gas of N classical point particles 
with repulsive two-body interactions ("Coulomb gas") 
moving in a common harmonic oscillator potential at in- 
verse temperature (3 = 1. The particles will tend to keep 
apart as much as is consistent with the overall harmonic 
oscillator potential. This property leads to the "spectral 
stiffness" of the GOE discussed below. 

Besides the GOE there exist two more canonical 
random-matrix ensembles: The Gaussian unitary ensem- 
ble (GUE) and the Gaussian symplectic ensemble (GSE). 
The GUE is the Gaussian ensemble of Hermitian (but not 
necessarily real) matrices. This ensemble plays a role for 
systems which are not time-reversal invariant. If that 
invariance does not hold, the Hamiltonian matrix is Her- 
mitian but cannot in general be chosen real and symmet- 
ric. In nuclear physics, the GUE is used for tests of time- 
reversal invariance. It is also used as a theoretical testing 
ground because the calculation of ensemble averages over 
observables is typically simpler for the GUE than for the 
GOE. For these reasons, we briefly introduce the GUE in 
the next paragraphs. The GSE applies to systems with 
half-integer spin which are invariant under time reversal 
but which are not rotationally invariant. The GSE does 
not apply to nuclei directly (see, however, Lombardi et 
al., 1994) and is not discussed in this review. In addition 
to the three canonical ensembles introduced by Dyson 
(see Porter, 1965) and often distinguished by the label 
(3 with $ = 1 for the GUE, (3 = 2 for the GOE, (3 = 4 
for the GSE, there exist seven more random-matrix en- 
sembles which are defined in terms of invariance require- 
ments (Altland and Zirnbauer, 1997) ; their construction 
is based upon Cartan's classification of Lie groups (see 
Chevalley, 1946). Some of these ensembles play a role in 
the low-energy behavior of quantum field theories and 
find application in lattice QCD, the discretized form of 
quantum chromodynamics. Others relate to the scat- 
tering of electrons off the interphase between a normal 
conductor and a superconductor. 

For the GUE the independent variables are the real and 
the imaginary parts of the complex elements H^ E of the 
Hamiltonian (a matrix of dimension N). The invariant 
measure has the form 

d[ F GuE ] = n d[KF«™]d[^« uE ] n dH^ . o) 

With this definition, the equation for the probability den- 
sity of the GUE is similar to Eq. ^ and reads 

V(H GVE ) d[ff GUE ] 

= A/- cxp[ - _^Trace(ff GUE ) 2 j d[H GVE ] .(10) 

The GUE is invariant under unitary transformations of 
Hilbert space. The real and imaginary parts of the matrix 
elements are uncorrelated random variables with equal 
Gaussian probability distributions centered at zero. The 
factors in the exponent are chosen in such a way that the 



second moments have the values 
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(11) 



Level repulsion in the GOE is linear, see Eqs. ([7]) and 
(I19p . This is a consequence of orthogonal invariance. In 
the GUE, the transformation to eigenvectors and eigen- 
values as new integration variables involves a unitary 
transformation U and yields 



P(i7 GUt )d[ff G H = A/" dWexp 



N 
2A2 



E^ 



xl[(E p -E a ) 2 Y[dE„ .(12) 

p<a v 

Here dli denotes the Haar measure of the unitary group 
in N dimensions. Instead of the factor \E p — E a \ occur- 
ring in Eq. (0, Eq. ([12]) contains the factor (E p - E a ) 2 . 
As a result, level repulsion for the GUE is quadratic. 
That difference between GOE and GUE is easily under- 
stood: In the GOE, the coupling of any pair of levels is 
described by a single parameter, the real coupling ma- 
trix element. For two levels to have a small spacing, the 
value of that parameter must be small. In the GUE, the 
coupling is described by two parameters, the real and the 
imaginary parts of the coupling matrix element. For two 
levels to have a small spacing, both parameters must be 
small, and the probability of small spacings is reduced 
accordingly. 

As concerns the GUE analog of Eq. ©, the form of 
the volume element in matrix space leads to a different 
inverse temperature [3 = 2, while for the GSE we have 
(3 = 4. The "Dyson parameter" (3 with (3 = 1,2,4 is 
often used to label the three canonical random-matrix 
ensembles GUE, GOE, GSE. 



C. Properties of the GOE 

1. Average Level Density 

A central property of the GOE is the mean level density 
p(E) , a function of the energy E. It is defined as 



p(E)=J2t(E-E„) 



(13) 



and, for N — > oo, given by 



**>-&A-(§)' 



(14) 



The average spectrum extends from — 2A to +2A. This 
confinement of the spectrum to a finite stretch of the en- 
ergy axis is a consequence of the Gaussian cutoff (or, for 
that matter, of any other sufficiently strong cutoff factor) , 
and of the factor N in the exponent of Eq. @. When 



plotted versus E/(4X), p(E) has the shape of a semicircle. 
(That shape is specific for the Gaussian cutoff). That is 
why Eq. (Tl"4"|) is often referred to as "Wigner's semicircle 
law" . The factor N on the right-hand side of Eq. (TT4")) 
ensures that p(E) is normalized to the total number of 
levels. The mean level spacing d(E) is defined by 



d{E)=p-\E) 



(15) 



and tends to zero as N — > oo because we fit a spectrum 
of N eigenvalues into a finite energy interval of length 4A. 
At the center of the spectrum, we have d(0) — irA/N. 



2. Universality 

The form of the spectrum is due to the Gaussian cut- 
off factor. That form is obviously totally unrealistic: 
Hardly any real physical system possesses such a spec- 
trum. While reality and symmetry of the matrices -ff M „ 
reflect time-reversal invariance and are thus a conse- 
quence of quantum theory, the Gaussian cutoff is not, al- 
though the arguments of Rosenzweig and Porter (1960) 
and of Balian (1968) lend some plausibility to its use. 
The Gaussian cutoff is preferred from a practical point of 
view, of course, because of the ease with which Gaussian 
integrals can be performed. But the GOE is physically 
interesting only if it furnishes information which is inde- 
pendent of the form of the cutoff factor. That property 
is guaranteed by the universality of the GOE. 

In using the GOE we are usually not interested in the 
overall shape of the spectrum. Interest rather focuses on 
local spectral fluctuation properties such as the NNS dis- 
tribution or correlations between level spacings. These 
are predicted in a parameter-free fashion. That means 
that all local spectral fluctuation properties are functions 
of a dimcnsionless parameter s, the ratio of the actual 
level spacing and the mean level spacing. Local spec- 
tral fluctuations characterize properties of the spectrum 
on an energy scale which in the limit N — > oo is negli- 
gibly small compared to the length 4A of the spectrum. 
On that scale, the spectral fluctuations are universal: As 
functions of the parameter s they have the same form 
for both the GOE and all non-Gaussian cutoff factors, 
as long as the latter are orthogonally invariant and con- 
fine the spectrum to a finite singly-connected piece of the 
energy axis (Hackenbroich and Weidcnmuller, 1995). 

Non-Gaussian cutoffs which obey that proviso modify 
the overall shape of the spectrum. In fact, for any given 
form of the spectrum it is always possible to find a cutoff 
factor such that the resulting random matrix ensemble 
has an average spectrum of that form. The local fluctu- 
ation properties are unaffected by such a choice: In the 
limit N — > oo, the local fluctuation properties separate 
from the global spectral properties and become universal. 



3. Ergodicity 

Theoretical predictions of the GOE are obtained as av- 
erages over the ensemble. How can we compare such pre- 
dictions in a meaningful way with data which, after all, 
are taken from a physical system with a single Hamilto- 
nian (and not from an ensemble of Hamiltonians)? That 
question is answered by the property of ergodicity of the 
GOE. Spectral data on a given system can be used to 
calculate spectral measures such as the mean level spac- 
ing or the NNS distribution as running averages over the 
spectrum. We denote such a running average by angu- 
lar brackets. We would like to ascertain that O = (O) 
holds true for all members of the ensemble and for all ob- 
servables O that describe local spectral properties. That 
equation cannot be proved in general because there is no 
way to evaluate (O) in the framework of the GOE. It is 
possible, however, to prove the slightly weaker statement 
(Brody et al, 1981) 



(O-(O)Y = 



(16) 



The proof is made possible because all terms on the left- 
hand side are ensemble averages. The statement says 
that for almost all members of the ensemble (with the 
exception of a set of measure zero; the measure being de- 
fined in Eq. ([T])) the running average of an observable O 
(calculated for a single member of the ensemble) is equal 
to the ensemble average of the observable. That prop- 
erty is referred to as ergodicity. The name derives from 
the formal similarity of the statement with ergodicity in 
classical statistical mechanics (equality of phase-space 
average and time-average along a single trajectory). 



4. Information Content of GOE Spectra 

Eq. ([3|) shows that in the GOE, every state in Hilbert 
space is coupled to every other one by a Gaussian- 
distributed random matrix element: In the GOE all 
states in Hilbert space are completely mixed with each 
other. Choosing the parameters N and A and drawing all 
independent matrix elements from the resulting Gaussian 
distribution generates a random GOE matrix. Diagonal- 
izing that matrix yields a GOE spectrum. By construc- 
tion, that spectrum contains no information beyond the 
input parameters N and A. In particular, the spectral 
fluctuations are void of physical information. If the spec- 
tral fluctuations of an experimental spectrum agree with 
GOE predictions, and if there is no further information 
on that system, then the spectral data alone cannot be 
used to extract any physical information on the system 
beyond the mean level density. 

That conclusion is also reached when we ask: How 
many pieces of spectral data are needed to determine 
the underlying Hamiltonian HI In the case of a GOE 
spectrum, counting shows that we need all N eigenvalues 
and all N orthonormal eigenfunctions to determine the 
N(N + l)/2 independent matrix elements of H. This 



must be compared with the usual dynamical approach to 
physical systems where the Hamiltonian is given in terms 
of a few (say n) parameters. Then n pieces of data suffice 
to determine the Hamiltonian. Further data can be used 
to check the consistency of the underlying theory. 



D. GOE Fluctuation measures 

1. Porter-Thomas distribution 

We recall that in the GOE eigenvalues and eigenfunc- 
tions are uncorrelated random variables. For N — ► oo, 
the projections of the eigenfunctions onto an arbitrary 
vector in Hilbert space have a Gaussian distribution cen- 
tered at zero. Therefore the squares ip 2 of such projec- 
tions have a % 2 -distribution with one degree of freedom. 
We introduce the variable 



v = i> 2 1 ip 2 



(17) 



The resulting distribution is also known as the Porter- 
Thomas distribution and has the form 



P{y) 



i 



: exp(-2//2) . 



(18) 



The function P{y) is given in terms of the mean value T = 
ip 2 . That parameter is an input parameter which is not 
predicted by random-matrix theory. The distribution 
can be checked experimentally: Transition probabilities 
of nuclear levels to a fixed final state and decay widths 
to a fixed channel are proportional to squares of matrix 
elements containing the nuclear wave functions. These 
matrix elements can be read as projections of the wave 
functions onto a particular vector in Hilbert space. 

It may happen that the mean value T undergoes a sec- 
ular variation. This is the case, for instance, for doorway 
states. Then it is necessary to "unfold" the fluctuations 
by scaling the intensities properly, see the end of Sec- 
tion |ILG] 



2. Nearest-Neighbor-Spacing Distribution and A3-Statistic 

It takes substantial theoretical effort to work out the 
spectral fluctuation measures in the GOE. That is not 
described here. We confine ourselves to introducing 
two fluctuation measures that have found wide applica- 
tion in the analysis of experimental data: The nearest- 
neighbor-spacing (NNS) distribution P(s) and the A 3 
statistic due to Dyson and Mehta. These are obtained in 
the limit N — > oo. Prior to using these measures for data 
analysis, it is neccessary to "unfold" the experimental 
spectra, see Section Hi. D. 31 

The NNS distribution P(s) depends on s, the ratio of 
the actual level spacing and the mean level spacing d. It 
cannot be given in closed form. An excellent approxima- 
tion due to Wigner is known as the Wigner surmise, 



P(s) = -sexp(-7rs 2 /4) 



(19) 



The linear increase with s for small s is due to GOE level 
repulsion as displayed in Eq. ([7]) . Universality shows that 
the Gaussian falloff is not related to the Gaussian cutoff 
factor defining the GOE and simply accounts for the fact 
that very large spacings are unlikely to occur. The exact 
expression for P(s) was first derived by Gaudin (1961). 
P(s) is displayed in Fig.SJ 



would grow linearly with L. In actual fact Y?JL) is, for 
large L, proportional to InL. The slow growth indicates 
that large spacings and small spacings do not follow each 
other at random but almost alternate and reflects the 
"stiffness" of GOE spectra, see the text below Eq. ((5J. 
For the three canonical ensembles, the number variance 
is shown in Fig. [5] The number variance is seldom used 
in nuclear physics because it fluctuates too strongly, and 
one uses the A3-statistic by Dyson and Mehta instead. 
The latter is defined by 
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FIG. 4 The nearest-neighbor-spacing (NNS) distribution of 
the GOE (solid line) versus s, the ratio of the actual level spac- 
ing and the mean level spacing. For the sake of comparison we 
also show the NNS distributions for the GUE (dashed line) 
and the GSE (dotted line). The parameter /3 is the Dyson 
index with 13 = 1, 2,4 for GUE, GOE, GSE, respectively. 



Aa(£) = 

min Qj6 -/ / dE" (Af(E')-a-bE 



(21) 



We integrate the ensemble average of the square of the 
difference between the number staircase function and the 
straight line (a + bE') over an energy interval, divide by 
the length L of that interval, and minimize the result 
with respect to the parameters a and b of the straight 
line. The angular brackets denote an average over the 
initial point Eq. It can be shown that A^(L) can be 
written as an integral over the number variance £|(Z). 
Therefore A3 is much smoother than So(L) and is better 
suited for data analysis. Similar to S|(L), A^(L) grows 
logarithmically with L. For large L, 



A 3 (L) « -4{lnL- 0.0687} 

7H 



(22) 



Similar to E 2 , the A3-statistic reflects the stiffness of 
GOE spectra and is often referred to as "spectral stiff- 
ness". Fig.[7]shows A 3 (L) versus L for the GOE. 



The NNS distribution describes the distribution of 
level spacings but does not contain information about 
their correlations. Such information is contained in an- 
other fluctuation measure, the A3-statistic. The number 
staircase function 



tf(E)= J dE' ^5(E'-E„) 

J —OO 



(20) 



counts the number of eigenvalues below energy E. With 
increasing energy, it increases by unity as E passes 
a (non-degenerate) eigenvalue and is otherwise con- 
stant. The number of eigenvalues in the energy interval 
[Eo, E + L] is given by n(E , L) = N{E Q +L)- N{E ). 
By definition of the mean level spacing d(E), we have 
n(E ,L) = L/d(E ). (We use that for N -> 00, d(E) is 
constant (independent of E) in any energy interval con- 
taining a finite number of levels). The number variance 
S 2 (L) = n 2 (Eo,L)— (n(E 0l L)) 2 is a fluctuation mea- 
sure which contains information about correlations be- 
tween level spacings. Suppose, for instance, that actual 
GOE spectra can be constructed by drawing spacings at 
random from the NNS distribution. In this case, Yig{L) 



3. Unfolding of Spectra. Purity and Completeness. 

In the limit N — ► 00, the average level density p(E) of 
the GOE in Eq. (fl"4"|) is constant in every energy interval 
containing a finite number of levels, and the same is true 
of the average level spacing d. In nuclei the situation dif- 
fers: The level density grows nearly exponentially with 
energy. In many cases, even a fairly short stretch of levels 
displays this fact: The spacings of the lowest-lying lev- 
els are consistently larger than those of the highest-lying 
ones. That fact distorts the spectral fluctuation mea- 
sures and must be taken into account prior to comparing 
data with GOE predictions. This is done by "unfolding" 
the spectra: The actual spectrum is modified such that 
the average level spacing is constant. GOE predictions 
relate to spectra consisting of levels with identical quan- 
tum numbers. Spectra obtained experimentally may be 
incomplete (i.e., miss levels (especially those with small 
or very large widths)), or not be pure (i.e., may con- 
tain levels with uncertain or incorrect quantum number 
assignments). It is important to know how lack of com- 
pleteness and/or purity affects the comparison of data 
with GOE predictions. 



w 




FIG. 5 The number variance versus the length L of the in- 
terval (L is in units of the mean level spacing) for the three 
canonical ensembles. Top curve: GOE; middle curve: GUE; 
bottom curve: GSE. The parameter /3 is the Dyson index, see 
Fig.H From Guhr et al. (1998). 



eigenvalues are again given by the values of the fitted 
staircase function taken at the original eigenvalues E^. 

How does the omission of levels affect spectral fluctu- 
ation properties? If a fraction / of levels is removed at 
random from a complete sequence (no missing levels, no 
levels with wrong quantum numbers) then the resulting 
correlation functions can easily be related to the corre- 
lation functions of the complete sequence (Bohigas and 
Pato, 2004). As a consequence, in the limit of large level 
numbers the information on spectral correlations is fully 
preserved provided / is known. In nuclei, levels with 
small widths are hard to detect and easily missed. But 
widths are theoretically predicted to be uncorrelated with 
the positions of eigenvalues. The removal of levels with 
small widths is thus a random process, and the analysis 
of Bohigas and Pato (2004) applies. The fraction / can 
be estimated using the Porter-Thomas distribution and 
the experimental detection efficiency. The admixture to 
a complete sequence of levels of levels with wrong quan- 
tum numbers can be treated similarly (Bohigas and Pato, 
2004). 



E. Discussion 



Unfolding requires the knowledge of the average level 
density p(E) for the data at hand. The situation is easy 
if a theoretical prediction for the average level density 
is available. This is the case, for instance, in billiards 
(where a point particle moving in two dimensions is scat- 
tered elastically on some surface) . Here the Weyl formula 
(see Baltes and Hilf, 1976) gives the average level den- 
sity in closed form in terms of the area enclosed by the 
surface and the length of the boundaries of that surface. 
Given p(E), the spectrum (or the spectra) are subse- 
quently unfolded by mapping the eigenvalues E^ onto 
new eigenvalues e M by the prescription 



AE p(E) 



(23) 



By construction, the new eigenvalues are dimensionless 
and have an average level spacing equal to unity. The 
En can be used to construct the NNS distribution and 
the A 3 -statistic. We observe that the right-hand side of 
Eq. (|2"3"|) is the average of the staircase function defined 
in Eq. (f20|) . The unfolded eigenvalues e M are the values 
of that function taken at E^. Usually, however, the ex- 
act form of the average level density is not known. If the 
data are obtained by numerical simulation of an ensem- 
ble (diagonalization of many matrices), the average level 
density is best found by numerically averaging over the 
ensemble. If we deal with an empirical spectrum of, say, 
several tens of levels, it is advantegeous to use the data 
to construct the staircase function rather than the level 
density (the representation of the latter in the form of a 
histogram depends on the bin width chosen), and to fit 
a low-order polynomial to that function. The unfolded 



The random-matrix approach described above is based 
on a few general principles: Invariance of the system un- 
der time-reversal (which leads to an ensemble of real 
and symmetric matrices), absence of a preferred direc- 
tion in Hilbert space (which makes the ensemble generic 
and implies orthogonal invariance), and confinement of 
the spectrum to a singly-connected finite stretch of the 
energy-axis (which is realized for the Gaussian cutoff as 
well as for many other cutoff factors). 

The GOE is universal: Cutoff factors different from 
the Gaussian cutoff but subject to the conditions formu- 
lated in Section III. C. 21 lead to different forms of the av- 
erage spectrum but to identical predictions for the local 
spectral fluctuation measures. GOE predictions of local 
spectral fluctuation measures are useful for the analysis 
of data because the ensemble is ergodic: For almost all 
members of the ensemble, the fluctuation measure calcu- 
lated as an average over the ensemble equals the result 
obtained by taking a running average over the spectrum 
of that member. 

It is rather amazing that the few principles just men- 
tioned lead to parameter-free quantitative predictions for 
the local spectral fluctuation measures. Eq. (fT5|) gives the 
distribution law relevant for transition matrix elements 
and decay widths. Eqs. (fl9|) and ([22)) present the two 
measures which are most frequently used for the analy- 
sis of spectral data. Needless to say, other fluctuation 
measures have also been worked out (see Brody et al., 
1981). 

Random-matrix theory predicts fluctuation properties 
in terms of mean values. In the examples treated so far, 
the input mean value has been the average of the squares 
of the projected wave functions as in Eq. (fT8"|) , or the local 



average level density p(E) (or the average level spacing 
d) as in Eqs. ([!§)) and P2j) . Both parameters must be 
determined from the data. Then the theory predicts the 
distribution of the relevant observables. In case the level 
density changes significantly over the length of the given 
spectrum, an unfolding of the spectrum must precede the 
comparison with GOE predictions. 



F. Random Matrices and Chaos 
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tion of the spectral fluctuation measures, their numerical 
accuracy had to be much better than the average level 
spacing. It was also important to make sure that no 
eigenvalue was missed. The Weyl formula (see Baltcs 
and Hilf, 1976) was used to unfold the spectrum. The 
number of eigenvalues was large enough to evaluate for 
the first time the A3-statistic (in addition to the NNS 
distribution which had been used before). The results 
are shown in Figs. [6] and [7] 



There exists a close connection between random matri- 
ces and quantum chaos. The latter term refers to quan- 
tum systems that are chaotic in the classical limit. In 
classical mechanics, chaos is a dynamical property char- 
acterized by the exponential divergence in time of trajec- 
tories starting in close-lying points of phase space. The 
phase space of a fully chaotic system is filled with such 
chaotic trajectories and is void of islands with regular 
dynamics. The analysis of chaos in classical conservative 
systems uses the existence of periodic orbits and was pi- 
oneered, among others, by Gutzwiller. The results are 
summarized in his book (Gutzwiller, 1995). Periodic- 
orbit theory has been successfully applied to many sys- 
tems with few degrees of freedom (including systems that 
are not fully chaotic). To the best of our knowledge, 
a similarly complete understanding of chaos in classical 
many-body systems with their high-dimensional phase 
space does not exist. Atomic nuclei pose the additional 
difficulty that the matter density is very high. Even in 
the classical limit it is not possible to neglect the fact 
that neutrons and protons are fermions. This fact leads 
to complications regarding periodic-orbit theory (see, for 
instance, Sommermann and Weidenmuller, 1993; Wei- 
denmiiller, 1993; Sakhr and Whelan, 2003). 

Since the late 1970s much effort has been devoted to 
identifying the dynamical properties of quantum systems 
that are fully chaotic in the classical limit (see Gutzwiller, 
1990; Haake, 2001). There were two lines of develop- 
ment: Some authors looked for signals of classical chaos 
in the time-evolution of wave packets, others focused at- 
tention on the fluctuation properties of eigenvalues and 
eigenfunctions of closed systems (which possess a dis- 
crete spectrum). Here we confine ourselves to the lat- 
ter development which has had repercussions in nuclear 
physics. Mounting numerical evidence from classically 
chaotic few-degrees-of-freedom systems due to McDon- 
ald and Kaufman (1979), Casati (1980), Berry (1981), 
and others culminated in the work of the Orsay group 
(Bohigas et at, 1984) on the Sinai billiard: A point par- 
ticle moving in two dimensions is scattered elastically by 
the interior surface of a square and by the exterior sur- 
face of a circle inscribed into the square. That system is 
fully chaotic and invariant under time reversal. Solving 
the Schrodinger equation for the Sinai billiard numeri- 
cally, the Orsay group accumulated a sequence of about 
1000 consecutive eigenvalues belonging to eigenfunctions 
of the same symmetry class. For a meaningful evalua- 
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FIG. 6 The A3-statistic for the Sinai billiard (open cir- 
cles), the GOE prediction (solid line), and the Poisson result 
(dashed line). From Bohigas et al. (1984). 

In the figures the data on the Sinai billiard are com- 
pared with GOE predictions and with the Poisson dis- 
tribution. The latter has an exponential form, is typ- 
ical for regular (or integrable) systems (Berry and Ta- 
bor, 1977), and is briefly explained below. The fig- 
ures show excellent agreement between the results for 
the Sinai billiard and the GOE predictions. This agree- 
ment led Bohigas et al. (1984) to formulate a con- 
jecture (the "Bohigas-Giannoni-Schmit (BGS) conjec- 
ture"): The spectral fluctuation properties of a quan- 
tum system which is fully chaotic in the classical limit 
coincide with those of the canonical random-matrix en- 
semble having the same symmetry. Massive numerical 
evidence on other few-degrees-of-freedom-systems has 
given extensive support to the conjecture. Attempts to 
prove the conjecture analytically have partly been based 
on periodic-orbit theory (which has been very impor- 
tant all along in the understanding of classical chaos) 
and the semiclassical approximation. Berry (1985) has 
established a connection between classical chaos and the 
A3 statistic. More recently and for a special fluctua- 
tion measure (the two-point function, a quantity inti- 
mately related to the A3-statistic), the conjecture has 
been proved or, to use a mathematically less demanding 
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corresponds to the existence of a complete set of quan- 
tum numbers that label every state. There is no level 
repulsion, and no correlation between levels. As shown 
by Berry and Tabor (1977), this case generically yields 
a spectrum where the spacings have an exponential dis- 
tribution ("Poisson spectrum"), see Fig. [7] The general 
case in classical mechanics is the one where phase space 
consists of islands of regular trajectories separated by 
domains filled with chaotic trajectories. There does not 
seem to exist a generic description of the distribution of 
level spacings for that case. The Brody distribution de- 
fined in Eq. (|2"8)) below is a purely heuristic (and not the 
only) interpolation formula between the Poisson distri- 
bution and the Wigner surmise. 



G. Doorway States 



FIG. 7 The NNS distribution for the Sinai billiard (his- 
togram), the GOE prediction (solid line), and the Poisson 
result (dashed line). While the theoretical literature com- 
monly uses the label s for the actual level spacing in units of 
the mean level spacing, that quantity is usually referred to as 
x when RMT predictions are compared with data. We follow 
that usage here. Insert: The eigenvalues must be selected 
according to the symmetries of the eigenf unctions. From Bo- 
higas et al. (1984). 



term, demonstrated in Hcusler et al. (2007). The argu- 
ment uses generic properties of periodic orbits in classi- 
cally chaotic systems with few degrees of freedom. 

The semiclassical arguments used in Heusler et al. 
(2007) do not apply directly to chaos in classical or quan- 
tum many-body systems. This point was made at the 
beginning of this Section and is taken up again in Sec- 
tion IIII. C. 11 For the time being, we adopt the BGS 
conjecture also for many-body systems. Thus we speak 
of quantum chaos (or, briefly, of chaos) whenever in a 
many-body system such as the nucleus the spectral fluc- 
tuation measures introduced above agree with GOE pre- 
dictions. We have to keep in mind that the NNS distri- 
bution is a less safe indicator of quantum chaos than the 
A 3 -statistic. To see this, we consider (Rosenzweig and 
Porter (I960)) a real symmetric matrix with diagonal en- 
tries that are eigenvalues of a regular system with average 
level spacing d and with non-diagonal elements of typical 
strength v (see the example in Eq. ([5])). The parameter 
which rules the regular-to-chaos transition is v/d. As we 
increase that parameter from very small values, neighbor- 
ing levels begin to repel, and the NNS distribution of the 
GOE is approached. This happens before the long-range 
stiffness of the spectrum as manifest in the A3-statistic 
is attained. 

In classical mechanics, the case of complete chaos is a 
limiting case. Another limiting case is that of completely 
regular motion. In quantum mechanics, regular motion 



As explained in Section Hl.Bl above, the GOE yields the 
generic description of the spectra of bound quantum sys- 
tems. It offers the best first guess of spectral properties 
if we have no specific knowledge of the system except for 
the fact that it is invariant under time reversal. There 
are cases, however, where we possess some limited addi- 
tional dynamical information. We may then look for a 
theoretical description which takes that information fully 
into account but is otherwise generic. Doorway states are 
a case in point. Other examples for this type of approach 
are described in Section UlI.DI below. 

We consider a particular mode of excitation of the sys- 
tem. To be specific, we take the electric dipole operator 
D acting on the ground state \g) of an even-even nucleus 
in the long- wavelength limit. With H the Hamiltonian, 
we shift for simplicity the energy such that the ground 
state has energy zero, H\g) = 0. We choose the normal- 
ization of the dipole mode |0) = D\g) such that (0|0) = 1. 
Then the expectation value Eq — (0\H |0) gives the mean 
excitation energy of the dipole mode. In general the 
dipole mode is not an eigenstate of H . Therefore the 
variance of H with respect to the dipole mode, i.e., the 
expression AH 2 = (0\H 2 \0) — JSg = E^if^) 2 , does 
not vanish. (Here y, labels states with the same quan- 
tum numbers as, but orthogonal to the dipole mode). 
As a consequence the cross section for dipole absorp- 
tion possesses a large number of sharp lines, each oc- 
curring at an eigenstate of H. (Here we disregard the 
fact that for most nuclei, Eq is greater than the thresh- 
old for particle emission, so that the sharp lines actu- 
ally become more or less broad resonances). We ob- 
serve that (1/N)AH 2 = (1/A0]T; m >i(#Om) 2 represents 
the mean coupling strength of the dipole mode with the 
other states of the system. The ratio of that expression 
and the mean level spacing (taken at Eq) is a measure 
of the length of the energy interval over which the dipole 
mode is strongly mixed with other states. Within that 
interval, the peak heights of the dipole absorption lines 
are enhanced, and an average of the cross section for 
dipole excitation (taken with a Lorentzian weight func- 
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tion whose width is larger than the mean level spacing) 
displays a resonance at or near the energy E . This is 
the giant dipole resonance. Dipole absorption may be 
viewed as a two-step process: First the dipole mode is 
formed. The ground states of even-even nuclei have spin 
zero and positive parity and the dipole operator is the 
z-component of a vector. Therefore in these nuclei the 
dipole mode has spin one and negative parity. Subse- 
quently that mode decays into the eigenstates of H with 
the same quantum numbers. It is obvious that that same 
picture applies to many other nuclear reactions: neutrons 
impinging on a nucleus give rise to a single-particle mode, 
protons generate isobaric analog modes, etc. The picture 
can likewise be used to describe the distribution of sim- 
ple configurations like one-particle one-hole states over 
the eigenstates of the system. Early summaries of the 
doorway state idea may be found in Bohr and Mottel- 
son (1969) and Mahaux and Weidenmuller (1969). For a 
recent discussion which goes beyond the standard model 
described in the next paragraph, see Zelevinsky et al. 
(1996) and De Pace et al. (2007). 

We turn to the standard description of a doorway state 
within random-matrix theory. To this end, we define 
the following extension of the GOE. With Eq the mean 
excitation energy of the doorway state, and with Ho^ 
(where /i = 1 , . . . , N) the coupling matrix elements of 
the doorway state with the other states of the system, 
the Hamiltonian matrix H for the doorway-state model 
takes the form 



H 






(24) 



Here H^ is a GOE matrix of dimension N. With |0) 
the doorway state, we obviously have (0|J?|0) = E and 
(1/N)AH 2 = (I/AOEmC^V) 2 - Because of the distinct 
role of the doorway state, the doorway model of Eq. (124]) 
is not orthogonally invariant. It does possess that invari- 
ance, however, in the subspace of states carrying labels 
fj, > 1. That statement implies that ensemble averages 
of observables cannot depend on the individual coupling 
matrix elements Hq^ but depend only on the orthogo- 
nal invariant (1/iV) E u (^0m) 2 j the mean square coupling 
matrix element. 

Every random-matrix model that incorporates addi- 
tional dynamical information does so at the expense of 
complete orthogonal invariance. That invariance is ex- 
tremely helpful in working out the spectral properties of 
the GOE. Therefore non-invariant extensions of the GOE 
are often very difficult to handle. This is not the case for 
the model of Eq. (|2"4")1 because we add only a single state 
to the GOE. As a consequence, in the limit N — > oo 
both the spectral statistics and the Porter-Thomas dis- 
tribution of the Hamiltonian (j2~4"|) coincide with those of 
the GOE. The only distinct feature of the model is the 
strength function for the doorway state. It is defined 
as the ensemble average of E T I v-H 7 ")! 2 ^^ — E T ) and 
gives the probability per unit energy interval to find the 
doorway state admixed to the eigenstates \t) of H with 



eigenvalues E T . Put differently, the strength function 
is the average value of the level density weighted with 
the square of the overlap matrix element. In form, the 
strength function is closely related to the local density 
of states used in condensed-matter physics. Because of 
the ensemble average, the strength function is a smooth 
function of energy E and has the form 



J2mr)ME-E T ) = l 



r' 



2vr(£-£ )2 + (l/4)(ri)2 



Here 



r 1 = 2tt 
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0/-J 



P(E) 



(25) 



(26) 



is the spreading width of the doorway state, with p(E) = 
(l/d) the average level spacing of the GOE at energy E. 
Since (l/N)AH 2 /d = T l /2ir, the spreading width mea- 
sures the length of the energy interval within which the 
eigenstates of H carry significant admixtures of the door- 
way state. The notation for the spreading width with a 
downarrow is a reminder of the fact that the spreading 
width does not account for a decay process into some 
open channel (with the ensuing probability flux of parti- 
cles at large distance) , but for the mixing of a particular 
mode with other bound states. 

The Lorentzian form of the strength function applies 
approximately for Y* <C A where 2A is the radius of the 
semicircle. For small coupling, Eq. (|26p has the form 
of Fermi's Golden Rule. In the framework of random- 
matrix theory, the result ([26| is valid beyond the pertur- 
bative regime, however, and holds even if T^p(E) ^> 1. 
Integrating the left-hand side of Eq. (f25|) and using com- 
pleteness we obtain unity. That same statement ap- 
plies to the right-hand side, and the strength function 
is properly normalized. The Lorentzian form (|25[) does 
not apply when the energy Eo of the doorway state is 
close to one of the end points of the semicircle, or when 
the spreading width becomes very large (i.e., comparable 
with the radius of the semicircle), see Kota (2001), De 
Pace et al. (2007), and numerical examples given, for 
instance, in Zelevinsky et al. (1996). 

The spreading width T^ has a remarkable property 
which makes it a useful measure for the spreading of a 
doorway state. We ask: How does the mean square ma- 
trix element (1/^) Eu(%) 2 change with the average 
level density of the GOE states? That question arises 
in the nuclear context because doorway phenomena are 
encountered at various excitation energies and in nuclei 
with widely different mass numbers for which the nu- 
clear level density differs markedly. An intuitive answer 
is obtained by noting that a significant increase of the 
level density implies a significant increase in the com- 
plexity of the wave functions making up the GOE states. 
Therefore, each of the matrix elements Hq^ connecting 
the doorway state with the GOE states is strongly re- 
duced, and so is the mean square matrix element. But in 
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the expression (j26f for the spreading width this strong re- 
duction of the mean square matrix element is essentially 
compensated by the increase of the average level density. 
That compensation is exact in simple models (Brody et 
al, 1981) and is expected to hold to a good degree of 
approximation in realistic cases. An experimental veri- 
fication of this expectation comes from investigations of 
the spreading width for isospin mixing (Harney et al, 
1986) and will be discussed in Part 2 of this review. We 
conclude that in contrast to the exponential dependence 
of the average level density on excitation energy and mass 
number, the spreading width is expected to be a slowly 
varying function of these parameters and is thus a useful 
measure for the spreading of a doorway state. 

If the model (|24|) applies in reality, a doorway state 
has hardly any influence on spectral properties of the 
system: The average level density of the states that 
carry the same quantum numbers as the doorway state is 
unchanged, their spectral fluctuation properties are the 
same as for the GOE, and their partial widths for de- 
cays different from electric dipole decay to the ground 
state have the same Porter-Thomas distribution as for 
the GOE. The only difference to the pure GOE case is 
the Lorentzian enhancement ([25]) of the strength function 
for dipole absorption. This is the only trace left of the 
doorway state after we take account of its mixing with 
the complicated states. Dividing the partial widths for 
dipole absorption of the eigenstates of H by the value of 
the strength function (taken at the corresponding eigen- 
value) removes that trace and should yield quantities that 
have a pure Porter-Thomas distribution. 



III. APPLICATION OF RMT TO NUCLEAR SPECTRA 
A. General Remarks 

Nuclear energy levels are characterized by quantum 
numbers that reflect the symmetries of the nuclear 
Hamiltonian: Total spin (J) reflects rotational symme- 
try, parity (II) reflects invariance under mirror reflection, 
isospin (T) reflects proton-neutron symmetry. We ex- 
hibit the consequences of such symmetries for the appli- 
cation of RMT to nuclear spectra. 

The total Hilbert space is spanned by many-body wave 
functions that carry the quantum numbers J, II, T. These 
can be arranged in such a way that the matrix represen- 
tation of the nuclear Hamiltonian has block structure, 



H 








#J 2 n 2 T 2 



(27) 



Here {J^^} ^ {J 2 n 2 T 2 } ^ . . ., and each of the ma- 
trices H JUT couples only many-body states which carry 
the same quantum numbers. If in addition the nuclear 
dynamics is chaotic, then the BGS conjecture (see Sec- 
tion III.Fp implies that each of the matrices H JI1T is 



obey time-reversal invariance, the suitable ensemble is 
the GOE. In the framework of RMT, Hamiltonian ma- 
trices referring to different sets of quantum numbers are 
assumed to be uncorrelated. 

To compare RMT predictions on spectral fluctuations 
with data on nuclear energy levels, sequences of lev- 
els carrying the same quantum numbers are needed. 
The data are subject to three requirements: (i) The 
sequence(s) should be as long as possible, (ii) the se- 
quence^) should be pure (i.e., should not contain levels 
carrying quantum numbers which differ from those of the 
rest), and (iii) the sequence(s) should be complete (i.e., 
there should not be any levels that were not detected). 
The first requirement is needed to ensure that the run- 
ning average over the actual spectrum is as close as pos- 
sible to the running average over the complete spectrum, 
the latter by ergodicity (Section III.C.3[) being equal to 
the GOE ensemble average. The two other requirements 
guarantee that the statistical predictions of RMT can 
meaningfully be applied to the data, see Section III. D. 31 
and the discussion in Section IIII.DI Unfortunately, the 
number of nuclear data sets of sufficient quality to pro- 
vide detailed tests of RMT is fairly limited. This is pri- 
marily due to the requirements of purity and complete- 
ness imposed by the sensitivity of the standard fluctua- 
tion measures. 

In many cases nuclear levels are observed as narrow 
particle-unstable resonances, see Fig. [TJ Then a multi- 
level i?-matrix fit (Lane and Thomas, 1958) is used to 
determine the positions the levels would have if they were 
stable under particle decay. These positions are used to 
calculate level spacings and to test GOE predictions. (In 
R-matrix theory, the nucleus is thought to be enclosed 
by a ficticious boundary that lies some distance beyond 
the nuclear radius. Boundary conditions on that surface 
and the nuclear Hamiltonian jointly define a set of dis- 
crete states within the boundary. These states appear as 
resonances in the scattering matrix S. Approximations 
to the resulting formal expression for S serve as the basis 
of the fits to data). 

In comparing nuclear data with GOE predictions for 
the NNS distribution and/or the A3-statistic, one faces 
a difficulty: Both distributions are parameter free, and 
it is difficult to assess the significance of the usual tests 
for goodness of fit such as the % 2- test when one is far 
from these limiting cases. Therefore, one uses measures 
which interpolate between the GOE prediction and the 
case of a totally regular system. These do have free pa- 
rameters, and the goodness-of-fit tests are easily inter- 
preted. The NNS distribution has the form of the Wigner 
surmise (fT9|) for the GOE and is proportional to exp(— s) 
(Poisson distribution) for regular systems (Berry and Ta- 
bor, 1977). An expression that interpolates between both 
is the Brody distribution (Brody et al. (1981). It depends 
on a single parameter uj and is given by 



Puj(s) = (1 + uj)as u exp(-a.s 1+w ) 



(28) 



a member of a random-matrix ensemble. Since nuclei where s is the actual level spacing in units of the mean 
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level spacing and the constant a = [r{(2 + u>)/(l + 
co)}] 1+UJ is fixed by normalization. For u = (ui = 1), 
the Brody distribution equals the Poisson distribution 
(the Wigner distribution), respectively. For all a; > 0, 
the Brody distribution vanishes at s = 0. The Brody 
formula is only one of several formulas that interpolate 
between the Poisson and the Wigner distribution. An- 
other example is the Berry-Robnik distribution (Berry 
and Robnik, 1984) . 

For the GOE, the A3-statistic has the logarithmic de- 
pendence on the length L of the energy interval shown 
in Eq. (|2"2")1 , while for a regular system it is linear in L. 
A parameter-dependent measure for deviations from the 
GOE is obtained by considering a spectrum which is a 
superposition of k independent GOE spectra. For k 3> 1 , 
the A3-statistic approaches the linear dependence of the 
regular case. (Apparently this was first noticed by Gure- 
vich and Pevsner (1956)). The deviation from the GOE 
prediction is significant already for k = 2, see Fig. [TBI 
below. 

The predictions of RMT on fluctuation properties of 
nuclear wave functions can only be tested in terms of 
the distribution of matrix elements (either for decay into 
open channels or for electromagnetic transitions, weak 
interaction matrix elements not being numerous enough 
for such a test). Here the size of the sample again is 
important. In addition there is usually an experimental 
cutoff for small matrix elements so that only part of the 
Porter-Thomas distribution can be tested. 

Two properties of nuclei are central for tests of RMT. 
(i) In every nucleus, the average level density increases 
roughly exponentially with excitation energy. Thus, 
while typical level spacings near the ground state are 
several hundred keV, spacings of levels having the same 
spin and parity (a subset of all levels!) at neutron thresh- 
old in heavy nuclei are typically 10 eV, see Fig. [TJ For 
fixed excitation energy the level density increases like- 
wise with mass number A (save for corrections due to 
nuclear shell structure, see Section [iV.A.l[) . The require- 
ments on experimental energy resolution obviously in- 
crease with increasing level density and, in general, limit 
nuclear spectroscopy except for very fortuitous situations 
such as those leading to the data in Fig. [TJ (ii) Lev- 
els below the threshold for particle emission have only 
small widths (in comparison with the mean level spacing) 
due to beta- or gamma-decay. Above particle thresh- 
old, the total widths of the nuclear resonances increase 
rapidly with increasing excitation energy. This is because 
the number of open channels for particle decay increases 
rapidly (the number of states available for decay in the 
daughter nuclei increases roughly exponentially with ex- 
citation energy in these nuclei). As a consequence, iso- 
lated resonances as shown in Fig.[JJare observed only just 
above the lowest particle threshold. A few hundred keV 
above that threshold, resonances begin to overlap (the 
mean level spacing decreases, the average total width in- 
creases), and it is no longer possible to investigate spec- 
tral fluctuations. Rather, this is the domain of statistical 



nuclear reaction theory (see Part 2 of this review). Thus 
tests of GOE predictions in nuclear spectra are limited 
to the energy interval between the ground state and an 
energy somewhat above the first particle threshold. 

In describing the application of RMT to nuclear data, 
we first discuss the experimental methods that have been 
used to obtain the relevant spectral information (Sec- 
tion IIII.Bp . We then review the results of comparing 
the data with GOE predictions on spectral fluctuations 
( Section IIII. C[) . RMT can be extended to deal with vi- 
olations of symmetry or invariance. This is described in 
Section HTLDl 



B. Experimental Methods 

1. Neutron Resonances 

The early tests of RMT involved neutron resonances, 
as shown in Fig. [TJ The scattering of slow neutrons en- 
ables the study of individual resonances in a narrow win- 
dow of energies at high excitation energy in the com- 
pound nucleus - typically 5 to 7 MeV. At these energies 
the level density in medium-weight and heavy nuclei is 
very large. However, the angular momentum barrier for 
the incident neutron severely restricts the neutron's or- 
bital angular momentum and thus the spins J of the com- 
pound nuclear resonances which contribute to the scat- 
tering. For slow neutrons, only s- and p-wave resonances 
are usually observed. 

The experimental method of choice is a time-of-flight 
measurement. Longer flight paths allow for better en- 
ergy resolution (essential to resolve the resonances), but 
reduce the counting rate because of the smaller detec- 
tor solid angle. Thus high-intensity neutron sources are 
required. Today the most intense neutron beams are pro- 
duced at spallation sources. 

The most common experiment is a transmission mea- 
surement. The transmission of the neutron beam through 
a target with nuclei of mass number A determines the to- 
tal cross section for the n + A reaction. Neutron capture 
followed by 7 emission is also very helpful in determining 
the resonance parameters. Analysis of the resonance data 
is normally performed with the Lane and Thomas version 
of the Wigner-Eisenbud R-matrix formalism (Lane and 
Thomas, 1958). The classic monograph on neutron reso- 
nance reactions is by Lynn (1968). 

For a comparison with RMT predictions the levels in 
a sequence must have the same quantum numbers. Thus 
one key issue is to determine the spin J and parity n of 
each resonance. This is done using the angular momen- 
tum £ of the scattered neutron. For spin-zero targets all 
s-wave resonances have J n = l/2 + . For slow neutrons, 
the difference in penetrabilities for £ = and £ — 1 is so 
large that at first sight £ can be assigned by inspection 
- strong resonances are s-wave and weak resonances are 
p-wave. One normally formalizes this with a Bayesian 
analysis (Bollinger and Thomas, 1968), but this approach 
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is not reliable in the gray area between weak s-wave and 
strong p-wave resonances. Of course other experiments 
can be used to improve the spin and parity assignments. 
We mention neutron capture with high resolution 7-ray 
spectroscopy or with calorimeters (where for each cap- 
ture reaction the total number and individual energies of 
the emitted gamma rays are registered). However, these 
are very time consuming. 

In addition to the issue of spurious resonances (incor- 
rect spin or parity assignments), the other major prob- 
lem is missing levels. The missing levels are expected to 
be the weakest levels, hence the focus on signal-to-noise 
ratios and energy resolution in resonance measurements. 
The Porter-Thomas distribution predicts many weak res- 
onances. Assuming that distribution one can estimate 
the fraction of missing levels. Unfortunately this nearly 
universally used correction method can be misleading if 
non-statistical effects are present. An alternative ap- 
proach utilizing the NNS distribution was developed only 
recently (Agvaanluvsan et al., 2003). Bohigas and Pato 
(2004) extended the investigation of the effects of missing 
levels to other level fluctuation measures. 

Almost all of the neutron resonance data suitable for 
detailed comparison with RMT are obtained for spin- 
zero targets. Most of the neutron resonance data used in 
the early evaluation of RMT were obtained by Rainwa- 
ter's group at Columbia (Liou et al, 1972a and 1972b) 
for a number of nuclei with mass numbers A > 110. Due 
to experimental limitations, the number of resonances in 
each nucleus was never significantly larger than 200. 



high-quality proton resonance data are also for s-wave 
sequences, there are (in contrast to the neutron case) a 
few p-wave sequences that are considered pure and com- 
plete. With very few exceptions the level density becomes 
too great for this method to work much beyond mass 
number A = 60. Here the typical number of resonances 
of the same spin and parity is 50 or so. Thus the proton 
resonance data complement the neutron resonance data; 
the best results for each set are obtained in quite differ- 
ent mass regions. Almost all of the data used to compare 
with RMT are for spin-zero targets and from the Trian- 
gle Universities Nuclear Laboratory (TUNL) (Wilson et 
al, 1975; Watson et al., 1981). 

In order to observe (nearly) all of the resonances one 
needs very good beam-energy resolution and high beam 
intensity. These requirements seem contradictory. One 
approach is to accept the time-dependent energy fluctu- 
ations intrinsic to the accelerator and make a correction 
later. The most successful correction method solves the 
resolution-intensity impasse by using two beams. One 
(high-intensity) beam is used to perform the experiment: 
the other beam is used to generate a feedback signal that 
follows the beam energy fluctuations; this signal gener- 
ates a voltage difference which is applied to the target. 
Thus the time dependent energy fluctuations are can- 
celed. The method works well for a Van de Graaff accel- 
erator where most of the fluctuations have low frequency. 
The details are given by Bilpuch et al. (1976). 



2. Proton Resonances 



Due to the Coulomb barrier proton resonances cannot 
be studied near zero bombarding energy as neutron res- 
onances are. High-resolution proton resonance measure- 
ments are typically taken at bombarding energies corre- 
sponding to 60-70% of the Coulomb barrier. This has 
two major advantages. First, the Coulomb barrier pen- 
etrability serves to narrow the proton widths and makes 
possible the resolution of individual resonances in rather 
dense spectra. Second, the addition of Coulomb and nu- 
clear resonance scattering amplitudes leads to striking 
interference patterns that are used to identify the spin 
and parity of each resonance. The parity assignment is 
normally apparent by inspection, since the interference 
patterns for even and odd orbital angular momenta are 
quite different. This is important because the proton 
penetration factors for different orbital angular momenta 
do not differ as much as for the neutron resonances. As 
a result one usually observes s~, p-, and <i-wave and 
sometimes even /- and g-wave resonances. The primary 
difficulty consists in determining the J value of the pro- 
ton resonances. Additional experiments (inelastic scat- 
tering or capture) can resolve this problem, but as in the 
neutron case, these experiments are very time consuming 
and therefore are rarely performed. Although most of the 



3. Low-lying Levels 



The spectroscopy of low-lying levels (excitation ener- 
gies below 2 MeV or so) has always been a primary object 
of study in nuclear physics. Many different approaches 
have been used: various nuclear reactions including in- 
clastic scattering, pickup, and transfer reactions, 7-ray 
spectroscopy following (3 or a decay, etc. Almost all of 
these processes are quite selective. Therefore, one needs 
to use many different approaches to ensure that all levels 
(in some energy interval) are observed. One very pow- 
erful technique uses the neutron capture reaction. Neu- 
tron capture on nucleus A is followed by sequences of 7 
transitions which finally populate the ground state of nu- 
cleus A + l. The average neutron resonance capture tech- 
nique (Bollinger and Thomas, 1968) effectively averages 
over many neutron resonances and is non-selective; it 
also averages over Porter Thomas fluctuations, increas- 
ing the probability of observing weak transitions. Every 
low-lying state within some spin range is expected to 
be populated. The combination of neutron capture and 
direct reactions has led to a number of complete level 
schemes at low energies (von Egidy et al., 1986). 
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4. High-spin States 

The ground states of even-even nuclei have spin zero, 
those of even-odd, odd-even and odd-odd nuclei have 
small spin values. The excitation energy E(J) of the 
lowest state with given spin J generically increases with 
J; the function E(J) defines the "yrast line", see Fig. [8] 
High-spin states located near the yrast line are typically 
investigated via the collision of two heavy nuclei ( "heavy- 
ion collisions"). At non-zero impact parameter, the two 
colliding nuclei with mass numbers A\ and Ai typically 
carry a large angular momentum of relative motion. The 
high-spin intermediate complex formed by the collision 
(with spin values as large as 60?i or so) may decay by a 
sequence of 7 transitions (perhaps with intermittent neu- 
tron evaporation) to the ground state of a nucleus whose 
mass number is smaller than but close to A\ + A2. The 
7 rays emitted in this process are analyzed using large- 
scale gamma-ray detection arrays. This approach has 
generated a very large amount of data. These comprise 
many rotational bands (with many states in each) for a 
range of heavy nuclei. (Rotational bands are typical for 
deformed nuclei and are dealt with in Section HV.B.ip . 
The observed states have relatively high spin and rather 
large excitation energies, but are not far above the yrast 
line. Some aspects of the method of analysis are summa- 
rized by D0ssing et al. (1996), where further references 
may be found. 

Unfortunately, the difficulties in obtaining suitable 
data sets for comparison with RMT are many. One of- 
ten sees many rotational bands, but of course each band 
has only one state of a specific spin and parity. Thus 
one is forced to combine results from many bands and 
nuclei. Another serious issue is the problem of quan- 
tum number assignments. Within a given band with a 
well-known bandhead the assignments are reliable; as- 
signments based on interband transitions are more prob- 
lematic. Until now, levels up to a few 100 keV above 
the yrast line have been analyzed (Garrett et ai, 1997); 
higher-lying rotational bands cannot be individually re- 
solved. The evidence here points to regular motion, see 
Section IIII.CJ.3l Theoretical expectations are that at 
around 800 keV above the yrast line, the spectral fluc- 
tuations become chaotic, see Section liV.B. 21 It is to be 
hoped that with improved resolution (perhaps attainable 
with the next generation of large-scale detectors), spec- 
troscopic data in that interesting energy region will be- 
come available. 



5. Complete Level Schemes 

The ideal is a complete scheme which begins at the 
ground state and extends into the neutron or proton reso- 
nance region, with perfect quantum number assignments 
to each level. Obtaining such a complete level scheme is 
exceptionally difficult at best, and impossible in medium- 
weight to heavy nuclei. The level densities are simply too 



great. For very light nuclei, on the other hand, the level 
density is small and the total number of states is not 
sufficient for a detailed statistical analysis. Only for nu- 
clei in the mass range between 20 and 40 or so does the 
level density have suitable values. These are essentially 
the nuclei belonging to the 2sld-shell, see Section HV.A.ll 
below. 

It might seem that complete spectra might best be 
measured by using a variety of reactions as done for the 
spectroscopy of low-lying states. However, this approach 
meets practical difficulties. Most reactions are not only 
selective, but also provide information only in a limited 
energy range. In the approach that was successfully used 
for two nuclei - 26 A1 and 30 P - the properties (quantum 
numbers, positions, widths) of a number of proton res- 
onances were determined. The proton capture reaction 
was then measured for these resonances which had dif- 
ferent quantum numbers. The method essentially guar- 
antees that all the levels below the proton separation 
energy are observed. Quantum numbers are assigned 
to the observed levels using high-resolution gamma-ray 
spectroscopy, including angular distributions of primary 
and secondary gamma rays. The experimental proce- 
dure for 30 P is described in detail by Grossmann et al. 
(2000). The general approach is summarized by Mitchell 
and Shriner (2001). With the help of the neutron cap- 
ture reaction, a "nearly complete" level scheme below 4.3 
MeV excitation energy was measured in 116 Sn (Raman 
et ai, 1991). 

The analysis of the complete spectra must allow for 
isospin-symmetry breaking and is dealt with in Sec- 
tion ITTlEIJ 



6. Low-lying Modes of Excitation 

The emphasis here is not on complete level schemes 
in a restricted energy range (as in Section IIII.B.3[) , but 
rather on phenomena that relate to levels of fixed spin 
and parity and seem linked to the concept of a door- 
way state, see Section III.G1 The primary example of a 
doorway state is the giant electric dipole resonance which 
manifests itself in a large resonance-like structure in the 
absorption cross section for 7 rays. The isobaric analog 
resonances provide another classic example of a doorway- 
state phenomenon. The detailed analysis of both types of 
resonances involves nuclear reaction theory and is dealt 
with in Part 2 of this review. 

There are several other interesting excitation modes at 
lower energy - including the low-lying isovector magnetic 
orbital dipole or scissors mode with J n = 1 + (Bohle et 
ai, 1984; Richter, 1995) and the electric pygmy dipole 
resonance with J n = l - (both spin assignments apply- 
ing to even-even nuclei) . Nuclear resonance fluorescence 
measurements have been used to generate extensive data 
sets of 1 + and 1~ states; the method provides a unique J 
value of 1 and a probable parity assignment, while mea- 
surements with a polarized photon beam provide a defini- 
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tive parity assignment. Work at Darmstadt has been fo- 
cused on both the scissors mode (Enders et al, 2000) and 
the electric pygmy dipole resonance (Enders et at, 2004). 



7. Summary 

Fig. [5] summarizes in a qualitative fashion the domains 
of excitation energy E x and mass number A where spec- 
tral fluctuations have been investigated. The figure is 
largely self-explanatory; suffice it to say that the proton 
resonances around A — 50 are measured above threshold 
but below the Coulomb barrier. The figure shows that 
high-spin states are measured at comparatively high ex- 
citation energies. This statement has to be taken with 
a grain of salt, however. For an even-even nucleus with 
mass number A around 160 or so, Fig. [5] shows schemat- 
ically the excitation energy of the lowest state with spin 
J versus J (the "yrast" line). The high-spin states that 
were analyzed so far lie up to a few 100 keV above the 
yrast line and, for that value of </, represent low-lying ex- 
cited states. Sequential decay of a rotational band close 
to the yrast line by repeated emission of 7 rays reduces 
both J and the overall excitation energy while the dis- 
tance to the yrast line remains essentially the same. 
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FIG. 9 The domains of excitation energy E x and mass num- 
ber A where spectral fluctuations have been investigated are 
shown for four classes of states (resonances, low-lying states, 
high-spin states, and states that belong to a complete spec- 
trum. Complete spectra are known for three nuclei only, 26 A1, 
30 P, and 116 Sn). The letter S denotes the particle threshold. 
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FIG. 8 The yrast line, i.e., the excitation energy E(J) of 
the lowest state with spin J versus J for an even-even nu- 
cleus with mass number around 160 (schematic) . The shaded 
area indicates the domain where spectroscopic information is 
available, see D0ssing et al. (1996). The letter S denotes the 
particle threshold. Depending on A, S typically lies between 
5 and 8 MeV. 



C. Tests of Fluctuation Measures 



The fluctuation measures described in Section III.DI 
have been applied to nuclear data obtained with the ex- 
perimental methods summarized in Section IIII.BI We 



emphasize again that due to the sensitivity of these mea- 
sures, the quality of the data sequence (the degree of 
purity and completeness) is of paramount importance; 
only a very small fraction of all nuclear data can be used 
for such tests. 

We mention in passing another measure, the "correla- 
tion hole" . As a test for spectral fluctuations of the GOE 
type, it has been applied much more widely in molecular 
physics (Leviandier et al. 1986, Guhr and Weidcnmiillcr 
1990b, Lombardi et al. 1994) than in nuclear physics (Al- 
hassid and Whelan 1993). Let 1 — Y%(b) denote the proba- 
bility of finding two levels at a distance b. For completely 
uncorrelated (Poissonian) spectra one has ^2(6) = for 
all b while GOE level repulsion implies 1^(0) = 1. The 
Fourier transform of the spectral autocorrelation function 
(a function of time i) depends on the Fourier transform 
of Yi and is sensitive to the difference between regular 
and chaotic motion. For chaotic motion it displays a 
"correlation hole" at t = 0. 



1. Neutron and Proton Resonances 

Although the fluctuation measures described in Sec- 
tion lll.Dl wcre proposed to describe neutron resonances in 
the 1950s, even in the early 1960s there were no neutron 
data of sufficient quality to provide an adequate test of 
RMT. For example, Dyson and Mehta (1963) considered 
the best available neutron resonance data and concluded 
the data were such that the RMT "model" was neither 
proved nor disproved. They exhorted experimentalists to 
improve the data quality. 
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By the early 1970s the high-quality neutron resonance 
data from the Columbia group( Liou et al, 1972a, 1972b) 
were available and seemed to confirm the predictions of 
the GOE version of RMT. However, due to the limited 
number of resonances (of order 100) for each nucleus, 
these results were considered suggestive but not defini- 
tive. Haq et al. (1982) and Bohigas et al. (1983) com- 
bined neutron resonance data sets from a number of nu- 
clei. This was made possible by scaling level spacings 
in units of the mean level spacing (the GOE fluctuation 
measures depend on that scaled parameter only). The 
authors also included some of the proton resonance data 
in their analysis. A complication (relative to the neu- 
tron data) was here the much larger energy range needed 
in order to obtain a reasonable sample size. This larger 
energy range required a correction (unfolding) of the ex- 
perimental data in order to transform to a new set of lev- 
els with constant mean level spacing, see Section III.D.3I 
The analysis included also suitable spectra for other than 
s-wave proton resonances. Using all these data, the au- 
thors obtained a set of 1407 levels that they labeled the 
nuclear data ensemble (NDE) . The analysis of the NDE 
is done by combining an energy average (for every nu- 
cleus) with an ensemble average (over all nuclei that are 
included in the NDE). Both the NNS and the A 3 statistic 
for the NDE agreed very well with the GOE predictions, 
see Figs. [TU] and QT] A number of later tests with other 
measures also agreed well with GOE predictions (see, for 
instance, Lombardi al., 1994) even though the presence 
of non-statistical effects can never be excluded (Koehler 
et al, 2007). 




FIG. 10 The NNS distribution for the nuclear data ensemble 
(histogram) and the GOE prediction (solid line) . From Bohi- 
gas et al., 1983. (By the time that paper was published, the 
nuclear data ensemble had grown to 1726 spacings). 

The analysis by Haq et al. (1982) and by Bohigas et 
al. (1983), and in some of the subsequent papers marks 
a turning point in the history of applications of RMT to 
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FIG. 11 The A3-statistic for the nuclear data ensemble (data 
points) and the GOE and GUE predictions (solid lines). The 
dashed lines estimate the finite-range-of-data errors. From 
Haq et al. (1982). 



nuclear spectra. As a result of these analyses it became 
generally accepted that proton and neutron resonances in 
medium-weight and heavy nuclei agree with GOE pre- 
dictions. With the later recognition of the connection 
between spectral RMT fluctuations and quantum chaos 
(see Section Hl.Fp . the term chaos began to be used by 
nuclear physicists. 

As mentioned in Section III.F1 chaos in classical many- 
body systems has not been investigated as thoroughly 
in terms of periodic-orbit theory as in classical few- 
degrees-of-freedom systems, not to speak of the com- 
plications due to the exclusion principle. Therefore, the 
connection between classical chaos and RMT is less well 
established and the use of the term "chaos" is some- 
what more tentative in nuclei. By the same token, the 
use of semiclassical periodic-orbit theory in nuclei has 
been basically limited to independent-particle motion. 
In that domain, it has been very successful. We mention 
early applications by Strutinsky (1966, 1967, 1968) whose 
"shell-correction method" is reviewed in Section HV.C.ll 
by Balian and Bloch, 1970, and recent work by Bohigas 
and Leboeuf, 2002; Leboeuf and Roccia, 2006; Roccia 
and Leboeuf, 2007. 

According to RMT, the eigenvalues and eigenvectors 
are uncorrelated random variables, see Eq. ©. In nu- 
clei, this prediction was tested by Bohigas et al. (1983). 
The correlation coefficient was found to be 0.017 ±0.029, 
the error reflecting the finite number of data points. In 
microwave billiards, the test yields 0.02 ± 0.05 (Alt et 
al., 1995). Another test (also giving agreement with the 
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GOE) has been reported in molecules by Lombardi and 
Seligman (1993). 



2. Low-lying Levels 

After the success of RMT in describing the fluctua- 
tion properties of highly excited (resonance) states, it was 
natural to attempt to extend such analyses to low-lying 
states. Although there is an enormous amount of exper- 
imental information available for states near the ground 
state, for most nuclides the quantum numbers are known 
for only a very limited number of low-lying states. In 
particular, complete and pure sequences of levels with 
the same spin and parity are typically very short. There- 
fore data from several or many nuclei must be combined 
to generate a sufficiently large ensemble. Moreover, the 
shortness of the available sequences precludes the study 
of other fluctuation measures than the NNS distribu- 
tion. An extensive data set was compiled by von Egidy 
et al. (1986). An initial analysis of a subset of these 
data was performed by Abul-Magd and Weidenmullcr 
(1985). A more extensive analysis of this same data 
set was performed by Shriner et al. (1991). The spac- 
ing distributions and their cumulative sums are shown in 
Fig. [T2j The nuclei are grouped into classes according to 
mass number A. The size of each class was determined 
by the data available. Obviously, the cumulative sums 
have smaller fluctuations. While the agreement with 
the Wigner distribution looks satisfactory in most cases, 
clear deviations occur for 150 < A < 180 (rare earth 
nuclei) and for 230 < A (very heavy nuclei). In both 
ranges of mass numbers, sizable nuclear deformations oc- 
cur and cause rotational motion, see Section lTV.B.ll The 
rotational model is integrable and the motion therefore 
regular. The same statement holds for other forms of 
so-called collective motion, see Section IIV.B1 

The spacing distributions were fit with the Brody dis- 
tribution (Brody et al., 1981). The overall trend of the 
Brody parameter u> was to decrease with increasing mass 
number A - for the lightest region (A = 25 - 50) the 
average value of w was about 0.7, while for the heavi- 
est mass region (A — 225 - 250) the value of u> was 0.2. 
Various theoretical works have attempted to explain this 
behavior, including Bae et al. (1992) and Yoshinaga et 
al. (1993). To exhibit the connection of the NNS dis- 
tributions with the degree of collectivity, attention was 
focused on the behavior of the 2 + and 4 + states because 
these states play a prominent role in collective rotations 
and vibrations of the nucleus. Two types of nuclei were 
considered: nuclei with approximately spherical ground 
states and nuclei with strongly deformed ground states, 
see Section IIV.BI The transition between both classes 
was studied. Depending on the model chosen, the mo- 
tion is chaotic or regular in one but not in the other limit. 
Fig. 1131 shows that the experimental results are striking. 
2 + and 4 + states in strongly deformed nuclei have NNS 
distributions that agree with the Poisson distribution, 



while the corresponding states in spherical nuclei have 
spacing distributions that agree with the Wigner distri- 
bution. Unfortunately, the limited amount of data and 
the corresponding large uncertainties preclude a more de- 
tailed assessment of the effects of collectivity. 

Another approach by Abul-Magd et al. (2004) focuses 
on specific states (the lowest 2 + states in even-even nu- 
clei) and collects all complete sequences of low-lying 2 + 
states from the nuclear data tables. The sequences are 
short in most cases. Nuclei are classified by the ratio 
-R4/2 of the excitation energy of the lowest 4 + state over 
that of the lowest 2 + state. That ratio is a well-known 
measure of collectivity. The size of each class is chosen 
such that it contains a sufficient number of sequences 
for a meaningful statistical analysis. The NNS distri- 
butions are analyzed using a measure different from the 
Brody distribution and obtained by superposing a num- 
ber of uncorrelated GOE sequences, each of mean frac- 
tional level density /. The parameter / serves as a fit 
parameter and is referred to as "chaoticity parameter". 
In its dependence on R4/2, this parameter displays deep 
minima when R4/2 equals 2.0, 2.5, and 3.3. These val- 
ues correspond to the dynamical symmetries of a specific 
collective model, the interacting boson model described 
in Section IIV.B.U Whenever one of these symmetries 
prevails, the motion of the nucleus is integrable and thus 
regular. 

In summary, there is evidence that the nuclear dynam- 
ics in the ground-state region is partly chaotic and partly 
regular. The regular features are dominant whenever col- 
lective motion with a high degree of symmetry applies. 



3. High-spin States 

The only extensive analysis of the statistical proper- 
ties of high-spin states was performed by Garrett et al. 
(1997). The data set comprised energy levels in deformed 
nuclei in the range of proton numbers Z = 62 - 75 and 
mass numbers A = 155 - 185. The spin values ranged up 
to J = 40. The levels were at high excitation energies, 
but only up to several hundred keV above the yrast line, 
see Fig. |5J The authors found that the NNS distribution 
agreed best with the Poisson distribution. This is con- 
sistent with the results for low-lying states in deformed 
nuclei, see Section HV.B. 21 

An observed deficiency of small spacings is not well 
understood. It is possible that right above the yrast line 
a symmetry related to the K quantum number is partly 
broken. This quantum number measures the projection 
of the nuclear spin onto the body-fixed symmetry axis, 
see Section llV.B.ll This might lead to level repulsion at 
small distances. Hopefully the analysis of similar data 
will shed more light on that question. We return to the 
general issue of symmetry breaking and its influence on 
spectral fluctuation properties in Section IIII.D.ll below. 
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FIG. 13 Comparison of the NNS distributions versus x, the 
level spacing in units of the mean level spacing, for 2 + and 
4 + states in strongly deformed (left panel) and in spherical 
nuclei (right panel). Adapted from Shriner et al. (1991). 



FIG. 12 Upper panels: NNS distributions (histograms) are 
compared with the Wigner distribution (solid lines) and the 
Poisson distribution (dashed lines) for several ranges of mass 
numbers. Lower panels: same for the cumulative distributions 
(number of spacings smaller than x). (Cf. the caption of 
Fig. [7J| . The approximate number of levels in each region is: 
A = 0-50 (N = 150), A = 50-100 (N = 50), A = 100-150 (N 
= 270), A = 150-180 (N = 450), A = 180-210 (N = 60), A> 
230 (N = 190). From Shriner et al. (1991). 



4. Analysis of Low-Lying Modes of Excitation 

At low excitation energies, one observes several modes 
of excitation. Statistical measures have been used in or- 
der to identify the character of the mode in the case of 
the scissors mode (Enders et al., 2000) and of the pygmy 
dipole resonance (Enders et al, 2004). 

For the scissors mode (a low-lying isovector magnetic 
dipole mode) data were generated by nuclear resonance 
fluorescence measurements (Enders et al., 2000). The 
spectra of 13 heavy deformed even-even nuclei with neu- 
tron numbers in the 82 - 126 range (corresponding to a 
major shell) were used to generate an ensemble of 152 
scissors-mode states with spin/parity 1 + , all in a range 
of excitation energy between 2.5 and 4.0 MeV. In each 
nucleus, the sequence of states used in the analysis was 
required to contain a minimum number of 8 states. After 
unfolding, the ensemble was analyzed with the standard 
RMT fluctuation measures. The data agreed very well 
with Poisson statistics. The authors examined carefully 
the effects of missing levels on the spacing and width dis- 
tributions and concluded that missing levels can be ruled 
out as a cause of this behavior. They concluded that the 
levels of the scissors mode are excited by a common mech- 



anism. The levels are collective but it is not possible to 
identify a common doorway. It seems that the underlying 
microscopic mechanism is not yet fully understood. 

The electric pygmy dipole resonance is so named be- 
cause of its small strength relative to the giant electric 
dipole resonance. In heavy nuclei, the pygmy resonance 
is located at excitation energies around 5 to 7 MeV. En- 
ders et al. (2004) studied the statistical properties of 
this mode in four isotones, all with neutron number 82. 
They created an ensemble of 184 1 _ states in the ex- 
citation energy range of 4 to 8 MeV, along with their 
dipole transition strengths to the ground state. After 
unfolding, the spectral fluctuations (strength and spac- 
ing distributions) are close to Poissonian. Because of a 
significant number of missing levels the analysis is rather 
involved in this case, however, and the authors conclude 
that the weak correlations found point to GOE behavior 
of the complete spectra. That conclusion is reinforced 
by an extensive comparison with spectra calculated us- 
ing a particular nuclear model, the quasiparticle phonon 
model. These agree with the data but yield GOE behav- 
ior for the full spectra (including the levels missing in the 
data). The fundamental mode of excitation is collective. 
Many other states exist at the excitation energy where it 
occurs. These fragment the doorway state and produce 
a correlated spectrum of GOE type. 

In comparing their results for these two modes, the 
authors conclude that the key reason for the apparently 
different statistical behavior is the difference in excitation 
energy. The higher-energy mode (electric pygmy dipole) 
is in a region of greater level density. This results in 
correlated spectra. 
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5. Eigenvector Distribution 

The GOE predicts a Gaussian distribution for the pro- 
jections of the eigenvectors and the Porter-Thomas dis- 
tribution for their squares, see Section Hl.D.ll While both 
the NNS distribution and the A3-statistic are strongly 
affected by missing levels and/or impure sequences, the 
Porter-Thomas distribution (a probability density and 
not a correlation function) is expected to be less sensi- 
tive to missing or wrongly assigned levels. The experi- 
mental data normally used are the reduced widths (for 
resonances) or the reduced transition strengths (for elec- 
tromagnetic transitions). Early neutron resonance data 
appeared to agree with the Porter-Thomas distribution; 
a frequently quoted example is for neutron resonances on 
232 Th measured by Rainwater's group at Columbia Uni- 
versity (see Garg et al., 1964). The analysis by the Orsay 
group of the nuclear data ensemble included a test of the 
Porter-Thomas distribution for the widths (see Bohigas 
et al, 1983), and used a total of 1182 measured widths. In 
addition to a direct comparison with the Porter-Thomas 
distribution, a search for the best x 2 distribution was also 
done. Very good agreement with the GOE prediction was 
found. 

Further and more detailed attempts to confirm the 
Porter-Thomas distribution have run into the follow- 
ing difficulty. For a set of Gaussian-distributed am- 
plitudes {7^} the second and the fourth moments are 

related by (7 4 ) = 3 (7 2 ) where the angular brack- 
ets denote the running average. According to Harney 
(1984), the error (square root of the variance) of the ra- 
tio R = ^7 4 ) /3 (7 2 ) is i/8/(3fc), where k is the number 
of data points. This is a rather large value. For the often 
quoted 232 Th data, for example, k = 171, and the error 
of R is 0.125. Thus this excellent data set only confirms 
the Gaussian nature of the amplitude distribution at the 
12 % level. 

To overcome this problem, a larger data set seemed 
useful. An ensemble of 1117 reduced widths was formed 
with TUNL data (Shriner et al, 1987). With y = 7 2 /(7 2 ) 
and P(y) the Porter-Thomas distribution, the result for 
s/yP(y) is shown in Fig. 1141 The visual agreement with 
the GOE prediction is striking. However, the value of R 
for this ensemble turned out to be 1.26. The problem 
is that a relatively small number of non-statistical large 
widths has a major impact on R because R depends on 
the fourth moment of the amplitudes. 

To overcome that difficulty, a different measure which 
is less sensitive to a few non-statistical amplitudes was 
needed. The normalized linear correlation coefficient 
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for two data sets {xi} and {x'j} is expected to be a sen- 
sitive measure of correlations, as it combines informa- 
tion on both the magnitudes and the phases of the data 
points. To test whether two sets of amplitudes {a^} and 
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FIG. 14 The distribution ^/yP{y) for y = 7V7 2 for 1117 
reduced widths. From Shriner et al. (1987). 



{a,;} follow the Gaussian distribution, one calculates the 
amplitude correlation coefficient p(a, a') and the width 
correlation coefficient p(w,w') where w stands for the 
square of the amplitude. (Since it is impossible to mea- 
sure the absolute sign of an amplitude, p(a,a!) is calcu- 
lated with the assumption that (a) = = (a')). The 
Gaussian distribution predicts p 2 (a,a') = p(w,w') and 
this is the relation which is tested. 

To work out the correlation coefficients it is necessary 
to measure partial width amplitudes including their rel- 
ative phases. That was done using the inelastic decay of 
proton resonances. For example, for a spin-zero target 
with mass number A and a 2 + first excited state, a pro- 
ton resonance in the nucleus with mass number A+l and 
spin/parity 3/2~ can decay to the 2 + state by emitting 
the proton with angular momentum/spin P1/2 or ^3/2- 
The relative phase of the two decay amplitudes is de- 
termined by measuring the angular distribution of the 
inelastically scattered protons and the subsequent de- 
excitation 7 rays. Application of this approach is de- 
scribed in detail in a review by Mitchell et al. (1985). 

The relation p 2 (a,a') — p(w,w') can be checked for 
each data set. However, the sample size for each spin 
value in each nuclide was small, and it was necessary to 
combine the data sets. To this end, an ambiguity in the 
definition of the correlation coefficient was used. Instead 
of writing the amplitudes in terms of the decay channels 
P1/2 and P3/2, any two linear combinations of amplitudes 
obtained by orthogonal transformations from the original 
set can be used. The correlation coefficient depends on 
the chosen representation. There is always a represen- 
tation in which the amplitude correlation coefficient is 
zero. Each of the data sets was individually transformed 
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to that representation. The resulting width correlations 
and their uncertainties were combined for all of the data 
sets in order to obtai n a final value for the width correla- 
tion. The result was p(w,w') = -0.01 ± 0.03, (Shriner et 
at, 1987, 1989) , in agreement with the prediction based 
on the Gaussian distribution. This is the most sensitive 
test of the Gaussian assumption. 



D. Violation of Symmetry or Invariance 



the same quantum numbers J and II (which we omit) 
but different isospin quantum numbers T\ and T%. The 
Hamiltonian has the form 



Two of the symmetries mentioned in Section IIII.AI 
hold only approximately: Isospin symmetry is broken by 
charge effects, and parity conservation is violated by the 
weak interaction. Is it possible to extend RMT so as to 
account for such symmetry breaking? And how do the 
resulting statistical measures compare with data? We 
answer these questions here for the case of isospin sym- 
metry breaking. Parity violation is a very weak effect 
that has so far received a statistical analysis only in the 
framework of nuclear reaction theory, see Part 2 of this 
review. We apply the results to the complete spectra of 
26 Al and 30 P. 

We have also assumed that nuclei obey time-reversal 
invariance. One of the most precise tests in nuclei of that 
assumption applies RMT. We describe how a measure 
for violation of time-reversal invariance is derived in the 
framework of RMT and applied to data. 



1. Model for Isospin Violation 

Isospin symmetry is violated in nuclei by charge- 
dependent effects such as the Coulomb interaction be- 
tween protons, the neutron-proton mass difference, or 
the mass differences between charged and neutral pions. 
Needless to say, the isospin-violating interaction is small 
compared to the strong force. The breaking of isospin 
symmetry manifests itself differently in different ranges 
of mass numbers. In nuclei with mass numbers around 
40 or more, it leads to the occurrence of fragmented iso- 
baric analog resonances. The typical features of this phe- 
nomenon relate to nuclear reaction theory and are not 
dealt with here (see Part 2 of this review). In some 
light nuclei with mass numbers smaller than 40 or so, 
the ground state has isospin T — but the density of 
states with T = 1 in the ground-state region is roughly 
the same as that of states with T = 0. The charge- 
dependent forces mix states with T — and T = 1. It is 
for these nuclei that the following random-matrix model 
applies. 

In order to account for isospin violation (and for sym- 
metry breaking in general) , the Hamiltonian (|27[) must be 
modified. The matrix elements of the isospin-violating 
interaction couple states with different T-quantum num- 
bers. For simplicity of presentation we consider two diag- 
onal blocks only. The isospin-breaking interaction con- 
serves parity and total spin so these two blocks carry 
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With Ni (N 2 ) the dimensions of the two block-diagonal 
matrices, the running indices in the first (second) block 
are p, v — 1, . . . , iVi (p, a = N\ + 1, . . . , Ni + N 2 ), respec- 
tively. We deviate from our earlier systematic notation 
and denote the coupling matrix elements connecting the 
two blocks by V. 

We use Eq. (|3"0|) to define a random-matrix model for 
symmetry breaking (Rosenzweig and Porter, 1960). We 
assume that the matrices H Tl and H T2 are each mem- 
bers of a GOE and are uncorrected. For simplicity we 
assume that the two GOEs have identical semicircle radii 
2Ai = 2A2 = 2A and equal dimensions Ni = N 2 = N al- 
though in practice it is necessary to take iVi j^ N 2 in 
order to account for the fact that the level densities for 
states with different isospins differ. The (real) matrix el- 
ements of V are assumed to be Gaussian-distributed ran- 
dom variables with zero mean value and a common sec- 
ond moment V 2 . They are not correlated with each other 
or with the elements of either of the two GOEs. Strictly 
speaking, the symmetry-breaking interaction also con- 
tributes to the diagonal blocks H Tl and H T2 . However, 
such contributions can be incorporated in the random- 
matrix description of these blocks and therefore do not 
appear explicitly in the model. 

To define the strength of \" 2 we recall Eq. (4). That 
equation seems to suggest that we put V 2 — a 2 X 2 /N, 
with q 2 < 1 to account for the weakness of the 
symmetry-breaking interaction. That is not correct, 
however, and we must choose V 2 = a 2 X 2 /N 2 , with a a 
strength parameter which is independent of N in the limit 
N — > 00. Indeed, without symmetry breaking, the spec- 
tra of H x and H 2 are uncorrelated. A weak symmetry- 
breaking interaction induces level repulsion and stiffness 
among levels with different isospin. That happens when 
the matrix elements of V are of the order of the mean 
level spacing d = nX/N, or when V 2 is of the order of 
A 2 /TV 2 . Hence the ratio of the strength of the symmetry- 
breaking interaction (average of the square of the matrix 
elements) over that of the symmetry-conserving interac- 
tion vanishes asymptotically as 1/N for N — > 00, since 
the mean level spacing likewise vanishes in that limit. 
Conversely, symmetry violation becomes detectable in 
spectral fluctuation measures when the matrix elements 
of the symmetry-breaking interaction are of the order of 
the mean level spacing. This condition was already men- 
tioned below Eq. ([5]) and is met by the isospin-violating 
matrix elements in light nuclei, see Section HlI.D.2l bclow. 
The same conclusion is reached when we consider the vi- 
olation of time-reversal invariance in Section IIII.D.4I be- 
low. Apparently Pandey (1981) was the first author to 
note the relevance of such small parameters for violations 
of symmetry and/or invariance. 
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The parameter V 2 suffers from the same shortcoming 
as discussed in Section lil.Glfor the mean-square-matrix 
element of a doorway state: V 2 changes strongly with ex- 
citation energy and/or mass number. A much smoother 
measure of symmetry breaking is the spreading width, 
see Eq. (|26p . There are two possible definitions for T+, 



2. Complete Level Schemes 
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where p T is the level density for the states with isospin T. 
It is largely a matter of convenience which of these defi- 
nitions is used. To interpret the spreading widths, let us 
consider without loss of generality the first of Eqs. (|3ip. 
We use a basis in which both H Tl and H T2 are diago- 
nal. We first assume that the mean level spacing l/p Tl of 
states with isospin T\ is significantly larger than T^ Tl and 
that T^ Tl is significantly larger than l/p T2 . Then the ar- 
guments of Section Hi. Gl apply, each state with isospin T\ 
acts as an isolated doorway state, and T^ is the average 
width of the probability distribution for finding the door- 
way state mixed into the eigenstates of the full system. 
We expect that that interpretation remains qualitatively 
valid also when the inequalities l/p Tl 3> T^ 3> p T2 are 
violated. In other words, I^ Tl (T^ T2 ) is a measure of the 
width in energy with which every eigenstate of H Tl (of 
H Tl , respectively) is spread over the eigenstates of the 
full system. 

The random-matrix ensemble (|30p lacks the overall or- 
thogonal symmetry of the GOE. By construction, the en- 
semble is invariant, of course, under orthogonal transfor- 
mations of the first Aq (the last N%) states, respectively. 
Nonetheless, the analytical treatment of symmetry viola- 
tion in RMT is much more difficult than treating a single 
GOE. While analytical results for the ensemble ([50)1 are 
not available, replacing each of the two block-diagonal 
GOEs in Eq. (|30p by a GUE, and considering the com- 
plex matrix elements in the non-diagonal blocks as un- 
corrected Gaussian-distributed random variables with 
zero mean value and common variance V 2 , one arrives 
at a tractable problem (Guhr and Weidcnmuller, 1990a). 
Although belonging to a different symmetry class, the re- 
sulting ensemble is expected to possess features which are 
qualitatively similar to those of the ensemble (|30| . Ad- 
ditional information is generated by a numerical analysis 
of the ensemble (|30|) . 

As V 2 /d 2 oc a 2 increases from zero, the spectral fluc- 
tuations change from those of two uncorrelated GUEs to 
those of a single GUE. From V^/d 2 = to V^/d 2 > 0, the 
change is discontinuous: Level repulsion amongst states 
with different T sets in suddenly. For the local spec- 
tral fluctuation measures considered by Guhr and Wei- 
denmiiller (1990a) the case of a single GUE is attained 
when V 2 /d 2 « 1. We expect this statement to apply 
likewise to other fluctuation measures, and to hold sim- 
ilarly for the ensemble in Eq. (j3"0)) . The expectation is 
confirmed by numerical simulations. 



As described in Section lIlI.B.5[ complete level schemes 
were determined for the nuclides 26 A1 and 30 P. These 
odd-odd N = Z nuclei are particularly interesting be- 
cause here the densities of states with T = and with 
T = 1 arc almost equal, starting from the ground state, 
while in most nuclei the states with higher T are shifted 
toward higher excitation energies. Thus these nuclei are 
ideal to study the effect of isospin-symmetry breaking. 
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FIG. 15 The four panels show the A3-statistic for the states 
in 26 A1 obtained by taking into account only the quantum 
numbers indicated at the top of each panel. 

A qualitative test for a conserved symmetry is to con- 
sider what happens when that symmetry is neglected. As 
an example, in Fig. [T5] the A3-statistic is shown for the 
states in Al. Only the quantum numbers shown at the 
top of each panel are taken into account in evaluating 
A3. We note that ignoring the good quantum number J 
leads to a major increase in A3. Ignoring T, on the other 
hand, leads to a very small change of A3. This seems 
to suggest that isospin is not a good quantum number, 
in apparent contradiction to other evidence that isospin 
is only slightly broken (at about the 3 % level) in this 
nuclide. 
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The point is that (as shown in Section HlI.D.l[) . the im- 
pact of symmetry breaking on A 3 depends on the ratio 
of the symmetry-breaking matrix element to the mean 
level spacing d and is thus enhanced by small spacings. 
Hence, a small degree of symmetry breaking can have a 
large effect on the statistical measures. Somewhat fortu- 
itously, the strength of symmetry breaking in 26 A1 is such 
that the A3-statistic lies between the values for a single 
GOE and for two GOEs, see Fig. [16] More precisely, 
the root-mean-square value of the symmetry-breaking 
Coulomb matrix elements in 26 A1 is a little smaller than 
but of the order of the mean level spacing (Guhr and 
Weidenmiiller, 1990a). The results were consistent with 
other experimental determinations. Experimental results 
on isospin violation for 30 P were almost identical with the 
26 Al results (Shriner et al., 2000). 



0.4 






_-"' 


0.3 




^"l" -±^xT^ 


, 






,{'O^C~- 


—--•"" . 








IJ.J 






0. 1 


/ 




0.0 







FIG. 16 The spectral rigidity A3 versus L for 75 levels with 
T — and 25 levels with T = 1 measured in Al (dots with 
error bars). The excitation energies lie between zero and 8 
MeV. The lower (upper) dashed line is the prediction for a 
single GOE (for two uncorrelated GOEs with fractional den- 
sities 3/4 and 1/4, respectively). The solid line is the result 
of a numerical simulation incorporating symmetry breaking. 
The strength of the symmetry-breaking interaction was fitted 
to the data. From Guhr et al. (1998). 

Although these results were consistent with theoreti- 
cal expectations, they were not considered definitive due 
to the small sample sizes. Definitive results on symme- 
try breaking were provided by measurements of acous- 
tic resonances in quartz blocks (Ellegaard et al., 1996) 
and of electromagnetic resonances in coupled microwave 
billiards (Alt et al., 1998). In these measurements the 
strength of the symmetry breaking force could be effec- 
tively varied and much larger sample sizes were obtained. 
The agreement with the RMT model of Eq. ([30]) was ex- 
cellent. 

The third nucleus with a "nearly complete" level 
scheme (for excitation energies below 4.3 MeV) is 116 Sn. 



In that nucleus, the NNS distribution was studied (Ra- 
man et al., 1991). Only sequences with a minimum of five 
levels with the same spin and parity assignments were 
included in the analysis; there were 6 such sequences. 
The histogram for the NNS distribution was fit with the 
Brody parametrization (|28|) . The fit gave u> = 0.51±0.19, 
similar to the best fit in 26 A1 which gave w = 0.47 ± 0.14. 
In 26 A1 the deviation from GOE predictions is due to 
isospin symmetry breaking. The cause for the same phe- 
nomenon in 116 Sn is not clear. The nucleus 116 Sn is men- 
tioned in the present Section only because of its "nearly 
complete" level scheme. 



3. Transition Strengths 

After the studies of isospin symmetry breaking in the 
complete spectra of 26 A1 and 30 P (see Section lrn.D.2j) . 
the effect of symmetry breaking on the eigenvectors in 
the same systems was also explored. Although there 
was no formal proof, heuristic arguments predicted that 
the Porter-Thomas distribution would not be changed 
by symmetry breaking. The central point of the argu- 
ment was the complexity of the wavefunctions of initial 
and final states of the transitions. 

The TUNL group (Adams et al, 1998; Shriner et al, 
2000) used the reduced electromagnetic transition prob- 
abilities B in these nuclei as the data set. To elimi- 
nate issues of scale, subgroups of the transitions were 
considered. The transition probabilities were classified 
according to multipolariy (electric dipole El, electric 
quadrupole E2, magnetic dipole Ml) and to the isospin 
difference AT — or 1 between initial and final states. 
The parameter y = B / B turns out to have a very 
large dynamic range. Therefore it is convenient to use 
z = log 10 y. In terms of that variable, the upper panels 
of Fig. [TTlshow the Porter-Thomas distribution as solid 
lines. The lower panels give the integrated distributions. 
The histograms show the data for 30 P. The details of the 
analysis are given by Adams et al. (1998), and Shriner 
et al. (2000). Similar results were obtained for 26 A1. It 
is obvious that the data do not agree with the Porter- 
Thomas distribution. 

It was some years before this somewhat unexpected 
result was formally analyzed and explained. Barbosa et 
al. (2000) used basically the same approach as for the 
description of symmetry breaking on the level statistics, 
see Section IIII.D.ll We only sketch the central point. 
We consider electromagnetic transitions of a given mul- 
tipolarity (typically Ml or E2). The model (001) has to 
be extended because the electromagnetic transition op- 
erator does not, in general, conserve either T or J and 
connects states with different spins as well as different 
isospins. For V 2 — the eigenfunctions of H are eigen- 
functions either oiH JlTl or ofH j2T2 . For fixed J x , J 2 the 
transition matrix elements belong to one of three classes: 
(i) Those coupling two states with isospin T\\ (ii) those 
coupling two states with isospin T%; and (iii) those cou- 
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4. Test of Time-reversal Invariance 




_l 
1 


U 


V 




n 


-1- 


\^*\ .: 



z-log 10 (B/<B>) 



-4 -S S 

log, n (B/<B>) 



1 -4 -S 2 
z=log 10 (B/<B>) 




s-1ob„(B/<B>) 



-4 -3 
z-lo£ 10 (B/<B>) 



z-loi 10 (B/<B>) 



FIG. 17 Comparison between the Porter-Thomas distribu- 
tion (written in terms of the variable z = log 10 y) (solid lines) 
and the data (histograms) for several multipole transitions 
and isospin differences in 30 P. From Shriner et al, 2000. 



pling two states with different isospins. The distribution 
of matrix elements within each class is expected to be ap- 
proximately Gaussian, but the three Gaussian distribu- 
tions may have different heights and widths. Moreover, 
transitions with different multipolarities behave differ- 
ently. Therefore the squares of the transition matrix el- 
ements in Fig. [17] cannot have a simple Porter-Thomas 
distribution even for V 2 = 0. As V 2 increases from zero, 
the Gaussian distributions get mixed. Details are only 
accessible numerically. With a somewhat different ap- 
proach Hussein and Pato (2000) also predicted a devia- 
tion from the Porter-Thomas distribution. Neither group 
attempted to fit the experimental data in detail. Such an 
analysis is still missing. We conclude that while the ef- 
fect of isospin symmetry breaking on the level statistics 
is a generic phenomenon that can be accounted for com- 
pletely by a simple extension of RMT in terms of a single 
parameter (the spreading width) , the effect of symmetry 
breaking on the electromagnetic transition strengths in- 
volves additional dynamic elements (classification of the 
transition matrix elements). 

While the size of the data set in 26 A1 and 30 P is 
limited, experiments using coupled microwave billiards 
(Dembowski et al., 2005) yielded data with much better 
statistics. The distribution of transition strengths devi- 
ates from the Porter-Thomas distribution. The theory 
of Barbosa et al. (2000) was extended to this case by 
Dietz et al. (2006). 



Because of the anti-unitarity of the time-reversal op- 
erator, the modeling of a violation of time-reversal in- 
variance in RMT is fundamentally different from that of 
a broken symmetry as discussed in Section IIII.D.ll As 
explained in Section fll.Al time-reversal invariance allows 
us to choose the Hamiltonian matrix as a real and sym- 
metric matrix. The matrix ensemble that models such 
systems is the GOE. If time-reversal invariance does not 
hold, the Hamiltonian matrix is Hcrmitian but cannot, 
in general, be chosen real and symmetric. The matrix 
ensemble that models such systems is the GUE, see Sec- 
tion III.B1 To describe the violation of time-reversal in- 
variance in RMT, we need to construct an ensemble that 
interpolates between the GOE and the GUE. This is done 
as follows. Every Hermitian matrix can be written as the 
sum of a real symmetric matrix and of i times a real anti- 
symmetric matrix. A stochastic model for a Hamiltonian 
with some violation of time-reversal symmetry is then 



H = 



1 



y/1 + (l/N)a 2 



H 



GOE 



+ N- 1/2 aiA 



(32) 



Here H GOE stands for the GOE, and the elements of the 
real antisymmetric matrix A are Gaussian-distributed 
random variables with zero mean value and a second mo- 
ment given by 



A 2 



-A^i/A.po- — \dfipdi/cr dfi<jdi>p) 



(33) 



The elements of A and of H GOE are uncorrelated. The 
real dimensionlcss parameter a describes the strength of 
the violation of time-reversal invariance. For a = we 
deal with the GOE and for a = TV 1 / 2 , with the GUE. 

For tests of the violation of time-reversal invariance 
in nuclei, the central feature of the GUE is level repul- 
sion at small distances (scaled spacing s < 1). In the 
case of the GOE, level repulsion leads to a linear de- 
pendence of the NNS distribution for small s, see the 
Wigner surmise (TlT)|) . In contradistinction, the NNS dis- 
tribution for the GUE increases quadratically with s for 
small s. A test for a violation of time-reversal invariance 
in nuclei is, therefore, based on a detailed examination 
of the NNS distribution at small spacings (French et al., 
1985). As explained below Eq. © and, in a different 
context, in Section IIII.D.ll the mixing of levels (and, 
thus, the spacing distribution) are sensitive to very small 
mixing matrix elements which are of the order of the 
GOE mean level spacing d — ttX/N. For the perturba- 
tion N~ 1 ' 2 aiA in Eqs. (j32l |33|) . the root-mean-square 
matrix element has the value aX/N. This is comparable 
to d for awl which explains the choice of factors A 1 / 2 
in Eq. ([3"2")l . As in Section [III. D.ll the mixing parame- 
ter N~ x / 2 a, vanishes asymptotically as 1/y/N. The NNS 
distribution is sensitive to a violation of time-reversal 
symmetry when the relevant matrix elements are of the 
order of the average nuclear level spacing. The analy- 
sis was done (French et al., 1985) for the nuclear data 
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ensemble (see Section fill. C.l[) . where typical spacings d 
are of the order of 10 keV, and yielded an upper bound 
of about d/10 for the time-reversal non-invariant matrix 
element of the nuclear Hamiltonian. From this result 
French et al. (1985) inferred an upper bound of about 
1 percent for the time-reversal non-invariant part of the 
nucleon-nucleon interaction. 



IV. CHAOS IIM NUCLEAR MODELS 

Sections [TTI and HTT1 were based almost entirely on con- 
cepts of RMT and made very little use of the wealth of 
information on the dynamical behavior of nuclei. We fill 
this gap in the present Section. We discuss the two lead- 
ing nuclear-structure models which describe phcnomcno- 
logically the dynamics of nuclei: The shell model (which 
mostly applies to spherical nuclei), and the collective 
model (which mostly applies to nuclei with surface de- 
formations) . We give a brief introduction to both mod- 
els which, in their simplest form, are fully integrable and 
thus give rise to regular motion. We present evidence 
that both models also allow for chaotic motion. We then 
turn to a number of specific applications of RMT which 
incorporate dynamical aspects. 



A. Spherical Shell Model 

1. The Nuclear Shell Model 

In the elementary version of the nuclear shell model 
nucleons move independently under the influence of a 
common mean field. Attempts to introduce a mean field 
into nuclear theory date back to the 1930s. These at- 
tempts were unsuccessful at the time, partly because they 
lacked an essential ingredient (the strong spin-orbit cou- 
pling) and partly because of the success of Bohr's idea 
of the compound nucleus which depicted nuclei as sys- 
tems of strongly interacting particles. The introduction 
in 1949 of a central single-particle potential with strong 
spin-orbit coupling by Haxel et al. (1949) and Mayer 
(1949) changed that situation. That model was very 
successful in the description of nuclear properties in the 
ground-state domain and shifted attention away from 
the compound-nucleus picture. It gave rise to a burst of 
spectroscopic activity which lasted for many years and 
thoroughly validated the model. 

The sequence of single-particle levels of the nuclear 
shell model is shown in Fig. [IS] The scheme is fundamen- 
tally that of the harmonic oscillator in three dimensions 
with excitation energies nhuj and n — 0,1,2,.... The 
integer n defines the major shells. Individual levels are 
denoted by (n + 1), by the single-particle angular mo- 
mentum £ in spectroscopic notation, and by the single- 
particle total spin j obtained by vector-coupling the an- 
gular momentum operator £ and the spin s. (Here we 
do not distinguish neutrons and protons and take isospin 



as a good quantum number. The picture requires some 
modification for medium-weight and heavy nuclei which 
we do not address). The degeneracy of the single-particle 
states in each major shell, which is characteristic of the 
harmonic oscillator, is lifted because the single-particle 
potential does not have the shape of a harmonic oscilla- 
tor, and because of the presence of a strong spin-orbit 
coupling that pushes the states with highest spin in each 
major shell down into the next-lower major shell (except 
for the lowest shells where the spin-orbit interaction is 
not strong enough). Thus each major shell contains a 
number of subshells each of which is characterized by the 
quantum numbers (n + l,£,j). Within a given major 
shell the index n is redundant and will often be omitted. 
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FIG. 18 Level sequence in the nuclear shell model. 
Mayer and Jensen (1955). 
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For a nucleus with mass number A, the A nucleons fill 
the lowest shells in accord with the exclusion principle. 
The completely filled shells are considered as inert, and 
the spectroscopic properties of the low-lying states result 
from the m "valence nucleons" which partly fill the last 
shell (the "valence shell" ) . This scheme accounts for the 
strong binding energies of nuclei with closed shells, and 
for the ground-state properties of nuclei differing from 
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closed-shell ones by the addition or removal of a single 
nucleon. For nuclei with more than one valence nucleon 
(or more than one hole in the valence shell) , an extension 
of the model is called for because in such cases strong 
degeneracies occur in the model which are not observed in 
reality. For instance, two valence nucleons in the lds/2 - 
subshell of the 2sld-shell can be coupled to total isospin 
T = or T = 1. For T = (T = 1), the possible 
states have odd (even) total spin values J ranging from 
J = to J = 5. All of these states are degenerate in the 
elementary shell model. 

To remedy that situation, the elementary shell model 
is viewed as a mean-field theory which takes account of 
most (but not all) of the nucleon-nucleon interaction. 
The remaining "residual interaction" must be included 
to obtain quantitative agreement with data. From the 
point of view of many-body theory, the residual interac- 
tion is an effective interaction. The residual interaction is 
usually assumed to be a two-body interaction (although 
there is evidence (Pieper and Wiringa, 2001) that three- 
body forces are needed in some cases to obtain good fits 
to the data), to be time-reversal invariant, and to con- 
serve spin and parity. In the spirit of the shell model it 
is also assumed that the residual interaction is weak (it 
does not significantly mix the many-body states belong- 
ing to different major shells), so that it effectively acts 
only amongst the m valence nucleons in the valence shell. 

The input for this model (the full shell model or, in 
brief, the shell model) consists of the single-particle en- 
ergies E£j in the valence shell, and of the matrix elements 
of the residual interaction in that shell. The residual 
two-body interaction V,. cs is completely characterized by 
a finite number of two-body matrix elements. The two- 
body states \j1j2star) are obtained by coupling any two 
single-particle states (iiji) and ((2J2) in the valence shell 
to total spin s and total isospin t with z-components a 
and t. The antisymmetrized two-body matrix elements 
of V les have the form {J3JiSt\V Ies \jiJ2st). Our notation 
implies conservation of spin and isospin and takes ac- 
count of the fact that the values of the matrix elements 
do not depend on er and r. Conservation of parity im- 
poses an additional constraint not explicitly displayed in 
our notation. For brevity we will refer to these matrix 
elements by the symbol v a . The index a enumerates all 
allowed and distinct (i.e., not connected by symmetry) 
two-body matrix elements in the valence shell. For the 
2s Id-shell, the range of a is 63 while for the 2pl/-shell, 
it is 195. 

In second quantization, the shell-model Hamiltonian 
governing the m valence nucleons is then given by 



Ti 



Ij err 

+ 4 ^2 (JzjiSt IVrcsl.7lJ2.si) 



jlJ2J3J4St 

xYAt 



-4, 



J3J4St(TT jlJ2St(7T 



(34) 



Here a^ aT creates a nucleon in a state (£j) with spin z- 
component a and isospin z-component r while A- star 
creates a pair of nucleons in the state \j1j2staT}. The 
operators A^ are straightforwardly obtained by vector- 
coupling products of two operators aL (TT and are not 
given explicitly. 

The eigenvalues and eigenfunctions of Ti are identified 
with the low-lying states of nuclei. For instance, nuclei 
pertaining to the 2sl<i-shell (the "sd shell nuclei") have 
mass numbers 17 < A < 39, and the number of valence 
nucleons has the range 1 < m < 23. Mass numbers 
A = 16 and A = 40 correspond to the closed-shell nu- 
clei 16 (the Is- and the lp-shells are filled) and 40 Ca 
(the 2sld-shell is filled, too). Likewise there are 2/lp- 
shell nuclei etc. (We disregard the lp-shell nuclei as 
they yield too few spectroscopic data for a meaningful 
statistical analysis). According to the shell model, the 
Hamiltonian Ti determines the spectral properties of the 
low-lying states of all nuclei pertaining to the same shell, 
at least in principle. This claim is subject to a number of 
provisos, (i) To get good fits to the spectra of all nuclei in 
a major shell it may be necessary to allow for a weak de- 
pendence of the parameters ehj and v a on the number m 
of valence nucleons. (ii) Non-valence-shell states may be 
pushed down by the residual interaction into the domain 
of low excitation energies ( "intruder states" ) and require 
special treatment. Such states are obtained, for instance, 
by lifting one or several nucleons from the valence shell 
into the next higher major shell, or from the inert core 
into the valence shell, or both, (iii) In its upper part, 
the spectrum of Ti cannot be expected to correspond to 
reality because the much more numerous non-valence- 
shell states dominate the actual spectrum and mix with 
the states in the valence shell. In the studies of chaos re- 
ported below, one disregards this fact and confines atten- 
tion to the valence shell. This is done in order to obtain a 
manageable numerical problem. There are strong reasons 
to believe that the results are universal and, thus, also 
hold when the mixing between major shells is taken into 
account (see, for instance, Ormand and Broglia, 1992). 
(iv) In the form of Eq. (|3"4"|) , the shell model applies to 
spherical nuclei (although it can be extended to weakly 
deformed nuclei). Chaos in deformed nuclei is treated in 
Section MM 

For purposes of orientation, we cite a few numbers 
taken from Bohr and Mottelson (1969) and Zelevin- 
sky et al. (1996). The spacing between major shells 
is approximately given by the harmonic-oscillator en- 
ergy Tito w 40 A" 1 / 3 MeV. The spin-orbit interaction is 
« —20 (l-s) A~ 2 / 3 MeV. The spacings of adjacent single- 
particle energies may be as large as a couple of MeV. In 
the sd-shell, for instance, the empirical values in 17 are 
e d 5/2 = -4.15 MeV, e s 1/2 = -3.28 MeV, e d 3/2 = 0.93 
MeV. With these figures and for m — 12 valence nucle- 
ons, the range of the elementary sd-shell spectrum would 
be about 42 MeV. (These figures are not completely rep- 
resentative since the egjS themselves are actually used 
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as fit parameters). The diagonal matrix elements of the 
residual interaction may be as large as 1 or 2 MeV in mag- 
nitude. This removes some of the degeneracies of the el- 
ementary shell model referred to above and stretches the 
spectrum further. The non-diagonal elements typically 
amount to several hundred keV in magnitude. These 
numbers suggest that the residual interaction is able to 
mix the states in different subshells of the valence shell, 
while the spectral properties at low excitation energy are 
essentially determined by valence-shell states (the ad- 
mixtures of non-valence-shell states are negligible). 

Not all parameters of the shell model can be de- 
termined equally well by a fit to spectroscopic data. 
That topic is discussed further in Sections IIV.A.2I and 
IV.B.3I below. To overcome that difficulty, one uses nu- 
clear many-body theory (i.e., variants of the Bcthc- 
Brueckner-Goldstone expansion) to calculate the two- 
body matrix elements of V res . The results serve as start- 
ing values for a fit to the data. In the fit, the two-body 
matrix elements themselves or some parameters on which 
they depend, are varied (Brown and Wildenthal, 1988; 
Honma et at, 2002). There is evidence that the nec- 
essary corrections mainly account for three-body forces 
(Caurier, 2005). 

Since TC conserves spin, isospin, and parity, the eigen- 
functions of TL are simultaneously eigenfunctions of total 
spin J, total isospin T, and parity n. One way to obtain 
the eigenfunctions of TL with given J, T, n for m valence 
nucleons consists in constructing the matrix of TL in a 
basis of m-body states carrying these quantum numbers 
and in diagonalizing that matrix. (This is not always the 
most efficient procedure numerically, but is used here for 
the sake of the argument). The m-body states needed 
for this procedure are obtained from the elementary shell 
model by constructing all Slater determinants of m va- 
lence nucleons in the valence shell. These determinants 
form classes, each class being defined by the set {mij} 
where mij is the number of nucleons occupying the sub- 
shell (£j). The set {mij} obviously forms a partition 
of m so that y\- mij = m where the sum runs over 
the subshells of the valence shell. The Slater determi- 
nants are antisymmetric by construction, but typically 
are not eigenstates of J and T. The m-body states with 
good J, T, and II and fixed {mij} are found as linear 
combinations of the determinants in class {mij} and are 
denoted by | JTII/i) in the sequel. (We suppress the mag- 
netic quantum numbers and the class index mij). The 
states | JTII/i) span a Hilbert space of finite dimension 
£>(J, T, II), which defines the range of the running in- 
dex (i. In the 2sld-shell and depending on the quantum 
numbers (J,T), T> ranges from a few to about 7000 in 
the middle of the shell (m = 12), while T> attains consid- 
erably larger maximum values already in the 2pl/-shell. 
This is why the 2sl<i-shell has been the preferred object 
for theoretical studies of chaos in nuclei. The actual con- 
struction of the states | JTII/i) is cumbersome, involves 
angular-momentum algebra, and is not given here (De 
Shalit and Talmi, 1963). In that basis, the matrix ele- 



ments of TL have the form 

def 



TL^JTU) = (JTUn\H\JTUiy) = S^^Etjmtj 



U 



22v a C flv (a; JTIl) 



(35) 



The first term on the right-hand side of the last of 
Eqs. (f35|) is obvious. The form of the second term fol- 
lows from Eq. (|34[) except that we have grouped together 
all matrix elements which are connected by symmetry. 
Save for this operation, the coefficients C fll/ (a; JTU) are 
matrix elements of the operator £ mT A] aHStmT A Jij2StmT 
taken between the states | JTII/i) and \JTTh>). By con- 
struction, the matrix TL^ V is real and symmetric. 

The form (j3"5|) displays explicitly the dependence of 
the matrix TL^, U {JTI\) on the input parameters Eij and 
v a of the shell model. We emphasize that the coefficients 
C^„(a; JTII) (which for fixed a, J, T, and II form a 
real and symmetric matrix in Hilbert space) are deter- 
mined entirely by the valence shell we are working in, 
by the coupling scheme we have used to construct the 
states | JTII/i), and by the two-body operator labeled a 
of which the matrix elements are taken. Except for a 
set of unitary transformations connecting the coupling 
scheme we have chosen with any other one, the matrices 
C t j tv {a) are uniquely determined. These matrices reflect 
the symmetries and invariances of the elementary shell 
model and are independent of the residual interaction ac- 
tually considered. In other words: Going from one resid- 
ual interaction with matrix elements v a to another one 
with matrix elements v' a , all we need to do is to replace 
the coefficients v a in Eq. ([33]) by the coefficients v' a , the 
matrices C^ u (a) remaining the same. We emphasize this 
simple fact because some properties of the shell-model 
Hamiltonian H^(JTH) are determined by the matrices 
Cn V {a) alone and are thus generic (i.e., largely indepen- 
dent of the choice of the residual interaction) . 

It was mentioned before and we emphasize again that 
the shell model accounts successfully for a vast amount of 
spectroscopic data in the ground-state domain of spher- 
ical nuclei. This is not the place to go into any de- 
tails. Suffice it to say that some basic features of V I0S 
are well established: Pairs of nucleons coupled to isospin 
T = 1 have a strong and attractive interaction ("pair- 
ing force"). That force leads to spin-zero ground states 
in even-even nuclei and favors the seniority coupling 
scheme (De Shalit and Talmi, 1963). In the particle-hole 
channel, the diagonal elements of V Tes for pairs of nucle- 
ons with isospin T = (neutron-proton pairs) also show 
strong attraction, especially for large angular momenta 
( "quadrupole-quadrupolc interaction" ) and favor nuclear 
deformations, especially for nuclei far (in mass num- 
ber) from closed-shell nuclei (see Bohr and Mottelson, 
1969). For the sd-shell, Wildenthal (1984) and Brown 
and Wildenthal (1988) have established the standard pa- 
rameters of the residual interaction ( "Brown- Wildenthal 
interaction") by fitting the 66 parameters (63 two-body 
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matrix elements and 3 single-particle energies) to more 
than 400 pieces of data. The Brown- Wildcnthal inter- 
action may optionally include the Coulomb interaction 
between protons. 



2. Chaos in the Shell Model 

Ever since the shell model was established, a strange 
dichotomy pervaded nuclear-structure theory. On the 
one hand, the shell model was extremely successful in de- 
scribing the properties of low-lying states in many nuclei, 
and the collective models in their various forms afforded 
a description of those nuclei which were not accessible to 
the shell model. Fits to data revealed the basic proper- 
ties of the residual interaction. Ever more sophisticated 
measurements widened the data basis. Theoretical ef- 
forts were directed at both the determination of Vi- es from 
the interaction between free nucleons via many-body the- 
ory, and at the technology to calculate nuclear properties 
from the shell-model Hamiltonian, including those of de- 
formed nuclei, and thus at understanding collective mod- 
els. The hope seemed justified that one day one would 
arrive at a complete understanding of the structure of 
atomic nuclei. In the 1970s and 1980s, most work in 
nuclear-structure theory was devoted to that vision. 

On the other hand, during the same period the evi- 
dence for the applicability of RMT to neutron and pro- 
ton resonance data grew and was definitively established 
in 1982, see Section lilLC.il But aside from fundamental 
symmetries, RMT lacks any dynamical structure what- 
soever. What did the success of RMT imply for the shell 
model? Is there an excitation energy (somewhere below 
neutron threshold but above the energies where shell- 
model calculations were so successful) beyond which the 
shell model fails to work and chaos takes over? Or is 
chaos possible (and perhaps even generic) in the frame- 
work of the shell model itself? 

In spite of a strong growth of chaos-related research 
in other fields of physics in the 1980s and early 1990s, 
these questions did not receive much attention at the 
time by the nuclear-structure community. The main 
reason was probably the lack of statistically significant 
experimental data. In addition, a strong commitment 
on the part of the community to understand nuclei on 
a fundamental level perhaps deflected attention away 
from the issue of chaos. Only a few theoretical pa- 
pers addressed that issue. Early work by Whitehead et 
al. (1978), Verbaarschot and Brussard (1979), Brown 
and Bertsch (1984) and Dias et al. (1989) addressed 
the validity of the Porter-Thomas distribution for shell- 
model eigenfunctions, see the comments below. The re- 
lation between spectral fluctuations and the shell model 
was addressed by Weidenmujllcr (1985). Meredith et 
al. (1988) and Meredith (1993) studied the Lipkin- 
Meshkov-Glick model: M. Fermions occupy a system 
with three non-degenerate subshells each containing M. 
degenerate single-particle states. In each subshell, the 



single-particle states do not carry any further quantum 
numbers. The two-body interaction acts only between 
particles in different subshells. All non-zero matrix el- 
ements are identical. The Hilbert space has finite di- 
mension. Symmetries reduce the problem to manage- 
able size even when M. is large. The system possesses a 
classical limit which is attained for coherent states when 
M. — > oo. Classical chaos can thus be studied in its de- 
pendence on the strength of the two-body interaction. 
As that strength is varied, a close correspondence is nu- 
merically established between the transition from regular 
to chaotic motion and that from Poisson to GOE statis- 
tics for the spectral fluctuation measures of the quantum 
system. The Bohigas-Giannoni-Schmit conjecture was 
verified for the first time for an interacting many-body 
system (albeit in the framework of a toy model without 
characteristic ingredients of the shell model such as con- 
served quantum numbers). In the 1990s, several papers 
(Aberg, 1990; Alhassid et al, 1990; Alhassid and Whelan 
(1991); Martinez-Pinedo et al. (1997)) addressed chaos 
in the collective model. This topic is reviewed in Sec- 
tion WM 

Ormand and Broglia (1992) reported a study of quan- 
tum chaos in the shell model for the sd-shell. That work 
displayed several features which were also discussed in 
the paper by Zelevinsky et al. (1996). We focus on 
that later, much more extensive paper. The authors un- 
dertook a thorough and systematic theoretical study of 
chaos in the sd-shell. They used the well-established 
Brown-Wildenthal (1988) form of the residual interac- 
tion, which is free of random elements. In most of their 
calculations Coulomb effects were neglected. They fo- 
cused attention mainly on nuclei in the middle of the 
shell (to = 12) where the dimensions T>{ J, T) of the shell- 
model matrices are manageable and yet sufficiently large 
for a statistical analysis. (We drop the label II since 
all states in the sd-shell have positive parity). The au- 
thors calculated spectra and eigenfunctions of the shell- 
model Hamiltonian numerically and compared the results 
both with GOE predictions and with thermal averages. 
The latter were considered because chaos was suspected 
by Percival (1973) to cause the many-body eigenfunc- 
tions all "to look the same" , so that concepts such as 
mean occupation numbers for single-particle states and 
the Fermi-Dirac distribution may make sense for observ- 
ables averaged locally over energy. 

The authors used ergodicity (Section III. C.3[) to com- 
pare GOE ensemble averages with running averages over 
the calculated spectra. To obtain statistically significant 
results, they used all eigenvalues and/or eigenfunctions 
pertaining to fixed m and to a pair (J, T) of quantum 
numbers. As pointed out before (Section IIV.A.1[) . cal- 
culations using only states from the valence shell, yield 
results which are unrealistic in the upper part of the spec- 
trum where non-valence-shell states actually dominate. 
Therefore the work of Zelevinsky et al. (1996) must be 
considered a case study of quantum chaos within the shell 
model in a restricted Hilbert space, rather than a realis- 
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tic calculation of nuclear properties. While we summarize 
the results of this work in the present Section, we defer 
a detailed analysis of how chaos is generated in the shell 
model to Section IV.BI 

The average level density (calculated from the actual 
spectrum by smoothing) differs from that of the GOE 
and has approximately Gaussian shape. That result was 
anticipated many years ago by Mon and French (1975), 
see also Section IV.AI The nearly Gaussian shape is an 
artifact, of course, and due to the restriction to a single 
major shell; the actual nuclear level density grows nearly 
exponentially with excitation energy. Results of primary 
interest in the present context relate to quantum chaos 
and concern the local spectral fluctuation measures. Af- 
ter unfolding of the spectra, the NNS distribution and the 
A 3 -statistic were found to agree well with GOE predic- 
tions in the middle of the shell and for sufficiently large 
values of the matrix dimension D(J,T). For m = 12 
and the states with J = 2, T = (X>(2, 0) = 3276) this is 
shown in Figs.fT9land[20l It is of interest to follow the on- 
set of chaos as the strength of V^ e s is varied. For V Tes = 
the NNS distribution is found to be close to Poissonian as 
expected. As the strength is increased, the distribution 
approaches the Wigner distribution and, within statis- 
tics, becomes indistinguishable from that distribution al- 
ready for strength values amounting to 20 % of the ac- 
tual one. For smaller values of m ( 24 Mg with m = 8) the 
NNS distribution at full strength value is still interme- 
diate between the Poisson and the Wigner distribution, 
even for sizable values of D ( J, T) (D(0, 0) = 1161). Here, 
inclusion of the Coulomb interaction reduces the proba- 
bility of very small spacings. The A3-statistic in Fig. [20] 
agrees with the GOE prediction except for large values of 
L where it rises above that prediction. That seems to be 
a systematic trend and is not fully understood at present. 
It perhaps signals a weakening of GOE-type correlations 
between eigenvalues with spacings of more than 7 MeV 
or so and might suggest that strong mixing is restricted 
to unperturbed shell-model configurations with energies 
within an energy interval of 10 to 20 MeV. We return to 
this point below. (A later evaluation of A^(L) for the 3~ 
states showed good agreement with the GOE prediction 
op to L ss 3000, however, see Zelevinsky (2007).) 

Extensive shell-model calculations have also been done 
for nuclei in the 2pl /-shell, especially by the Strasbourg- 
Madrid collaboration (Caurier et al., 2005). We are 
not aware, however, of an equally thorough analysis of 
the results with respect to spectral fluctuation measures 
and other indicators for quantum chaos as done for the 
sd-shell by Zelevinsky et al. (1996). Kota (2001) has 
analyzed the calculated spectrum of 4 + states in 48 Ca 
(Martinez-Pinedo et al., 1997; Caurier et al., 1999.) Us- 
ing the 1355 states located in the middle of the spectrum 
(out of a total of 1755 levels), he finds good agreement 
with the Wigner surmise for the NNS distribution, and 
with the GOE prediction for the A3-statistic. The latter 
is reported only for L < 60. Electromagnetic transition 
intensities and moments for A rs 60 nuclei as calculated 




FIG. 19 NNS distribution (histogram) and Wigner surmise 
(solid line) for the states with m = 12 and J = 2, T — in 
the sd-shell. From Zelevinsky et al. (1996). 



by Hamoudi et al. (2002) showed good agreement with 
the Porter-Thomas distribution. 

Investigation of the distribution of eigenfunctions of- 
fers additional insight and shows the limitations of spec- 
tral fluctuation measures in establishing chaos. Here the 
Porter-Thomas distribution is the standard measure. In 
early calculations (Whitehead et al., 1978; Verbaarschot 
and Brussard, 1979; Brown and Bertsch, 1984) for the 
sd-shell, it was found that the actual distribution of 
widths differed from the Porter-Thomas form: There 
was an overabundance of the largest and the smallest 
widths at the expense of those with values closer to the 
median. This was ascribed to incomplete mixing of the 
configurations especially in the wings of the unperturbed 
spectrum: Too little mixing leads to admixtures with too 
small amplitudes and leaves the original states too pure. 
The issue was followed up in more detail by Zelevinsky 
et al. (1996). We denote by Wk^ the square of the am- 
plitude with which an unperturbed shell-model config- 
uration labeled fi is admixed into an eigenstate of the 
shell-model Hamiltonian. The index k labels the eigen- 
states such that the associated eigenvalues Ek increase 
monotonically with k. The measure used was the "infor- 
mation entropy" Su, which is defined as 



s k = -J2 w ^ lnW ^ ■ 



(36) 



For complete mixing and large matrix dimension T> 3> 1 , 
wave-function normalization makes us expect W-kp, ~ 
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FIG. 20 The A 3 -statistic (dots) and the GOE prediction 
(solid line) for the states with m = 12 and J = 2, T = 
in the sd-shell. From Zelevinsky et al. (1996). 
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1/2? independently of k and //: On average, all configura- 
tions are equally mixed into every eigenstate. This corre- 
sponds to a maximum value of Sk = In 2?. That naive es- 
timate does not take into account the Porter-Thomas dis- 
tribution for the Wku,s. Doing so yields Sk = ln(0.482?) 
for the maximum value of Sk (Zelevinsky et al., 1996). 
We refer to that value as to the GOE limit. Incomplete 
mixing is bound to reduce that value. 

The measure (j36|) has the advantage that it allows for a 
detailed study of configurational mixing versus the eigen- 
function index k and thus offers more insight than af- 
forded by the Porter-Thomas distribution. This advan- 
tage has to be weighed against the shortcoming that the 
measure depends on the representation chosen for calcu- 
lating the Wfc^s. (Indeed, Sk is not invariant under or- 
thogonal transformations of the basis of states {/i}, while 
the Porter-Thomas distribution applies to every projec- 
tion of the eigenvectors of the GOE onto a fixed vector in 
Hilbert space). The "natural" representation chosen by 
Zelevinsky et al. (1996) for Sk is defined by the unper- 
turbed shell-model configurations. The implications and 
limitations of that choice are discussed by Zelevinsky et 
al. (1996). 

Fig. [5T] shows the exponential of Sk versus k for the 
same 3276 J = 2, T = states as used in Figs. []jj] and 
[201 The GOE limit (1578) is almost reached in the middle 
of the spectrum. In the wings, the mixing of states is 
fairly incomplete. That pattern is generic and reinforces 
earlier findings (Whitehead et al, 1978; Verbaarschot and 
Brussard, 1979; Brown and Bertsch, 1984; Dias et al., 
1989). 

The eigenfunctions of the shell-model Hamiltonian are 
orthogonal and normalized. In the GOE, these condi- 
tions impose weak correlations among the expansion co- 



FIG. 21 The exponential of the information entropy (|36|1 ver- 
sus the energy of the state for the states with m = 12, J = 2 
and T = in the sd-shell (dots). The GOE limit of 1578 
is indicated by the horizontal line. From Zelevinsky et al. 
(1996). 



efficients of the eigenfunctions in an arbitrary basis which 
disappear for N — * oo. The correlations for the expansion 
coefficients of the shell-model Hamiltonian in the basis 
of unperturbed shell-model configurations are found to 
be somewhat larger than those of a GOE with the same 
dimension, especially in the tails of the spectrum and for 
the 10 or 20 nearest eigenstates. 

In Section III.GI we introduced the concept of the 
strength function for a doorway state. In analogy to 
Eq. (|25p , the strength function for an unperturbed shell- 
model configuration /i with respect to the exact eigen- 
states of the shell-model Hamiltonian is defined by 
J2k WknS(E — Ek). Since no ensemble averaging is in- 
volved in that definition, it is useful to average over 
a group of neighboring (in energy) shell-model config- 
urations to get a smooth function. If the eigenstates 
were complete mixtures of shell-model configurations, 
the smoothed strength function would be constant (in- 
dependent of energy E). In actual fact, the strength 
functions are peaked with a full width at half maximum 
of about 20 MeV even for shell-model configurations k 
in the middle of the unperturbed spectrum. The shape 
of the strength function is Gaussian near its peak, but 
only falls off with an exponential tail in the wings. Such 
exponential decay has been studied by Lewenkopf and 
Zelevinsky (1994) and Frazier et al. (1996). 
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Another measure which shows that the eigenfunctions 
of the shell-model Hamiltonian are not perfect mixtures 
of the unperturbed shell-model configurations, is pro- 
vided by the pairing force, an attractive interaction be- 
tween like nucleons. That force corresponds to a par- 
ticular linear combination of the interaction operators 
Y,ar A huat<TT A hnst<TT appearing in Eq. J33J). For a com- 
pletely mixed system, a plot of the expectation values of 
that linear combination (taken with respect to the eigen- 
states of the shell model Hamiltonian) versus /i is ex- 
pected to fluctuate about a constant mean value. In 
actual fact, the plot shows a systematic enhancement 
(with respect to the mean value) of about 70 % in the 
ground-state domain at the expense of a corresponding 
suppression at the upper end of the spectrum, with small 
fluctuations. 

The complete mixing of the eigenfunctions due to chaos 
may be similar to thermalization. As a test, Zelevinsky 
et al. (1996) calculated the occupation numbers of the 
three single-particle states efe/2, Si/2, and d 3 / 2 in the ex- 
act eigenstates, and plotted the result versus /i. The data 
have small fluctuations. With the help of a properly de- 
fined temperature, the results can be very well fitted with 
the Fermi distribution. The fit parameters are "effective" 
single-particle energies. These differ by only a few 100 
keV from the input parameters (the single-particle ener- 
gies eij of the shell model) . This corroborates the picture 
of thermalization. 

In summary: If chaos is measured in terms of the usual 
spectral fluctuation measures (NNS distribution and A3- 
statistic), there is evidence that the residual interaction 
mixes the unperturbed shell-model configurations suffi- 
ciently strongly to produce chaos. This is true for nuclei 
with sufficiently many valence nucleons and in the mid- 
dle of the spectrum. Chaos is diminished for nuclei with 
a smaller number of valence nucleons (or of holes in the 
valence shell) . Closer inspection of the eigenvectors offers 
a more subtle picture. Even for nuclei in the middle of 
the shell, the mixing of the unperturbed shell-model con- 
figurations is not complete, especially in the tails of the 
spectrum. Small correlations between eigenvector com- 
ponents exist beyond the ones imposed by orthogonality 
and normalization. The strength functions of the unper- 
turbed shell-model configurations are not constant, but 
are peaked with a width of twenty MeV or so. The ex- 
pectation values of the operator of the pairing force are 
enhanced in the ground-state domain. These deviations 
from GOE properties may be related to the fact that the 
A3-statistic shows an upward bend for large values of L. 
All this shows that within a single major shell, chaos is 
not fully attained for realistic strengths of the interaction, 
although it would be difficult to detect such deviations 
experimentally (except for the behavior of A3). In spite 
of the deviations, thermodynamic concepts apply, and 
the single-particle occupation numbers follow the Fermi 
distribution function. 

We return to the questions posed at the beginning of 
this Section. From the evidence presented, it seems likely 



that calculations that would allow for the presence of 
non-valence-shell states would show that as the excita- 
tion energy increases, nuclear levels attain an ever greater 
similarity to GOE eigenstates. The approach to the GOE 
limit is probably somewhat faster for nuclei in the mid- 
dle between two closed shells than for nuclei near closed 
shells. Chaos thus seems a natural ingredient of the shell 
model, and not an exclusive alternative to regular mo- 
tion as seen in the ground-state domain. The results on 
the strength function for the unperturbed shell-model 
configurations and on the A3~statistic suggest that at 
excitation energy E the strong mixing (which is char- 
acteristic of quantum chaos) involves only those unper- 
turbed shell-model configurations which are located in 
an energy interval centered at E and 10 to 20 MeV wide. 

The very success of the shell model, i.e., the ability 
of the model to account for many basic features of nu- 
clei, implies that many-body states pertaining to differ- 
ent major shells are mixed only weakly. In that sense 
the residual interaction is weak: It removes the numer- 
ous degeneracies of the single-particle model in a manner 
that causes nearly complete mixing of states within a ma- 
jor shell. However, the residual interaction is not strong 
enough to destroy the overall shell structure defined by 
the existence of major shells. In that sense, nuclei are 
not fully chaotic systems. This point of view has been 
emphasized by Bunakov (1999). 

Chaos does not preclude the existence of regular fea- 
tures. These are seen in the ground-state domain where 
the mixing of unperturbed shell-model configurations is 
incomplete, and in the presence of collective modes of 
excitation which act as doorway states and are seen as 
giant resonances. Chaos is helpful because it allows for 
the description of average properties of excited nuclei in 
terms of concepts of equilibrium statistical mechanics. It 
remains to show by which mechanism chaos originates in 
the shell model, and to clarify whether chaos is a generic 
feature of the shell model or depends upon specific prop- 
erties of the residual interaction. We return to these ques- 
tions in Section IV. Bl 

An analysis similar to the one in sd-shell nuclei de- 
scribed above was carried out for the Ce atom by Flam- 
baum et al. (1994) with very similar results. This sup- 
ports the view that chaos is a generic property of self- 
bound many-body systems. 



3. Limits of Validity of RMT in Nuclei 

No real physical system can be expected to possess 
spectral fluctuation properties that coincide exactly with 
RMT predictions. Indeed, RMT is based upon a purely 
stochastic approach, and we must expect that at some 
point, system-specific features dominate spectral prop- 
erties. Which then are the limitations of RMT in nuclei? 

The answer to that type of question is known in sys- 
tems with few degrees of freedom (Berry, 1985) and in 
disordered systems (Imry, 2002). For chaotic systems 
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with few degrees of freedom, the limitations of RMT 
are connected to the shortest periodic orbit in the sys- 
tem. With r m i n the period of the shortest periodic or- 
bit, AE = h/T m i n defines the maximum energy interval 
within which RMT predictions can be expected to hold. 
In disordered systems, the period T m ; n is replaced by the 
diffusion time Tdiff, and the characteristic energy inter- 
val is given by the Thouless energy E c — h/rdis- The 
dimensionless Thouless conductance g is the ratio of ei- 
ther of these intervals and the average level spacing and 
gives the number of levels over which RMT predictions 
are expected to hold. 

In nuclei, the situation is not completely clear, and two 
different schools of thought exist. Bohigas and Leboeuf 
(2002) used a mean-field approach. It is argued that 
the mean-field motion is partly chaotic. The shortest 
periodic orbit at the Fermi energy is used to estimate 
the characteristic energy interval as AE = 77.5A~ 1//3 
MeV. The approach was worked out further by Olofs- 
son et al. (2006). However, the mean-field approach 
addresses single-particle properties only. It is relevant 
for one-body chaos but is not clearly related to chaos in 
a many-body system: Without two- or many-body in- 
teraction the eigenvalues of the many-body system are 
sums of single-particle energies and have a Poisson distri- 
bution irrespective of whether the single-particle motion 
is regular or chaotic. 

A different view was taken, for instance, by Bunakov 
(1999) and Molinari and Weidenmiiller (2006), who argue 
that the independent-particle model gives rise to regular 
motion, while the residual interaction causes mixing of 
shell-model configurations and, thus, chaos. The char- 
acteristic energy interval over which configurations are 
strongly mixed (and the range of energies over which 
RMT predictions apply) is then given by the spreading 
width. In Section IIV.A.2I it was shown that in shell- 
model calculations, the spreading width is found to be 
of the order of 10 MeV. While that number is not sub- 
stantially different from the estimate obtained within the 
chaotic mean-field scenario, the origins of both estimates 
are clearly very different. A test of these predictions is 
not possible using experimental data. As pointed out 
before, sufficiently long sequences of levels with identi- 
cal quantum numbers are not known. The shell-model 
calculations reported on in Section IIV.A.2I seem to lend 
substance to the second view. 



B. Collective Models 

1. The Collective Model of Bohr and Mottelson 

Nuclei can undergo shape deformations. This fact be- 
came obvious with the discovery of nuclear fission in 1939: 
A very heavy nucleus splits spontaneously into two frag- 
ments of about equal mass. The energy liberated in that 
process could be roughly understood on the basis of the 
liquid-drop model. In that model, nuclei are described as 



charged droplets of an incompressible fluid held together 
by surface tension. The latter mimics the attractive nu- 
clear force. 

A dynamical theory of surface deformations was devel- 
oped by Bohr (1951 and 1952) and Bohr and Mottelson 
(1952) in the early 1950s. In an expansion of the shape 
of the nuclear surface in terms of spherical harmonics 
Yip, attention is focused on the lowest (in £) non-trivial 
terms. These are the quadrupole terms Y^ (the term 
with £ = is ruled out because of volume conserva- 
tion and the terms with £ = 1 describe the motion of 
the center of mass). For small values of the five expan- 
sion parameters a-i^ with /i = —2, — 1, 0, +1, +2, the sur- 
face has the shape of an ellipsoid. The principal axes 
of that ellipsoid define an intrinsic ( "body-fixed" ) coor- 
dinate system. The five expansion parameters a^ can 
be transformed into the three Euler angles which spec- 
ify the location of the body-fixed system with respect to 
the laboratory system, and two parameters (commonly 
called (3 and 7) which specify the nuclear shape in the 
body-fixed system. The potential energy V(/3, 7) defines 
the static energy of quadrupolar nuclear shape deforma- 
tions. A dynamical theory is obtained by considering the 
parameters a^ (or their transforms) as dynamical vari- 
ables which obey bosonic commutation relations. The re- 
sulting Hamiltonian ( "Bohr Hamiltonian" ) has a number 
of parameters (the "masses" connected with the kinetic 
energy terms and the parameters specifying the poten- 
tial energy V(f3,jj). Depending on the values of these 
parameters, the theory predicts inter alias rotational 
motion and vibrational motion of nuclei. In even-even 
nuclei, rotational motion manifests itself in the occur- 
rence of rotational bands (sequences of levels with spins 
J = 0, 2, 4, . . .). The excitation energies (above the band 
head) of the states with spin J in the band are propor- 
tional to J (J +1). The eigenf unctions are obtained by 
projecting the deformed intrinsic state onto fixed angular 
momentum J. This is done using the Wigner T> functions 
and integrating over Euler angles. The electromagnetic 
transitions within the band are electric quadrupole (E2) 
and are strongly enhanced over simple single-particle 
(i.e., shell-model) estimates. The transition matrix ele- 
ments are proportional to the static quadrupole moment 
of the intrinsic deformed state. In Fig. [5H two such ro- 
tational bands are displayed. Vibrational motion mani- 
fests itself in (nearly) harmonic vibrations of the surface 
about its equilibrium shape and is characterized by an 
harmonic-oscillator-like spectrum. The electromagnetic 
transition matrix elements are also enhanced over single- 
particle estimates but not as much as in the rotational 
case. This "collective model" (so named because many 
nucleons partake in an orderly way in the motion) of Bohr 
and Mottelson has been extremely successful in account- 
ing for many spectroscopic data (Bohr and Mottelson, 
1969). In cases of pure rotational and pure vibrational 
motion we expect, of course, a predominance of regular 
features. This view is supported by the empirical evi- 
dence reviewed in Sections IIII.C.2I and IIII.C.31 
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FIG. 22 Two rotational bands in 174 Hf. Adapted from Bohr 
and Mottelson (1969) Vol. 2, p. 168. 



deformed single-particle potential. The parameters of 
that potential are obtained by fits to the data. (We men- 
tion in passing that deformed single-particle potentials 
with a pure quadrupole deformation do not give rise to 
chaotic single-particle motion. An additional octupole 
deformation is needed. Even then, chaos arises only in 
the oblate case (Arvieu et at, 1987; Heiss et al., 1994). 

The obvious difficulty with this approach is that by 
construction, i?Nilsson is not rotationally invariant. This 
is not a problem in principle. We recall that according 
to the collective model (see Section HV.B.ip each rota- 
tional band is due to a deformed "intrinsic state" . The 
wave function of each band member (which has definite 
spin J) is obtained by projecting the intrinsic state onto 
spin J. In the same sense, the many-body eigenstates of 
-Hnussoii (Slater determinants) are viewed as microscopic 
realizations of intrinsic states. Projection onto states of 
good total spin will generate from each such Slater de- 
terminant the members of a rotational band. 

However, in practice projection of Slater determinants 
is cumbersome, and an approximation (the "cranking 
model") is used. It is postulated that in the laboratory 
system, the deformed single-particle potential rotates (is 
"cranked" ) about some axis with fixed frequency to. That 
axis must be perpendicular to the symmetry axis of the 
potential (since in quantum mechanics, rotation about a 
symmetry axis is not possible). Under a coordinate trans- 
formation from the laboratory to the body-fixed frame of 
reference, the rotating deformed single-particle potential 
becomes a static deformed potential. As in classical me- 
chanics, that coordinate transformation induces an addi- 
tional term in the Hamiltonian, and the single-particle 
Hamiltonian of the cranking model is 



H, 



cranking 



= H 



Nilsson 



WJ 



(37) 



2. Onset of Chaos in Rapidly Rotating Nuclei 

Studies of chaos in the collective model (Aberg, 1990; 
Matsuo et al., 1997) have addressed the onset of chaos 
above the yrast line for states of high spin. (The yrast 
line is defined in Section IIII.C1 As a function of J, it 
is given by the energy of the lowest level with spin J). 
We review first the later, more extensive paper by Mat- 
suo et al. (1997) and then the earlier work by Aberg 
(1990). In both papers, very similar techniques were 
used and similar results were obtained. The starting 
point is a generalized shell model. Chaos in these ap- 
proaches is due to the residual interaction which mixes 
the basic shell-model configurations. The approach is 
thus similar to that of Section IIV.A.1I However, the 
shell-model used differs from the spherical shell model of 
Section H V.A.I I This "Nilsson model" takes into account 
the fact that rotational motion in medium-weight and 
heavy nuclei is due to deformations of the nuclear shape. 
Thus the single-particle Hamiltonian -ffNiisson of the Nils- 
son model contains, in addition to the kinetic energy and 
a spin-orbit coupling term, a non-spherical, elliptically 



Here lu points in the direction of the axis of rotation, and 
j is the total spin of the nucleon. The ^-direction of the 
body-fixed system is commonly assumed to coincide with 
the symmetry axis, and the direction of lu is assumed to 
coincide with the x-axis, so that the last term in Eq. (I37[) 
(comprising both Coriolis and centrifugal forces) takes 
the form uij x - 

The many-body solutions of the cranking Hamilto- 
nian are Slater determinants. Let |0) be the vacuum 
state, and let a\ be the creation operator for the cranked 
single-particle state labeled \i{w)). Then, \fi) = J^^ aj |0) 
is a Slater determinant of cranked states (all taken at 
the same frequency lu). The label y, represents a set of 
occupied orbitals. Different choices of /i correspond to 
the ground state and to the excited states of the crank- 
ing model. It is assumed that the deformed potential 
is the same for all of these states. This assumption is 
realistic up to excitation energies of several MeV for nu- 
clei for which the potential energy of deformation dis- 
plays a deep and stable minimum. The states \\i) de- 
pend on the cranking frequency lu. Taken as functions 
of w, the single-particle energies of the cranking model 
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display avoided crossings. At such crossings, the adia- 
batic single-particle wave functions change abruptly. To 
obtain states |//) which depend smoothly on w, a diabatic 
single-particle basis is used by Aberg (1990) and Matsuo 
et al. (1997). 

To generate states of (approximate) total spin J from 
a given state |/i), the consistency condition (J x )(u) = J 
is used to determine to. Here J x is the x-component 
of the total spin operator (the sum of the spins of all 
nucleons). In other words, the cranking frequency u> is 
adjusted so that the system rotates on average with the 
desired spin J. In Matsuo et al. (1997) the average (J x ) 
is the thermal average over all cranked single-particle 
states taken with a temperature T = 0.4 MeV. This tem- 
perature corresponds to the excitation energies of inter- 
est (about 2 MeV above yrast). The states with different 
spins generated in this way form a rotational band, with 
the state \fj) playing the role of the intrinsic state. This 
is seen by using a Taylor expansion of the expectation 
value (/il-ffcrankinglA 4 ) in powers of w. That same expan- 
sion is used to transform the intrinsic energies back to 
the laboratory system. 

The states of fixed J generated in this fashion are or- 
thogonal. They are mixed by the residual interaction. 
In Matsuo et al., (1997) a two-body interaction of the 
surface delta type was used (the interaction is confined 
to an infinitely thin layer of the nuclear surface) . The in- 
teraction strength was determined by previous studies of 
rotational nuclei. The resulting Hamiltonian contains the 
eigenvalues of the cranking model in the laboratory sys- 
tem as diagonal elements and the matrix elements of the 
residual interaction. To obtain a manageable problem, 
only the lowest 1000 states |/i) were used. The result- 
ing lowest 300 eigenstates of the total Hamiltonian were 
found to be rather stable against that truncation. These 
are used in the statistical analysis. They cover a region 
of excitation energy up to about 2.4 MeV above yrast. 

After unfolding, the spectra were binned. The NNS 
distribution and the A3-statistic show a gradual transi- 
tion from near Poissonian behavior in the lowest bin to 
near GOE behavior in the highest one. For the NNS dis- 
tribution, this is indicated by the Brody parameter (see 
Section IIII.A|) , which increases monotonically with exci- 
tation energy U and reaches values close to unity at the 
upper end of the spectrum. The dependence of the Brody 
parameter on U is the same for all spin values studied. 
The A3~statistic agrees with the GOE value only up to a 
maximum value of L. That value increases with increas- 
ing U. Even for the bin containing the 50 states with 
highest excitation energies, however, this maximum value 
is as small as 6. This shows that in the cranking model, 
spectral stiffness is a local phenomenon. We recall that 
similar features (although for much larger values of L) are 
found in the spherical shell model, see Section HV. A. 21 

Essentially the same model was used in Aberg (1990). 
However, rather than studying the mixing of rotational 
bands at fixed interaction strength as a function of excita- 
tion energy, the excitation energy was held fixed and the 



strength A of the interaction was varied. This procedure 
simulates the increase of level density with increasing ex- 
citation energy. The results for the A3-statistic were 
analyzed by writing A3 as the sum of the GOE expres- 
sion ([22]) (taken at a scaled length qL with < q < 1) 
and of a term linear in L (as for the Poisson distribu- 
tion) with relative weight (1 — q). The fit parameter q 
was determined as a function of the strength A of the 
two-body interaction. The result is shown in Fig. [23] 
The fragmentation of collective strength was also inves- 
tigated. The cranking frequency depends upon J and so 
do, therefore, the single-particle energies of the crank- 
ing model. As a result, states with different values of J 
are mixed differently by the residual interaction. This 
implies that the collective El transition strength (which 
for unmixed rotational bands only connects states within 
the same band) becomes fragmented. As a measure of the 
distribution of the reduced matrix elements for E2 decay 
connecting a mother state with spin J and the daughter 
states with spin J — 2, the standard deviation is used. An 
average value ~oei for the standard deviation is obtained 
by averaging over 50 mother states. Figure [221 also shows 
~5ei versus A, the strength of the two-body interaction. 
Also shown is 1/2 r rot (Aberg, 1990; Matsuo et al, 1997), 
the half width of the average distribution function of the 
matrix elements. That function changes from Gaussian 
to Brcit-Wigner form as A increases. This fact explains 
why as a function of A, 1/2 r rot displays a maximum. 
We observe that as functions of A, q and We2 behave 
very similarly. This shows that the fragmentation of the 
E2 transition strength is an indirect measure of chaos. 
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FIG. 23 Mixing parameter q for the A3-statistic (see text) 
(left-hand scale), average standard deviation 7je2 of the dis- 
tribution of reduced E2 matrix elements (right-hand scale), 
and half width at half maximum l/2r rot (right-hand scale) 
for the same distribution, all versus A, the strength of the 
two-body interaction which mixes the rotational bands. From 
Aberg (1990). 

The theoretical ideas of Aberg (1990) and Matsuo et al. 
(1997) were applied to data by Stephens et al. (2005). A 
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beam of 48 Ca ions of 215 MeV hit a target of 124 Sn. The 
resulting isotopes of Yb decay by gamma emission. Pairs 
of gamma quanta were measured in coincidence. Nu- 
merical simulations along the lines described above were 
compared with the measured spectra. The fits were used 
to determine the ratio of the strength A of the nucleon- 
nucleon interaction versus the mean level spacing. The 
ratio covers the range from 0.15 (nearly fully ordered) to 
1.5 (nearly fully chaotic), in an energy interval which is 
consistent with the theoretical work. 

Some authors (Aberg, 1992; Mottelson, 1992) have 
proposed that at excitation energies of a few MeV above 
the yrast line, the Coriolis term in Eq. (|3T[) may be rather 
weak. Then rotational bands would exist with completely 
mixed wave functions of the GOE type. Nevertheless, the 
collective E2 decay would take place entirely within each 
band. 

On the basis of a representation using both quasipar- 
ticles and phonons, order, chaos, and the order-to-chaos 
transition are discussed by Soloviev (1995). 



3. Chaos in the Interacting Boson Model 

The shell model is generally considered the fundamen- 
tal phenomenological nuclear-structure model (Caurier 
et a/., 2005). However, applications of the model have 
been restricted in practice to nuclei for which the num- 
ber of valence nucleons is not too large. This is the case 
for all nuclei in the lp-shell and the 2slri-shell and, more 
recently, for most nuclei in the 2pl /-shell. For yet heav- 
ier nuclei (and thus shells beyond the 2pl/-shell), the 
number of valence nucleons for nuclei near the middle 
of the shells is simply too large, and a shell-model cal- 
culation is prohibitively difficult. At the same time, it 
is in these mass regions that the collective model finds 
its most successful application. This fact has prompted 
many authors to look for a derivation of the collective 
model from the shell model by introducing suitable col- 
lective variables. We mention, in particular, the idea of 
using a boson expansion within the shell model to ob- 
tain a simplified description with built-in collective fea- 
tures. In this approach, pairs of nucleons form boson-like 
entities. None of these approaches was truly successful 
in providing a derivation of the Bohr Hamiltonian from 
the shell model within controlled approximations. How- 
ever, in the midst of these efforts a phenomenological 
approach emerged which has become eminently success- 
ful, the interacting boson model (IBM) of Arima and 
Iachello (1975). The IBM postulates the existence of s- 
bosons and ri-bosons. These bosons may be thought of 
as representing pairs of nucleons coupled to spin and 2, 
respectively. Alternatively, the ci-bosons may be viewed 
as the five quanta of quadrupole surface deformations. 
The s-boson is then an artifice which is used to simplify 
the mathematics. For each nucleus, the total number N 
of bosons is fixed in the IBM. In the limit N — * oo, the 
solutions of the boson Hamiltonian approach those of the 



collective model (Dieperinck et at, 1980). The connec- 
tion between the IBM and the shell model is discussed by 
Iachello and Talmi (1987). The IBM at large is reviewed 
by Iachello and Arima (1987). 

In the original form of the IBM (often referred to as 
IBM 1) no distinction is made between nucleon pairs 
formed of protons and pairs formed of neutrons. This 
is the form of the IBM for which an analysis with re- 
gard to chaotic properties has been performed (Alhassid 
et al., 1990; Alhassid and Whelan, 1991). To introduce 
the model, we defined s^ and dt as the creation opera- 
tors for the s-boson and the five ri-bosons, respectively. 
While s> is a scalar, the d], transform under rotations as 
the components of an irreducible tensor of rank 2. The 
modified annihilation operators s = s and cL = (— ) M cL 
have the same transformation properties. The symbol 
(tff x d) K denotes the irreducible tensor of rank k with 
spherical components k obtained by vector-coupling <v 
and d, and likewise for s^ and s. A quadrupole operator 
Q(X) is defined by Q( X ) = (d f x S+s t x d)^ + x {^ x d)^ 
where x is a real parameter. Coupling Q{x) with itself 
to an irreducible tensor of rank zero generates the scalar 
Q(x) ■ Q(x)- I n the spherical representation, the three 
components of the angular-momentum operator L are 
defined by \/l0(dt x d)£ } , with k = -1, 0, +1. With n d 
the number operator for ri-bosons, the IBM Hamiltonian 
studied in Alhassid and Whelan (1991) is 



H = E Q + c h d + c 2 Q{x) ■ Q(x) + ciL 2 



(38) 



The Hamiltonian conserves the total number N of bosons 
and depends on the four parameters Cn,Ci,C2, and x- 
(The energy Eq is irrelevant as it only shifts the over- 
all spectrum). Except for an overall multiplicative scal- 
ing factor, only three parameters are relevant for the 
structure of the spectrum. Aside from the energy of 
the s-boson (which is put equal to zero) and the en- 
ergy of the ri-bosons (given by cq with c > because 
the energy of the ri-bosons is generically found not to 
be below that of the s-boson), H contains two scalar 
boson-boson interaction terms. The parameter x a P~ 
pearing in Q is empirically restricted by the limiting sym- 
metries described below. Moreover, there is an isomor- 
phism which maps x onto — X- Therefore, x has the range 
— v7/2 < X < 0. The parameters en and ci have opposite 
signs (the quadrupole-quadrupole is attractive). The to- 
tal angular momentum is a good quantum number; the 
value of the parameter c\ only defines the positions of 
states with different angular momenta with respect to 
one another and is irrelevant when chaos is investigated. 
Within the model, quadrupole transition matrix elements 
are calculated with the help of the same quadrupole op- 
erator Q(x) that also appears in the Hamiltonian. 

The IBM is an algebraic model and therefore easy to 
solve. Fitting data is also easy because of the small num- 
ber of parameters in Eq. ([38]) and in other versions of 
the model. The fits obtained are very good and cover a 
wide range of nuclei within a major shell. The param- 
eters change slowly with mass number. The connection 



37 



between the IBM and the shell model is understood rea- 
sonably well although it has not been possible yet to de- 
rive the IBM parameters from the shell model. All this 
suggests that the IBM encapsulates real properties of nu- 
clei and explains the success and popularity of the model. 
Studies of chaos have been restricted to the form (f38|) of 
the model which applies to even-even nuclei. It is possi- 
ble, however, to extend both the Bohr-Mottelson model 
and the interacting boson model to nuclei with an odd 
number of protons and/or neutrons. 

Tests of chaos in the IBM face the same difficulty as do 
tests of chaos in the nuclear shell model using only a sin- 
gle major shell: The IBM is well established as a tool to 
describe low-lying states. At higher excitation energies 
we expect to find states which cannot be modelled by the 
IBM. In tests of GOE predictions one needs large data 
sets and typically uses the entire spectrum of states of 
the IBM pertaining to fixed quantum numbers, thereby 
exceeding the domain of applicability of the IBM to real 
nuclei. 

The IBM has the advantage that domains (in the space 
of parameters) of full integrability coincide with symme- 
tries of the Hamiltonian. Three such symmetries exist. 
They are found using group-theoretical arguments. Let 
|0) denote the vacuum state. The six single-boson states 
s^|0) and dl\0) span a six-dimensional space. Under the 
general unitary transformations in six dimensions, i.e., 
under the group U(6), that space transforms onto itself. 
The generators of the Lie algebra of U(6) are the opera- 
tors b\bj where hi stands for one of the boson operators 
s or d^. Since H is constructed from scalar quantities 
involving these operators, H has non-vanishing matrix 
elements only between pairs of states that belong to the 
same irreducible representation of U(6). Equivalcntly, 
the eigenstates of H form bases of irreducible represen- 
tations of U(6). We deal with bosons and thus require 
totally symmetric eigenstates. These belong to the fully 
symmetric irreducible representations which are charac- 
terized by the integer TV, the total boson number. The 
three symmetries of the Hamiltonian correspond to the 
three chains of subgroups of f (6) given by 



17(6) D 17(5) D 0(5) D 0(3) ; (I) 

17(6) D SU(3) D 0(3) ; (II) 

17(6) D 0(6) D 0(5) D 0(3) . (1/7) 
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Case (I) is realized for C2 = 0, case (II) for Co = and 
X = — V^/2, and case (III) for cq = and x — 0- Each 
of the three cases corresponds to a simple dynamical sit- 
uation. In case (I) we deal with vibrational nuclei, in 
case (II) with rotational nuclei, and in case (III) with 
7-unstable nuclei (instability with respect to the equilib- 
rium value of 7). Full integrability of the Hamiltonian 
suggests classical regularity and, in the quantum regime, 
Poisson statistics for the eigenvalues. 

Alhassid et al. (1990) and Alhassid and Whelan (1991) 
checked these expectations. Using coherent states and 
the limit N — > 00, one reaches the classical limit (sim- 



ilarly as in Meredith et at, 1988.) The behavior of the 
trajectories in phase space was studied for several values 
of the total angular momentum. For fixed L 2 the param- 
eter c\ is redundant. Moreover, the Hamiltonian can be 
rescaled. This leaves two parameters, r\ (defined in terms 
of the negative ratio of Co and Nc2, with the factor N 
originating in the classical limit of the quantum Hamil- 
tonian ([38])) with < 77 < 1 and x with — \/7/2 < x < 0- 
The results are shown in Fig. 1241 The three symmetries 
in relation (f3"5)) correspond to the three corners of the tri- 
angles shown: Case (I) corresponds to 77 = 1, case (II) 
to 77 = and x = ~V7/2, and case (III) to 77 = and 
X = 0. The three corners are surrounded by regions of no 
or weak chaos as expected. Chaos is strongest for small 
angular momenta and in the center and lower part of the 
triangle. The authors emphasize the existence of a nar- 
row strip of nearly regular motion connecting the upper 
and the left-hand corners of the triangle. 

The quantum spectra were analyzed with the help 
of the NNS distribution and of the A3-statistic. The 
Porter-Thomas distribution was checked using E2 tran- 
sition probabilities. As a result, the Bohigas-Giannoni 
Schmit conjecture (see Section lil.F[) was confirmed once 
again, this time for a system of interacting bosons. More 
specifically, all three fluctuation measures are close to the 
GOE values in the domains of parameter space where the 
classical motion is close to fully chaotic, and show signifi- 
cant deviations and a tendency toward Poissonian behav- 
ior in the domains where the classical motion is almost 
regular. To obtain good statistics, the results shown were 
for N — 25 bosons and J = 2 states. This is close to the 
angular momentum value £ = 0.1 shown in the upper 
part of Fig. [2U The authors stress that for smaller, more 
realistic values of N similar effects are seen. These IBM 
studies of chaos were extended to higher energies and 
spins by including broken-pair degrees of freedom (Vrete- 
nar and Alhassid, 1992). The regular behavior predicted 
for nuclei in the "arc of regularity" which separates vibra- 
tional and rotational motion (see Fig. [24)) , was recently 
confirmed experimentally (Jolie et at, 2004). 

In summary, collective models also display quantum 
chaos. Regular features dominate near the yrast line and 
at and near symmetries of the Hamiltonian. Chaos is 
strongest in regions of parameter space that are far re- 
moved from such symmetries. This is consistent with the 
analyses of data reviewed in Sections IIII.C. 21 and IIII.CT31 



C. Special Issues 



In Sections IIV.A.21 IIV.B.2I and IIV.B.3I we have shown 
that the dynamics of nuclei often produce spectral fluc- 
tuations of the GOE type and thus quantum chaos. In 
the present Section, the emphasis is different. We ad- 
dress a number of complex situations where we cannot 
test whether RMT applies. Rather, we postulate that it 
does and use RMT as a tool to model the physical sys- 
tem. This leads to results which are used in analyzing 
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FIG. 24 Domains of regular or chaotic motion in the 77 — x~ 
plane for two values of (,, the angular momentum per bo- 
son. The accessible parameter space has the shape of a trian- 
gle. Three domains in parameter space are separated by solid 
lines: Domains with nearly regular motion (fraction of chaotic 
phase-space volume < 0.3; white areas); domains with nearly 
chaotic motion (fraction of chaotic phase-space volume > 0.7; 
dotted areas); and the regions in between. From Alhassid and 
Whelan (1991). 



data. 



1. Decay out of a Superdeformed Band 

The potential energy of deformation V(/3, 7) intro- 
duced in Section HV.B. II differs from the deformation en- 
ergy of the liquid-drop model. The difference is due 
to the shell-structure of nuclei. Nuclei with a com- 
pletely filled major shell for neutrons and/or protons (see 
Fig. [T5|) are particularly stable. But major shells well sep- 
arated from each other in energy exist not only for spher- 



ical nuclei. As the nucleus is being deformed, the energies 
of the single-particle states change, and new major shells 
emerge for certain values of the deformation parameters. 
This statement applies quite generally for single-particle 
potentials including those with a spin-orbit force. More 
generally, at fixed excitation energy the density of single- 
particle states of the shell model shows considerable fluc- 
tuations as a function of deformation. The same is true 
of the sum of the energies of all occupied single-particle 
states and, thus, of the deformation potential V((3,j). 
Beyond the smooth dependence on deformation parame- 
ters due to the liquid-drop model, V(/3, 7) therefore dis- 
plays maxima and minima that reflect the deformation 
dependence of the filled single-particle states of the shell 
model. For those values of the deformation parameters 
and for those mass numbers where a major shell is filled, 
the nuclear binding energy is expected to have a maxi- 
mum and, consequently, the potential energy of deforma- 
tion V{f3, 7) is expected to have a minimum. 

Strutinsky (1966, 1967, 1968) realized that shell clo- 
sures and, more generally, maxima and minima of V(J3, 7) 
are generic phenomena that are not limited to small de- 
formations. His "shell-correction method" is based on 
the semiclassical periodic-orbit theory for independent- 
particle motion and allows the calculation of a correction 
to the deformation energy of the liquid-drop model. The 
correction is the difference between the ground-state en- 
ergy calculated from the actual single-particle level den- 
sity (with states filled up to the Fermi energy), and that 
obtained from a smoothed level density. The method of 
calculation and many of its results are summarized by 
Brack et al. (1972) and by Brack and Bhaduri (1997). 
A particular prediction was that in certain nuclei and 
as a function of deformation, V(/3, 7) would display a 
second pronounced minimum (in addition to the abso- 
lute minimum defining the nuclear ground state). This 
is schematically shown in Figure 1251 

Following earlier evidence, the discovery of "superde- 
formed" rotational bands in 152 Dy (Twin et al., 1986) 
and other medium-weight nuclei confirmed the existence 
of the second minimum. The moment of inertia of a very 
strongly deformed nucleus is larger than usual. A rota- 
tional band corresponding to an intrinsic state located 
in the second minimum can therefore be identified by 
its particularly large moment of inertia and the resulting 
narrow spacings of its members. The intensities of the E2 
gamma transitions within a SD band show a remarkable 
feature: The intraband E2 transitions follow the band 
down with practically constant intensity. At some point, 
the intraband transition intensity starts to drop sharply. 
It either disappears abruptly or is much reduced in the 
next intraband transition(s) and then disappears. This 
phenomenon is referred to as the decay out of a SD band 
(Ragnarsson and Aberg, 1986; Herskind et al., 1987). It 
is attributed to the mixing of the SD state(s) and the nor- 
mally deformed (ND) states of the same spin and parity 
located in the first minimum, see Fig. 1251 Calculations 
using the shell-correction method show that the barrier 
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FIG. 25 The deformation energy V(/3, 7) for fixed total spin 
Jo versus deformation (schematic). In some nuclei, V(/3, 7) 
may display a second minimum at large deformation. The 
level shown in the second minimum, a member of the SD 
band with spin Jo, may decay either via gamma emission 
to the next lower level in the SD band (not shown) or, via 
tunneling, mix with the ND states with spin Jo located in the 
first minimum. The latter decay dominantly by statistical El 
emission. 



separating the first and the second minimum of the de- 
formation potential depends on and decreases with de- 
creasing spin J. Decay out of the SD band sets in at a 
spin value Jo for which penetration through the barrier 
is competitive with the intraband El decay. Theoretical 
efforts aim at a quantitaive description of that process. 
We describe here one of the theoretical approaches which 
is based on the use of random-matrix theory. 

The first minimum of the deformation potential is 
typically a few MeV deeper than the second minimum. 
Therefore, the ND states populated by decay out of the 
SD band have a typical excitation energy of 3 to 4 MeV 
above yrast. These states are believed to decay largely 
via statistical emission of El gamma quanta. Such decay 
would only contribute to the background in the coinci- 
dence spectrometer used to detect the rotational bands. 
This view is consistent with the fact that no signal for 
the decay of the ND states has been found. In view of 
the total lack of spectral information on the ND states, 
a statistical model is used in the analysis of the data. It 
is assumed that the spectral fluctuation properties of the 
ND states can be modeled by the GOE. This is the view 
taken by Vigezzi et al. (1990a, 1990b) and Shimizu et 
al. (1992), and in most subsequent work on the subject; 
see, however, D0ssing et al. (2004). It is the aim of these 
works to gain information on the barrier separating the 
first and the second minimum through the analysis of the 
decay data. 

The quantity of central interest is the probability / out 
for decay out of the SD band. A plausible formula for 



I out was obtained by Vigezzi et al (1990a, 1990b) and 
Shimizu et al. (1992) as follows. When one disregards 
the coupling of both the SD states and the ND states 
to the electromagnetic field, the resulting model for the 
Hamiltonian has the form of Eq. (|24| . with Hq^ the ma- 
trix elements for penetration of the barrier separating the 
SD state |0) and the N ND states |/x) with fx = l,...,N. 
This is in fact a doorway model for the SD state. The dif- 
ference from the usual doorway situation is that because 
of the barrier, the spreading width T^ is small compared 
to or at best of the order of the mean level spacing d of the 
ND states. Let \m) with m = 1, 2, . . . , N + 1 denote the 
eigenstates of H and c m = (0|m) the amplitudes with 
which the SD state |0) is admixed into the eigenstates 
\m). It is assumed that E2 decay out of the next-higher 
state in the SD band populates the state \m) with proba- 
bility |c m | 2 , and that the widths for decay of the state \m) 
back into the SD band or by statistical emission of El 



gamma rays are given by \c m \ Eg and by (1 



1 )Vn, 



respectively. Here T$ is the width for electromagnetic 
decay within the SD band, and T^ is the common total 
decay width for El emission of the ND states (it has the 
same value for every state |/x)). We then have (Vigezzi 
et at., 1990a, 1990b; Shimizu et al, 1992) 



Z^ 



\Cm I 



(i-\ Cm \ 2 )r N 



(i-| Cm | 2 )iv + |c m | 2 r s 



(40) 



To compare with the measured decay probability out of 
a SD band, the ensemble average of this expression was 
numerically simulated by putting Eq in the middle of 
the GOE spectrum, by diagonalizing a large number of 
random matrices of the form (|24[) and by using the eigen- 
functions to calculate the coefficients c m . Repeating the 
procedure for different strengths of the barrier penetra- 
tion matrix elements, one determines the actual strength 
by a fit to the data. 

The arguments used to write down Eq. (|40p imply a 
perturbative treatment of the coupling to the electromag- 
netic field. As a result, J ut depends on two dimensionlcss 
parameters, the ratio Tn/Ts and the ratio T^/d which 
determines the coefficients c m . Thus, I ou t is indepen- 
dent of the fine-structure constant. This is physically 
implausible because once the SD state is populated, the 
competition is between the intraband decay and the pop- 
ulation of the ND states and thus between the electro- 
magnetic and the strong interaction. More importantly, 
the analyses of data using the approach of Vigezzi et al 
(1990a, 1990b) and Shimizu et al. (1992) yielded values 
of r^ that were about two orders of magnitude smaller 
than r^r, putting a perturbative treatment of Tjv into 
question. In Weidenmiiller et al. (1998) and Gu and 
Weidenmiiller (1999) the same statistical model (|2"4")l was 
used in a non-perturbative way. The amplitudes for in- 
traband decay of the SD state and for El decay of the ND 
states were added as diagonal terms on the right-hand 
side of Eq. (|M)l . An expression for the amplitude for de- 
cay out of the SD band was derived using that modified 
Hamiltonian. The ensemble average of that amplitude 



was calculated analytically. The resulting expression for 
the average probability for decay out of the Sd band in- 
volves a threefold integral and depends on the dimen- 
sionless parameters T^/d and Ts/r^. The result was 
compared with that of Vigezzi et al. (1990a, 1990b)and 
Shimizu et al. (1992) and the limits of validity of the 
latter approach were determined. Analytical approxima- 
tion formulas to the exact result were given to simplify 
the data analysis. Many data have been analyzed using 
either approach, see, for instance, Kruecken et al. (2001). 
and references therein. A simplified treatment was devel- 
oped by Stafford and Barrett (1999). 

A few figures taken from Kruecken (2001) may serve as 
examples for the results obtained. Superdeformed bands 
occur not only for mass numbers around 150 (where they 
were first discovered) but also, for instance, in the lead 
region (A rs 200). A case in point is the first SD band 
in 194 Hg. For spin J = 12h, one finds T s = 0.097 meV, 
Tat = 4.8 meV, and d — 16.3 eV. The analysis using 
Vigezzi (1990a, 1990b) yields T l = 37 meV and that 
using Gu and Weidenmuller (1999) yields 1^ — 25 meV. 
Calculating the rms barrier penetration matrix element 
v from r^ and using results for several values of J, one 
obtains a dependence on J roughly in agreement with 
the theoretical expectation v ex exp(— aJ). 

Decay out of a SD band continues to receive consider- 
able attention, see, for instance, Sargeant et al. (2005) 
and references therein. Our aim here was to show how 
RMT is basically used in analyzing data. 



2. Double Giant Dipole Resonance 
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Here we focus on an explanation originally due to Carl- 
son (1999a, 1999b) which uses the Brink-Axel hypothe- 
sis. The Brink-Axel hypothesis postulates the existence 
of a giant dipole resonance built not only on the nuclear 
ground state, but also on every excited nuclear state as 
well (Brink, 1955; Axel, 1962). With the above picture 
for the GDR as an oscillation of the center-of-mass of 
the protons against that of the neutrons, the hypothesis 
is very plausible. The hypothesis suggests an enhance- 
ment of the excitation cross section of the DGDR by the 
following mechanism. The GDR is a doorway state and 
mixes with the background states having the same spins 
and parities as the GDR. The mixing time is simply esti- 
mated as h/rK If that mixing time is small or at most of 
the order of the time it takes to excite the DGDR from 
the GDR, then each of the background states admixed 
to the GDR mode may get excited into its own GDR. 
The contributions from the excitation of the background 
states would add to that of the GDR and lead to an en- 
hancement of the cross section. With a GOE model for 
the background states as in Eq. (|24p . the average inten- 
sity for excitation of the DGDR in a collision between 
two heavy ions can be worked out using controlled ap- 
proximations (Gu and Weidenmuller, 2001). For realistic 
values of the parameters, it is found that the contribu- 
tion of the background states is significant, and that the 
results of the calculations agree well with data on the re- 
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^Pb at 640 MeV/nucleon. The DGDR 



is, thus, another case where RMT is successfully used to 
understand data. The use of RMT vindicates the Brink- 
Axel hypothesis. 



We recall the discussion of the Giant Dipole Resonance 
(GDR) at the beginning of Section III.GI Action of the 
dipole operator on the nuclear ground state induces the 
dipole mode. That mode can be viewed as an oscilla- 
tion of the center-of-mass of the protons against that 
of the neutrons. This simple intuitive picture is exact 
in the framework of two extremely opposite models of 
nuclear structure, the harmonic-oscillator independent - 
particle model, and a collective model using neutron and 
proton fluids. Therefore, the picture is expected to have 
general validity. If the oscillation is approximately har- 
monic, a repeated El excitation of the GDR should be 
possible and should lead to the double giant dipole res- 
onance (DGDR). The DGDR does indeed exist and was 
first observed in a number of nuclei in the 1990s. For a re- 
view, see Bertulani and Ponomarev (1999). Compared to 
predictions of the harmonic picture, the measured cross 
sections for excitation of the DGDR were found to be 
larger by factors ranging from 1.3 to 2. The widths of 
the DGDR were found to be about 1.4 times larger than 
those of the GDR. This last fact is in keeping with some 
but not all theoretical estimates. These findings have at- 
tracted much theoretical attention. Anharmonicities of 
the Hamiltonian and non linearities of the external field 
were studied as possible causes for the discrepancies. 



3. Damping of Collective Modes and Friction 

Collective motion is due to the coherent motion of 
many nucleons. Examples are the rotational or vibra- 
tional motion of nuclei (Section HV.B.1[) or the giant 
dipole resonance, an oscillation of the center of mass of 
the neutrons against that of the protons (Sections III.GI 
and HV.G2J) . Other examples are encountered in heavy- 
ion collisions and in induced nuclear fission. In the first 
case, the grazing collision of two heavy nuclei at ener- 
gies of a few MeV per nucleon (at the Coulomb bar- 
rier) leads to a process known as deep inelastic scatter- 
ing (Norenberg and Weidenmuller, 1980; Weidenmuller, 
1980): The kinetic energy of relative motion and the 
associated angular momentum are partly transformed 
into intrinsic excitation energy and spin of both frag- 
ments. At the same time, nucleons are transferred be- 
tween the two reaction partners. As a result, the reac- 
tion mainly produces pairs of fragments in highly excited 
states with masses similar to those of the incident nuclei, 
but with considerably smaller kinetic energies. Similar 
processes occur at higher incident energies (Cassing and 
Norenberg, 1985). The excitation energies of either frag- 
ment cannot be measured precisely. A phenomenological 
description of the process focuses attention on the collec- 
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tive degrees of freedom, i.e., on the relative coordinates 
of both fragments. These move under the influence of a 
conservative potential (due to the Coulomb force and to 
the overlap of both nuclei) and of dissipative forces. The 
latter account for the transfer of energy and angular mo- 
mentum into intrinsic degrees of freedom as well as for 
nucleon transfer between both reaction partners. In the 
simplest approximation, the relative motion is described 
classically and the dissipative forces are represented by 
a friction constant (Wilczynski, 1973). Induced fission is 
a somewhat similar process: Here the shape deformation 
leading to fission is identified as the collective degree of 
freedom. For fission to happen, that degree of freedom 
must overcome the fission barrier. Between the top of the 
fission barrier and the scission point, the fission degree of 
freedom gains kinetic energy. That energy is partly trans- 
formed into intrinsic excitation energy of the fissioning 
system. The process can again be described phenomeno- 
logically as involving dissipative forces acting on the fis- 
sion degree of freedom (Grange and Weidenmuller, 1980). 
The number of neutrons emitted prior to fission serves as 
a measure of the strength of that dissipation (Gavron et 
at, 1986). In nuclear physics, dissipative processes typ- 
ically arise in the context of nuclear reactions. On the 
other hand, the theoretical treatment usually does not 
refer to scattering processes. That is why we deal we 
dissipation here and not in Part 2 of this review. 

The description of dissipative forces due to the interac- 
tion of a quantum system with its environment is a stan- 
dard topic in nonequilibrium statistical mechanics. If the 
/ collective degree(s) of freedom q v with v — 1,2,...,/ 
can be treated classically, the phenomenological descrip- 
tion of the collective motion may use an equation of the 
type 



dt dq v dq v 



(41) 



Here C(q v ,q v ) is the classical Lagrangian, t is the time, 
and F v are the friction forces. To get the most general 
description, stochastic Langevin forces must be added to 
the friction forces. An equivalent description is in terms 
of a Fokker-Planck equation. If the classical approxi- 
mation is not justified, Eq. (|4"Tj) must be replaced by its 
quantum analog which typically involves the density ma- 
trix. Often it is justified to characterize the friction forces 
by a friction constant and the Langevin forces by a diffu- 
sion constant. We refer to these constants as to transport 
coefficients. Calculating transport coefficients is the aim 
of a microscopic approach to these processes. 

The canonical approach is that of Caldeira and Legett 
(1983), who described the environment (the "heat bath") 
as an infinite set of harmonic oscillators. The nuclear 
case differs from the standard situation in statistical me- 
chanics in several respects. First, the intrinsic degrees of 
freedom cannot be viewed as an infinitely extended heat 
bath. The total energy of the system (intrinsic plus col- 
lective degrees of freedom) is conserved, and the energy 
content of either subsystem is roughly the same. In other 



words, the nucleus is a mesoscopic system for which the 
thermodynamic limit is not appropriate (although some 
statistical description may still apply). Second, the cou- 
pling between the collective and the intrinsic degrees of 
freedom is typically very strong. The energy content of 
the interaction part of the Hamiltonian is not small com- 
pared to that of either subsystem. As a consequence, the 
usual Markov approximation often does not apply (Brink 
et al., 1978). These facts must be taken into account 
as one works out the transport coefficients appearing in 
Eq. (|41[) or in the quantum analog of that equation from 
a microscopic approach. 

The microscopic approach itself also encounters some 
difficulties. First, the definition of the collective degree(s) 
of freedom in terms of the coordinates of the partici- 
pating nucleons poses problems. Second, in nuclei dissi- 
pative processes typically involve the transfer of several 
MeV or even several ten MeV from the collective to the 
intrinsic degrees of freedom, perhaps together with the 
transfer of several ten units of angular momentum. The 
level densities in the fragments at the resulting excitation 
energies are very high. A detailed microscopic descrip- 
tion of the intrinsic degrees of freedom is therefore not 
possible. Rather, one employs stochastic models of the 
random-matrix type. In such models, the coupling ma- 
trix elements between collective and intrinsic degrees of 
freedom are Gaussian-distributed random variables (in 
accord with the Porter-Thomas distribution in Eq. [|18p). 
The mean coupling strength and the density of intrinsic 
states are estimated using simple nuclear models. As a 
result, the force acting on the collective degree(s) of free- 
dom is random, and the transport coefficients are defined 
as ensemble averages involving that random force. 

Calculating the transport coefficients as ensemble av- 
erages poses physical questions and technical difficulties 
quite different from those encountered in calculating the 
average of the strength function J2 T l(0| T )| 2 ^(-^ — E T ) in 
Section III.GI That calculation leads to Eq. ([23]) . In a 
time-dependent picture, Eq. (|25|) shows that the ampli- 
tude of the collective mode decays exponentially due to 
mixing with the background states. The decay width is 
given by the spreading width of Eq. (f2"6"|) . (To see that it is 
the amplitude of the collective mode which undergoes ex- 
ponential decay we note that the strength function can be 
expressed identically in terms of the Green's function as 
-(l/7r)3(0|(£+ -iTr 1 ^)). In the case of transport co- 
efficients, it is necessary to average probabilities (squares 
of amplitudes). Quantum mechanics implies that total 
sum of all occupation probabilities is conserved in time. 
Therefore, transport coefficients describe the redistribu- 
tion in time of occupation probabilities caused by dissipa- 
tive forces. Friction, for instance, implies that collective 
states at lower energy become occupied at the expense of 
states at higher energy. Averaging expressions that de- 
pend on squares of amplitudes is technically harder than 
averaging amplitudes. 

The description of dissipative processes in heavy-ion 
reactions and of fission in terms of RMT has received 
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much theoretical attention. Aside from the reviews 
(Norenberg and Weidenmiiller, 1980; Weidenmuller, 
1980) and in addition to the papers cited above, we men- 
tion the work of Agassi et al. (1977) and of Hofmann 
and collaborators (1977, 2001). In a very general con- 
text, the interaction of collective and microscopic degrees 
of freedom has been studied in a series of papers by Bul- 
gac, Kusnezov, and collaborators (see Bulgac and Kus- 
nezov (1996) and references therein). There are other 
approaches to nuclear dissipation which do not use RMT 
but emphasize the chaotic aspects of single-particle mo- 
tion, see Blocki et al. (1995) and references therein. 



4. Fluctuations of Binding Energies 

The masses of atomic nuclei are keys to the under- 
standing of many physical and astrophysical processes. 
For this reason it is important to construct reliable the- 
oretical models for the values of the nuclear masses or, 
equivalently, for the binding energy B{A) as function of 
mass number A. (We suppress for simplicity the addi- 
tional dependence of B on neutron number N or proton 
number Z). This function is also needed to predict the 
masses of (as yet) unknown nuclei. 

The standard approach to a global modeling of the 
function B(A) starts out from the liquid-drop model of 
the nucleus and considers in addition shell corrections 
(see Section HV.C.ip . as well as corrections due to the 
pairing force (see Section HV.A.2[) . The latter lead to 
an odd-even staggering of B(A). The resulting "semi- 
empirical mass formula" contains about 30 parameters 
and is fitted to a large number of data. Years of painstak- 
ing work have culminated in a best fit (Moller et al., 1995) 
that reproduces the data points very well but not exactly. 
The overall difference (root-mean-square value taken in 
a limited window of mass values) is of the order of 0.5 
MeV and decreases with increasing A. Other approaches 
(Samyn et al., 2004; Duflo and Zuker, 1995) have led to 
similar differences. The figure of 0.5 MeV is obviously 
very small and of the order of 5 x 10~ 4 in comparison 
with the total binding energy which for medium-weight 
and heavy nuclei is of the order of GeV. Nevertheless, 
that small value has attracted considerable attention. 

Bohigas and Leboeuf (2002) have suggested that there 
are two types of contributions to the shell correction for 
the nuclear binding energy. The first one is due to the 
regular motion of nucleons in the mean field and is taken 
into account in terms of the Strutinsky shell correction 
method or, in the fits just mentioned, in terms of the 
semiempirical mass formula. The second one is due to the 
(partly) chaotic motion of nucleons within the nucleus. 
In Bohigas and LeBoeuf (2002) that part is also evalu- 
ated within the mean-field approximation. The authors 
argue, however, that their final result for the fluctuations 
of the chaotic part is of much more general validity, and 
may be interpreted as arising from the residual interac- 
tions. 



With p(E) the density of single-particle states, the 
shell correction to the binding energy has the form 
J dE Ep(E). Fluctuations of the binding energy are 
characterized by the variance of that expression, the av- 
erage being taken over a window of mass numbers and 
indicated by angular brackets. That variance can be ex- 
pressed identically in terms of the form factor K (t) as 
{B\A)) - {{B{A))f = (n 2 /(2^))f ° c drK(r)/r i . The 
form factor is essentially the Fourier transform of the 
density-density autocorrelation function. For the chaotic 
contribution to the variance, K(t) is approximated by 
the random-matrix result K(t) = It. This expression 
applies for values of r below the Heisenberg time, and 
for t > T m i n where r m i n is the period of the shortest pe- 
riodic orbit in the system while K(t) = for r < r m i n . 
The value of T m - nl is estimated using free motion at the 
Fermi velocity over a distance typically given by the nu- 
clear radius. All this yields for the root-mean-square 
value of the fluctuation <j(A) of the binding energy the 
result a{A) = 2.78A" 1 / 3 MeV. That result is in reason- 
able agreement with the deviations of the data from the 
semiempirical mass formula. 

The interpretation proposed by Bohigas and Leboeuf 
(2002) of the deviations of measured binding energies 
from the semiempirical mass formula as being due to 
chaotic motion has caused a lively discussion that cannot 
be reviewed here. It implies that throughout the valley 
of stability, the binding energies of all nuclei are corre- 
lated (Molinari and Weidenmuller, (2006); Olofsson et al. 
(2006)). Such correlations, while expected on dynamical 
grounds, are surprising from a statistical viewpoint. 



V. RANDOM-MATRIX MODELS INSPIRED BY 
NUCLEAR-STRUCTURE CONCEPTS 

In Section IIVI is was shown that in the framework of 
standard nuclear models such as the shell model or the 
cranking model, chaos is a generic property of nuclei. 
However, these models differ very much from the GOE, 
the standard ensemble to model chaotic systems. To mo- 
tivate the study of the random-matrix ensembles treated 
in the present Section, we elucidate the difference be- 
tween these nuclear models and the GOE, more precisely: 
Between these models and a generic realization of the 
GOE. In doing so, we focus attention on the shell-model 
Hamiltonian (|34[) ; what will be said applies mutatis mu- 
tandis likewise to the Hamiltonian of the cranking model 
in Section HV.B. 21 For the sake of simplicity we disregard 
total spin and isospin. A more complete discussion for 
the sd-shell including these quantum numbers may be 
found in Section IV. Bl 

Without spin and isospin, the many-body eigenstates 
of the single-particle part of the shell-model Hamilto- 
nian (|34p are Slater determinants. These are multiply 
degenerate. The degeneracies are lifted by the residual 
interaction V rcs - The way this happens depends on the 
rank k of V vcs . If V rcs is a pure two-body interaction 
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(k = 2), then V ICS has non-vanishing matrix elements 
only between those pairs of Slater determinants which 
differ in the occupation numbers of at most two single- 
particle states. In the Hilbert space spanned by Slater 
determinants, the matrix of V les then has many zeros. If 
Vres also contains three-body forces (k = 3), the number 
of zeros is reduced because V^ es now connects all Slater 
determinants which differ in the occupation numbers of 
at most three single-particle states. The systematic ap- 
pearance of zeros in the matrix representation of V les is 
completely removed only when V Ies contains m-body in- 
teractions where m is the number of valence nucleons. 
For a generic realization of the GOE, on the other hand, 
zeros do not appear systematically. This is why it is 
sometimes said that the GOE contains interactions of 
arbitrary rank. 

When the first evidence for the validity of a RMT ap- 
proach to nuclei became available in the 1960s, the ques- 
tion arose how that difference between the GOE and a 
realistic residual interaction would affect nuclear proper- 
ties. Would the spectral statistics be the same? Later, 
other nuclear properties (like the average level density) 
came into the focus of RMT. How would these depend 
on the rank of V^-cs? To answer these questions, sev- 
eral non-canonical random-matrix ensembles were intro- 
duced. French and Wong (1970) and Bohigas and Flores 
(1971) defined and studied what became known as the 
two-body random ensemble (see Section IV. Bp . A few 
years later, Mon and French (1975) defined the embedded 
/c-body random ensembles where the rank of the residual 
interaction is a free parameter (see Section fV-Aj) . Quite 
recently, constrained ensembles (where certain matrix el- 
ements of the GOE are suppressed) were introduced in 
Papenbrock et al. (2006) (see Section |V.C| . The present 
Section is devoted to a review of these ensembles. The 
analytical treatment of the first two ensembles is much 
more difficult than that of the canonical ensembles of 
RMT because they lack the orthogonal invariance of the 
GOE. Therefore, a complete analytical theory of these 
ensembles does not exist. The analytical treatment of 
the constrained ensembles is yet in its infancy. 



A. Embedded Ensembles 

For the fc-body embedded ensembles of random ma- 
trices (introduced by Mon and French, 1975) we confine 
ourselves to a brief account of the main features of the 
orthogonal case and refer to the reviews (Brody et al., 
1981; Kota, 2001; Benet and Weidenmiiller, 2003) for 
further details. 

The embedded ensembles dispose of all the complexi- 
ties due to the couplings of angular momentum and spin 
but retain the symmetries imposed by the exclusion prin- 
ciple. One considers m fermions in I > m degenerate 
single-particle states (which carry no further quantum 
numbers) that interact via a random &;-body interaction 
with k < m. To obtain an understanding of the transition 



from the two-body ensemble to the GOE, the parameter 
k is allowed to range from k = 1 to k = m, although 
k = 2 is the most interesting (i.e., realistic) case. The 
case k = 1 is integrable and, therefore, somewhat excep- 
tional. 

Labeling the single-particle states with a running in- 
dex j = 1, . . . , I, we introduce the usual fermionic cre- 
ation operators a, and define the creation operators for 
k fermions by 






IT 



(42) 



where a stands for the set ji < J2 < ■ ■ ■ < jk ■ The cor- 
responding annihilation operators are defined by ip k = 

(ip k y . The random A:-body interaction is 
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The coefficients Vk-a-y are Gaussian-distributed random 
variables with zero mean values and a common second 
moment. The matrix elements Vk-a-y are real symmet- 
ric with a second moment Vk-a^Vk-a'-y' = v 2 {S a a'S Ty > + 
SajiSya'). The Kronecker deltas stand for the string 
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. Sj k ji etc. The interaction Vk lifts the de- 



generacy of the many-body states; the parameter v 2 sets 
the scale for the spectrum. Without loss of generality we 
may put v 2 — 1. 

The fc-body embedded ensembles are then defined in 
terms of the three parameters k,l,m with k < m < I. 
The ensembles are referred to as EGOE(fc) and jointly as 
EGE(fc). In canonical RMT, universal results are ob- 
tained in the limit of infinite matrix dimension. For 
EGE(fc), the same limit is obtained by taking I — > oo. 
For fixed k this can be done by imposing constraints on 
the ratio m/l. Brody et al. (1981) define the dilute limit 
as I — > oo, m — ► oo, m/l — ► 0. 

Central questions in the theory of the embedded en- 
sembles are: (i) What is the shape of the spectral den- 
sity? (ii) What are the spectral fluctuation properties? 
(hi) Are these properties universal (i.e., independent of 
the assumed Gaussian distribution)? (iv) Are the spec- 
tra ergodic? (v) Can the embedded ensembles be usefully 
applied to real nuclei? If so, what are their predictions? 
Ideally, questions (i) to (iv) should be answered in the 
limit I — > oo. 

Partial answers to these questions have been obtained 
with a variety of methods. The moments method (Mon 
and French, 1975) evaluates moments of Vk in the limit 
I — > oo. The distribution of Hamiltonian matrix elements 
being Gaussian, products of Hamiltonian matrices are av- 
eraged by Wick contraction involving all pairs of matri- 
ces. In the dilute limit k <C rn <C I, m/l — ► for I — » oo, 
only pairs of neighboring matrices are taken into account 
("binary correlation approximation"). Numerical meth- 
ods give some insight although the extrapolation to large 
matrix dimension may pose problems. The study of the 
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second moments of the many-body matrix elements of 
Vk reveals a duality symmetry between EGOE(fc) and 
EGOE(m-fc). 

For the spectral density, there is a gradual transition 
from the semicircular shape (which is attained for fc = m) 
to the Gaussian shape (which applies for k <C m). The 
transition sets in at 2fc = m. Less is known for the level 
statistics. It is clear that for k = m the spectral fluctu- 
ations are Wigner-Dyson-like and that for fc = 1 they 
are Poissonian. The cases 1 < fc < m have been much 
debated without firm analytical conclusion. The numer- 
ical evidence points toward Wigner-Dyson statistics for 
fc > 2. Universality has not been addressed. Ergodicity 
has been proved for some observables in the limit I — > oo, 
but the non-ergodic contributions disappear very slowly 
with increasing I. 

We turn to the applications of EGOE(2) to nuclei. 
These are based upon the binary correlation approxi- 
mation (Kota, 2001). Moreover, it is stipulated that 
the suppression of spin and isospin quantum numbers 
does not limit the predictive power of EGOE(fc) in nu- 
clei. With that assumption, the EGOE(2) result for 
the spectral density explains why shell-model calcula- 
tions with a large number of valence nucleons yield ap- 
proximately Gaussian spectra. But a Gaussian spectrum 
results also for fc = 1 irrespective of m, so that spec- 
tra of approximately Gaussian shape are generically ex- 
pected. EGOE(fc) also predicts the distribution of tran- 
sition strengths. For an operator O causing a transition 
from energy £!j to energy Ef, the transition strength dis- 
tribution is defined as the ensemble average of the trace 
of OU(V 2 - E f )OS(V 2 - Ei). This expression equals the 
square of the transition matrix element multiplied with 
the densities of the initial and final states. The average 
can be worked out and yields a bivariate Gaussian dis- 
tribution. The comparison with shell-model calculations 
in Kota (2001) shows good agreement. Similarly, tran- 
sition strength sums can be worked out in closed form. 
The same is true for the average occupation numbers of 
single-particle states which were referred to already in 
Section IIV.A.2I and which also agree well with results of 
the shell model. The method can be extended to cases 
where the Hamiltonian is the sum of a single-particle op- 
erator and V 2 . It is then possible, for instance, to predict 
the onset of chaos versus the strength of V 2 , again in good 
agreement with numerical calculations. 

In summary we see that EGOE(2) yields a number of 
results that are in good agreement with the shell model. 
That means, in turn, that EGOE(2) is capable of making 
predictions that can reliably be used when shell-model 
calculations are not available. All this is in stark contrast 
to the GOE and is possible only because EGOE(2) takes 
account of an essential aspect of the shell model. On 
the other hand, it is very difficult to obtain analytical 
results for EGOE(fc) which are not based upon the binary 
correlation approximation. Perhaps most importantly, 
there still is no definitive analytical result on the spectral 
fluctuation properties of EGOE(fc). 



The embedded ensembles have recently been general- 
ized to cover particles with spin (see Kota, 2007). 



B. Two-body Random Ensemble 

The two-body random ensemble (TBRE) addresses 
the questions: Which nuclear properties obtained by di- 
agonalizing the shell-model Hamiltonian (f3"5)) are generic 
(i.e., hold for most two-body interactions), and which 
are specific properties of a given interaction? And how 
does the residual interaction mix the states so as to pro- 
duce chaos in nuclear spectra? To this end, the TBRE 
uses the actual form (|35[) of the shell-model Hamilto- 
nian but replaces the matrix elements v a of the actual 
two-body interaction by Gaussian-distributed real ran- 
dom variables with zero mean value and a common sec- 
ond moment (multiplied by the factor 2 for the diago- 
nal elements, see Eq. (|4"5j0 . In contrast to the embed- 
ded two-body random ensemble (see Section IV. A[) , the 
TBRE obviously does take into account spin and isospin 
quantum numbers. It is time-reversal invariant. More- 
over, the TBRE explores the properties of the residual 
two-body interaction uniformly in the space spanned by 
the variables v a . Statements derived for the TBRE apply 
to almost all two-body interactions with the exception of 
a set of measure zero. 

The TBRE was introduced by French and Wong (1970) 
and Bohigas and Flores (1971). The authors were mainly 
interested in the spectral fluctuation properties of the 
TBRE. The numerical results reported by French and 
Wong (1970) and Bohigas and Flores (1971) showed that 
the NNS distribution and the A 3 -statistic of the TBRE 
agree with those of the GOE. These results were con- 
firmed by later numerical studies and showed that chaos 
is a generic property of the TBRE. Interest in the TBRE 
was revived in 1998 by the work of Johnson et al. (1998, 
1999). These authors showed that in even-even nu- 
clei, the TBRE predicted ground states with spin zero 
and positive parity much more frequently than corre- 
sponds to their statistical weight, in spite of the fact 
that the matrix elements v a are random. That finding 
caused substantial theoretical activity (see the reviews by 
Zelevinsky and Volya, 2004 and Zhao et al, 2004) includ- 
ing studies of bosonic systems interacting via a random 
two-body interaction. Here we confine ourselves to the 
nuclear-physics aspects of the TBRE. Papenbrock and 
Weidenmuller (2004, 2005, 2006, 2007) investigated the 
mechanism by which the TBRE mixes the shell-model 
configurations. 

The TBRE contains the non-degenerate single- 
particle energies (first term on the right-hand side of 
Eq. (|35[) ) as non-random parameters. The mixing of 
shell-model configurations in the TBRE depends on the 
mean strength of the two-body matrix elements mea- 
sured in units of the spacing of the single-particle levels. 
Complete mixing occurs only when that ratio is large, 
see Section IIV.A.21 For simplicity and clarity, theoreti- 
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cal studies of the TBRE focusing on chaos often assume 
that the single-particle energies are degenerate (in which 
case they can be put equal to zero by a shift of the energy 
scale) . We follow that custom here and study the mixing 
of shell-model configurations in a "pure" TBRE. It has 
to be borne in mind, however, that the non-degeneracy 
of the single-particle levels tend to weaken the mixing 
found in that model. 



1. Comparison GOE - TBRE 

Before discussing specific properties of the TBRE, it is 
instructive to compare the TBRE with the GOE. In the 
GOE and with N » 1 the matrix dimension, the num- 
ber of independent random variables is N(N + l)/2 and, 
thus, large compared to N. In the TBRE, on the other 
hand, the number J\f of independent random variables is 
typically small compared to the dimension T>( J, T, II) of 
the Hamiltonian matrix. We recall that in the sd-shell, 
we have Af = 63 while typically T> « 10 3 in the middle 
of the shell and for low values of J. In the p /-shell, the 
corresponding figures are Af = 195 and T> w 10 4 to 10 5 . 
The complete mixing of the basis states and the ensu- 
ing validity of Wigner-Dyson statistics for the spectrum 
cannot be achieved by such a small number of random 
variables alone. In an essential way it is also due to the 
matrices C^ v (JTH;a) appearing in Eq. (|35[) . Therefore, 
studies of the TBRE must focus on the structure of these 
matrices, see Section IV. B. 21 

The GOE is mathematically accessible and formally 
attractive because it is orthogonally invariant, univer- 
sal, and ergodic (see Section III.CI) . By construction, the 
TBRE is much more realistic than the GOE (if one in- 
cludes the non-degenerate single-particle states of the 
shell model in the Hamiltonian) but probably lacks all 
these properties. It is not orthogonally invariant because 
the matrices C^(JTH;a) are fixed by the shell model. 
An orthogonal transformation would change the chosen 
representation of these matrices in Hilbert space but 
would leave every realization of the TBRE unchanged. 
It is not clear whether the TBRE is universal (i.e., yields 
results which do not depend on the assumed Gaussian 
distribution of the matrix elements v a ). We are not aware 
of any paper addressing that question. The TBRE is not 
ergodic because the limit of infinite matrix dimension 
cannot be taken in a meaningful way (except for the case 
of a single j— shell where j — > oo is a meaningful limit 
that has not been explored yet) . In spite of these short- 
comings, the TBRE has attractive features, see below. In 
Section III.CI it was pointed out that GOE spectra carry 
no information content. The TBRE produces spectra 
with Wigner-Dyson level statistics. At the same time, 
the TBRE does carry information content because the 
number of random variables is small compared to typi- 
cal matrix dimensions. Ideally it takes Af data points to 
completely determine the values of the random variables 
in the TBRE; that number is typically small compared 



to the number of eigenvalues pertaining to fixed values 
of J, T, and n. Again, this shows the important role 
played by the matrices C^(JTU; a) in the TBRE. These 
matrices are fixed by the geometry of the shell model it- 
self. At the same time, they are obviously very important 
for the strong mixing of the shell-model configurations. 
The choice of the residual interaction only determines the 
particular linear combination of the C's that forms the 
shell-model Hamiltonian H^iJTU) in Eq. ([35). 

Properly speaking, the TBRE is not a single ensemble 
but a set of ensembles. Indeed, in every shell a given set of 
matrix elements {v a } determines for all values of A per- 
taining to that shell the Hamiltonian matrices H^lJTH) 
for all values of J, T, and n. Taking the v a as random 
variables implies that all these matrices become Gaussian 
random-matrix ensembles. Since all these ensembles de- 
pend upon the same set {v a } of random variables, they 
are correlated. Such correlations are the hallmark of the 
TBRE. Correlations of this type do not occur naturally 
within the GOE. 



2. Structure of the Matrices C M „(JTn;a) 



The discussion in Section IV.B.ll has revealed the cen- 
tral role played by the matrices C^„ (jni; a) in the 
TBRE. As emphasized in Section HV. A. li these matrices 
are completely determined by the coupling scheme cho- 
sen to construct the many-body states \JTH/j,). They 
depend upon vector-coupling coefficients and on coef- 
ficients of fractional parentage and, thus, are given in 
terms of group-theoretical concepts. (The coefficients of 
fractional parentage give the decomposition of a state 
| JTUfi) for to nucleons in terms of the same states for 
(to — 1) nucleons). The matrices C M „ (J r TH; a) are deter- 
mined by the intrinsic symmetries of the shell model and 
are known completely. At the same time, these matrices 
are highly complex, and there is no analytical theory yet 
to describe their structure. We are confined to describ- 
ing briefly some of their properties and refer for details 
to Papenbrock and Weidenmiiller (2007). 

To see how the non-diagonal elements of V ras mix the 
unperturbed configurations, we use the coupling scheme 
described above Eq. (|3"5|) and classify the many-body 
states in terms of the partitions {mij} of to. For to = 12 
there are 41 such partitions in the sd-shell. Each par- 
tition consists of a string of three non-negative integers 
which give the number of nucleons in each subshell. In 
that basis, the matrix elements of V ras attain block struc- 
ture. Each block is labeled by a pair of partitions {mtj}. 
The blocks come in four classes: (i) The partition {mgj} 
is unchanged (diagonal blocks arising from those two- 
body matrix elements which do not change partition al- 
though they may change the actual many-body basis 
state); (ii) the partition {m^-} is changed by adding unity 
to one of its elements at the expense of another element 
(the residual interaction lifts a single nucleon from one 
subshell to another); (iii) the partition {mtj} is changed 
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by adding unity to two of its elements or two to one of 
its elements at the expense of one or two other elements 
(the residual interaction lifts two nucleons into different 
subshells); (iv) the partition {rngj} is changed by moving 
more than two nucleons into different subshells. These 
blocks are empty because a two-body interaction cannot 
move more than two nucleons. 

The block structure of the matrices C^ ( JTH; a) (also 
displayed in Zelevinsky et at, 1996) is here shown in 
Fig. [26] for m — 12 and for the states with J = 0, 
T = of the s<i-shell. Blocks in classes (i), (ii), (iii), 
and (iv) are shown in black, light gray, dark gray, and 
white, respectively (color online: red, green, blue, and 
white, respectively), ft is obvious that the structure of 
the matrix is fundamentally different from that of a typ- 
ical GOE matrix where all states are coupled with each 
other. That fact implies that complete mixing of the 
shell-model configurations can never be the result of a 
single matrix C Mt ,(JTTl;a) alone. A linear combination 
of the type appearing in Eq. ([35]) is definitely needed, ft 
is also obvious that the block structure shown in Fig. [2S] 
is not restricted to the sd-shell and is a generic property 
of the shell model. 
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FIG. 26 (Color online) Block structure of the matrix of the 
shell-model Hamiltonian for the sd-shell. Blocks in black, 
light gray, dark gray (color online: red, green, blue) indicate 
matrix elements which change the partition by zero, one or 
two units, respectively, as explained in the text. White areas 
do not carry matrix elements. From Papenbrock and Wei- 
denmiiller (2005). 

The structure of individual matrices C^(JTH;a) is 
governed by "geometric chaos" (Zelevinsky et al., 1996). 



The many-body states \JTHfj,) may be constructed by 
coupling first a pair of nucleons to given intermediate val- 
ues of spin and isospin, by then vector-coupling a third 
nucleon to the resulting pair etc. etc. There are obviously 
very many different paths leading to the same total spin 
J and isospin T; their number determines the dimension 
T>( J, T, II) of the resulting Hilbert space. Each path cor- 
responds to a different set of vector-coupling coefficients. 
Supposing that the vector-coupling coefficients are pseu- 
dorandom numbers, we expect that the states \JTHfi) 
also behave randomly. Several tests (Zelevinsky et at, 
1996; Zelevinsky and Volya, 2004) support this idea of 
"geometric chaos" . It is, therefore, reasonable to expect 
that the elements of the matrices C^ v {JTIi\a) behave 
randomly in those blocks where they do not vanish identi- 
cally. This view is obviously not restricted to the sd-shell 
but applies likewise to every major shell. Further insight 
into the mixing properties of the matrices C '^(JTH; a) 
is obtained by counting the number of matrices which 
yield non-vanishing contributions to a matrix element of 
the shell-model Hamiltonian 7i M „(JTH) and by calcu- 
lating the inverse participation ratios for these matrices 
(Papenbrock and Weidenmuller, 2005). All this evidence 
points to strong mixing of shell-model configurations by 
the matrices C ^{JTH; a) . 

Although we are very far from a complete theory of the 
TBRE, the evidence presented makes it plausible that 
chaos is a generic property of the shell model. The ma- 
trices C Atl/ ( JTII; a) are both, the fundamental building 
blocks of the TBRE, and the agents for complete mixing 
of the shell-model configurations. 



3. Another Representation of the Shell-Model Hamiltonian 

Can we gain further insight into the structure of 
the matrices C M „(JTIl;ev)? Is it possible to define 
a quantitative measure for the information content 
of shell-model spectra? Why are spin zero ground 
states dominant in the TBRE? These questions are 
answered by transforming the interaction part of the 
shell-model Hamiltonian H^JTTL) in Eq. ([35"]). This 
is done (Papenbrock and Weidenmuller, 2004) by di- 
agonalizing for each set of quantum numbers {JTH} 
the real, symmetric and positive-scmidefinite matrix 
S ap = X>- 1 (jrn)Trace[C(jrn,a)C(jrn;/3)]. We de- 
note by s 2 a > the eigenvalues of S a p, by s a > 
the roots of these eigenvalues, and by O a p the 
eigenvector belonging to the eigenvalue s\. We de- 
fine the new random variables w a = ^2a O a j3V a 
and, for s a > 0, the matrices B^(JTn;a) = 
(l/ s a)S/3 O a pC,j. v {JTTi;{}). By construction, the ma- 
trices B^ v (JTT\]a) are orthonormal with respect to the 
trace, X>- 1 (JTn)Trace[S(jrn, a)B(JTU; j3)} = S a/3 . 



47 



The shell-model Hamiltonian (l35l) takes the form 






Z^ 



<j a s a B py (JTIi;a) 



(44) 

In the last sum over a, only terms pertaining to non-zero 
eigenvalues s„ appear. 

The form (|44]) of the shell-model Hamiltonian is in- 
structive for several reasons. First, by adding further ma- 
trices we could enlarge the set {B( JXTI; a)} of orthonor- 
mal matrices to a complete set of V{ JTU)[D( JTTI)+l]/2 
orthonormal real symmetric matrices. The linear com- 
bination of all these matrices with random coefficients 
would be equivalent to the GOE. The number of matrices 
B(JTU) is very much smaller than T>{JTH)\D{JTU) + 
l]/2, however. This shows once again that the TBRE is 
very different from the GOE; it is a constrained ensemble 
in the sense of Section lV.CI Second, being obtained from 
the v a by an orthogonal transformation, the new random 
variables w a have the same Gaussian distribution as the 
former. However, not all w a but only those pertaining 
to non-vanishing eigenvalues s 2 a do appear in the sum 
in Eq. (|44|) . At least one (and often several) eigenvalues 
s„ always vanish. That implies that one or several lin- 
ear combinations of the v a can never be measured (nor 
do they affect the shell-model spectrum). Third, those 
w a that do appear in Eq. ([44)1 are multiplied with the 
root factors s a . The matrices B(JTH) are orthonormal. 
Therefore, the difficulty of determining the w a in Eq. (|4~4"|) 
from data increases with increasing smallness of the fac- 
tors s a . These root factors are derived from the matrix 
S a p and reflect intrinsic properties of the shell model. 

In Papenbrock and Weidenmuller (2004) the distribu- 
tion of the roots s a has been worked out for the case 
of a single j-shell and for two nuclei in the sd-shell. In 
Fig. [27] we show s a (J) versus J for n — 6 identical nu- 
cleons in the j — 19/2-shell. Very similar results have 
been found in all other cases (different nucleon numbers 
in the same shell and two nuclei in the sd-shell). All 
roots are smooth functions of total spin J. One root is 
substantially bigger than the rest. The associated two- 
body operator is the monopole operator and essentially 
determines the centroid of the shell-model spectrum. 
These features can be understood semi-analytically (Pa- 
penbrock and Weidenmuller, 2004). Identically vanishing 
eigenvalues are related to conserved quantum numbers. 
For instance, the matrix representation of the two-body 
operator J 2 — J(J + 1) (with J the total spin operator) 
can be written as a linear combination of the matrices 
C(JTH;a); that linear combination vanishes identically, 
and similarly for total isospin T. 



4. Preponderance of Spin Zero Ground States in the TBRE 

The discovery by Johnson et al. (1998) of the prepon- 
derance of spin zero ground states used a specific version 
of the TBRE. That version favors neither a particle- 
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FIG. 27 (Color online) The square roots of the eigenvalues 
of the matrix S a p(J) (see text) for the j = 19/2-shell with 6 
identical nucleons versus total spin J. From Papenbrock and 
Weidenmuller (2004). 



particle nor a particle-hole representation, and the au- 
thors refer to it as to the random quasiparticle ensemble. 
Following Johnson et al. (1998), we denote the two-body 
matrix elements (J3J4st\V les \jiJ2St) appearing in Eq. (j34|) 
by Vp'p where p and p' label the two-body states \j1j2st) 
and \j3j4st}. The V p > p have zero mean values and second 
moments given by 



Vp'pVc* 



(2s + l)(2t+l) 



[0 pad p' a' 



v pa' v p' <j\ 



(45) 



Here v 2 is a constant, i.e., independent of all the quan- 
tum numbers. The factor l/(2s + l)(2i + 1) guarantees 
that this a random quasiparticle ensemble. The single- 
particle energies ey in Eq. ([55]) are neglected. 

In Johnson et al. (1998) it was found that for sev- 
eral nuclei in the sd-shell, the probability of finding a 
ground state with spin zero lies between 2/3 and 3/4 al- 
though the total fraction of spin zero states is less than 
10 percent. This result is known as the preponderance of 
spin zero ground states. It holds for the random inter- 
action (|45|) . That interaction does not possess a strong 
pairing force (the agent usually held responsible for the 
preponderance of spin zero ground states in actual nu- 
clei). Other regularities were also observed (Johnson et 
al, 1998, 1999; Zhao et al., 2004b). We confine our- 
selves to the preponderance of spin zero ground states. 
A large number of theoretical papers is devoted to this 
phenomenon. We confine ourselves to the two successful 
explanations that have been offered for the phenomenon 
and refer the reader to the reviews: Zelevinsky and Volya 
(2004); Zhao et al. (2004a) for further references. 

The method used by Zhao and Arima (2001) and re- 
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fined in later papers (Zhao et al., 2002, 2004a) is based 
upon a simple counting procedure. The authors put one 
of the A/" different matrix elements of the residual interac- 
tion equal to (—1) and all others to zero and calculate the 
spectrum. The procedure is repeated Af times, each time 
with a different non-zero matrix element. Let Mj be the 
number of times the ground state is found to have spin J. 
We obviously have J2 j Nj = N '. The probability of find- 
ing a spin zero ground state is then estimated as A/o/A/". 
Comparing the results with an average over many diag- 
onalizations of the TBRE, the authors find good agree- 
ment for a number of cases (four to six fermions in single 
j-shells and two j-shell systems, boson systems). 

The analysis of Papenbrock and Weidenmiiller (2004) 
uses a two-step argument. (For simplicity we replace 
the quantum numbers J, T, II by the single symbol J) . 
First, the square of the width <jj of the spectrum of lev- 
els with quantum number J is defined as the normalized 
variance of the shell-model Hamiltonian. Since shell- 
model spectra have approximately Gaussian shape (Mon 
and French, 1975), <jj has a direct physical interpreta- 
tion. From Eq. (|44| we obtain a 2 j = J2 a w a(J) s a(J)- 
The argument Jonio and on s serves as a reminder that 
both quantities depend on J. That dependence is weak, 
however, at least for the eigenvalues, see Fig. [57] More 
importantly, the spectral widths depend via the w a on 
the same random variables v a and are, therefore, strongly 
correlated: They tend to be all large or all small for a 
given realization of the v a . The variances are biggest 
when many terms contribute almost uniformly to the 
sum over a; Fig. [28] shows that this is the case for <7n. 
Altogether the correlations favor either ern or the spectral 
width of the largest spin. Second, the quantity of interest 
is Rj, the largest (or smallest) eigenvalue of the spectrum 
of levels with spin J. (There is no distinction between 
the two because the v a have random signs). Numeri- 
cal calculations show that the linear relation Rj = rjaj 
holds with rj practically constant (i.e., independent of 
the realization of the v a ). For the states with J = 
this is shown for 6 identical nucleons in the j — 19/2 
shell in the inset in Fig. [28] The dependence of rj on J 
is largely determined by T> (J), the dimension of Hilbert 
space. Typically, T> (J) decreases with increasing J and 
shows odd-even staggering. Both features are present in 
rj, see Fig. [28] These tendencies of rj suppress the com- 
petition of the spectral width belonging to the largest 
spin and favor i?o over all other Rj, see Fig. [55] The ex- 
planation carries over to other cases including nuclei with 
an odd number of nucleons. In the nuclei 20 Ne and 24 Mg, 
it leads to a semiquantitative agreement with numerical 
TBRE calculations (Papenbrock et al., 2006a). 

The method of Papenbrock and Weidenmullcr (2004) 
was improved by Yoshinaga et al. (2006). The au- 
thors considered, for instance, a single j-shell with iden- 
tical nucleons. The energy E(J) of the lowest state of 
spin J was written as E(J) = P(J) _1 Trace(7i) — Qjcrj- 
This equation differs from that of Papenbrock and Wei- 
denmiiller (2004) by the inclusion of the trace of H and 
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FIG. 28 (Color online) 6 identical fermions in a j = 19/2 
shell. Inset: Spectral radius Ro versus spectral width <jq (data 
points) and the linear fit (line). Bottom: Scaling factor rj 
versus J. Top: Probability that the ground state has spin J 
(points) ; probability that spin J has the largest spectral width 
(solid line); probability that the product r.ja.j is maximal 
(dashed line). From Papenbrock and Weidenmiiller (2004). 



by the fact that an analytical form for the function <I>j 
was proposed. The trace of 7i vanishes upon taking 
the ensemble average but, for each realization, fluctu- 
ates around zero. Inclusion of the trace in the equa- 
tion for E(J) removes the scatter of the points around 
the best linear approximation to rj shown in Fig. [55] 
The function $j was fitted and given analytically as 
$j = ^/0.99 \nV(J) + 0.36. Very good agreement is ob- 
tained for m = 4 to m = 6 fermions in several single 
j-shells and systems with two j-shells between TBRE 
results and the predictions based upon this approach. 
The method works also for bosons. 



5. Correlations between Spectra with Different Quantum 
Numbers 



As pointed out at the end of Section lVJ3.il the TBRE 
causes correlations between spectra carrying different 
quantum numbers (A, J, T, II) but belonging to the same 
major shell. From the point of view of the shell model, 
this is not surprising: Switching from one realization of 
the TBRE to another is tantamount to using a differ- 
ent residual interaction; such a different choice of Vi- es 
is bound to affect the spectra of all nuclei in the shell. 
However, the existence of such correlations is surprising 
from the point of view of RMT and exceeds the tradi- 
tional framework of the theory. Moreover, the statistical 
analysis of nuclear data has always assumed the absence 
of correlations between observables with different quan- 
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turn numbers. That assumption is put into question by 
the TBRE. 

In Papenbrock and Weidcnmuller (2006) correlations 
in sd-shell nuclei were displayed between states in the 
same nucleus carrying different quantum numbers and 
between states with identical quantum numbers in dif- 
ferent nuclei. The normalized cross-correlation functions 
had maxima of around 6—10 per cent. Similar figures 
result by calculating correlations in actual shell-model 
spectra by averaging over 17 sd-shell nuclei (Papenbrock 
and Weidenmuller, 2006). 



6. Summary 

The studies of the TBRE reviewed above suggest that 
chaos is a generic property of the shell model. The strong 
mixing is due to the matrices C^(JTJI;a) which, for 
every major shell, are determined by the geometry of 
that shell, by the number of valence nucleons, and by 
the quantum numbers { JTII}. A complete theory of the 
TBRE would have to be based on the analysis of these 
matrices; such an analysis is not available yet. As a con- 
sequence, there is no analytical proof yet that TBRE 
spectra do obey Wigner-Dyson statistics. TBRE and 
GOE are very different; in the TBRE there exist spec- 
tral correlations which are totally absent in the GOE. 



C. Constrained Ensembles 

The spectral fluctuation properties of both, the embed- 
ded ensembles and the TBRE, are known only through 
numerical simulations. Such studies are necessarily con- 
fined to matrices of small dimensions N. A generic an- 
swer (valid for any N) can only be obtained analytically. 
A renewed attempt at such an answer was made in Pa- 
penbrock et al. (2006). The authors introduced and stud- 
ied constrained Gaussian random-matrix ensembles. We 
review here the unitary case. 

Starting point is a complete basis of A*" 2 orthonormal 
Hermitian matrices B a in the A-dimensional Hilbert 
space. Orthonormality is defined in terms of the trace, 



(B a \B p )= Trace(B a B ) = 8 a t 



(46) 



We note that in contrast to Section IV.B.3I we do not 
include the matrix dimension in the definition (|46|l . Any 
Hermitian matrix H can be expanded in terms of that 
basis, 
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H — 2_^ h a B a 



(47) 



Q = l 



Taking the N 2 complex coefficients h a in Eq. (|47| as un- 
corrected Gaussian-distributed random variables with 
zero mean value and a common second moment, one ob- 
tains the GUE. A constrained ensemble is obtained by 



requiring a subset {h q } of the expansion coefficients h a 
in Eq. (|4"T|) to vanish identically. (It was shown in Sec- 
tion lV.B.3l that the TBRE can be viewed as a constrained 
ensemble, and it is easy to see that the same statement 
holds for the embedded ensembles). The resulting ensem- 
ble is not unitarily invariant, however. Unitary invari- 
ance is restored by integrating the constraining condition 
Yl q S((B q \H)) over the unitary group, see Section Hl.BI 
This yields for the probability density of the matrices H 



W(H) ex exp 



N 
2A 2 " 



(H\H) 



J 6u(j[5((UB q W\H))\ . (48) 



The first factor on the right-hand side is the same as in 
Eq. (fi"0|) . The second factor represents the constraints. 
It is unitarily invariant and, therefore, does not affect the 
eigenvector distribution. The eigenvalue distribution of 
the constrained ensemble differs from that on the right- 
hand side of Eq. (|12|) by the additional factor 
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F(E 1 ,...,E N ) = J dU\T[6({B q \UEU*)) 



Here E is the diagonal matrix of eigenvalues. 

In the GUE, quadratic level repulsion is a consequence 
of the Vandermonde determinant IIp<cr \E p — E a \ appear- 
ing on the right-hand side of Eq. (fT2"]) . Level repulsion 
is lifted by the constraints if and only if F{E\, . . . , En) 
is singular whenever two eigenvalues coincide. Using 
Fourier transformation, one replaces each of the delta 
functions by an integral over a plane wave. The function 
F(Ei, . . . , En) can then be written as a Harish-Chandra 
Itzykson Zuber integral (Harish-Chandra 1957, Itzykson 
and Zuber, 1980). Inspection of the result of the integra- 
tion yields a sufficient condition for F(E\, . . . , En) not 
to be singular. In its simplest form, that condition reads 



Nq< 



N{N - 1) 



(50) 



Here Nq is the number of constraints. It is quite re- 
markable that the sufficient condition depends only on 
Nq. 

The result can be extended to the GOE and the GSE. 
Moreover, other properties (spectral radius, distribution 
of matrix elements) of the constrained ensembles can 
be investigated and are seen to differ from those of the 
canonical ensembles. Relations can be established be- 
tween the ensemble defined by the constraints on the 
set {h q } of matrix elements, and the one defined by the 
complementary set (all h a but the set {h q } are con- 
strained). Unfortunately, the condition ([50)) does not 
cover the physically interesting cases of the TBRE and of 
the embedded ensemble with two-body interactions. The 
spectral fluctuation properties of these ensembles remain 
an open theoretical problem. 
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VI. SUMMARY AND CONCLUSIONS 

Comparison of GOE predictions with data reviewed in 
Section HTT1 shows that there is good agreement on spec- 
tral fluctuations not only near neutron threshold but also 
in the ground-state domain. This is true also for the 
three nuclei with complete level schemes. Extensions of 
the GOE are used to successfully describe isospin mixing, 
and to test time-reversal invariance. Deviations from 
GOE predictions are found up to several 100 keV above 
yrast and in case of symmetries. An enlargement of the 
data basis would be highly desirable. There is hope that 
the next generation of large-scale gamma detectors will 
contribute pertinent information on the interesting en- 
ergy interval 500 keV to 1 MeV above yrast. 

Studies of the spherical shell model, of the Nilsson 
model, and of the interacting boson model reviewed in 
Section [IVI show a joint tendency towards strong mixing 
of the unperturbed configurations and towards GOE fluc- 
tuation properties. This resolves the apparent dichotomy 
between Bohr's picture of the compound nucleus de- 
scribed in Sect ion HTland the independent -particle model. 
Typically, the GOE limit is not fully attained, however. 
For the sd-shell, mixing is strongest for nuclei in the 
middle of the shell but even here the spreading width 
of the unperturbed shell-model configurations located in 
the center of the spectrum is 20 MeV or so and, thus, 
smaller than the range of the spectrum. 

Random-matrix ensembles patterned after nuclear- 
structure concepts were reviewed in Section [V] The em- 
bedded two-body random ensemble possesses strong pre- 
dictive power for average properties such as the shape of 
the spectrum or the distribution of transition strengths. 
This is true even though that ensemble disregards con- 
served quantum numbers like spin or isospin. The two- 
body random ensemble is closest in structure to the shell 
model. The strong mixing of shell-model configurations 
is seen to be a generic feature of the shell model. At 
the same time, the two-body ensemble possesses features 
that go beyond the GOE such as the correlations between 
spectra with different quantum numbers. It would be in- 
triguing to find such correlations experimentally. 

Following the Bohigas-Giannoni-Schmit conjecture, 
we have identified chaos in nuclei with the occurrence of 
spectral fluctuations of the GOE type. Closer inspection 
has shown, however, that chaos in nuclei differs from that 
in few-degrees-of-freedom systems. The evidence pre- 
sented shows that chaos is caused by the residual two- 
body interaction of the shell model. Shell structure is 
a typical feature of fermionic many-body systems. The 
residual interaction strongly mixes the shell-model con- 
figurations within a major shell, but leaves the overall 
shell structure largely intact. Evidence for this assertion 
comes not only from the success of shell-model calcula- 
tions in the ground-state domain, but also from the suc- 
cess of the optical model for elastic scattering and from 
the existence of distinct maxima in the neutron strength 
function. (The last two features are mentioned with- 



out explanation for the sake of the argument and will be 
treated in part 2 of this review) . The existence of regular 
features in nuclei also attests to the incomplete mixing of 
the shell-model configurations. We mention the regular- 
ities in the ground-state domain caused by the pairing 
force, collective motion, and doorway states. 

The evidence presented in this review strongly sup- 
ports the view that chaos is a generic property of nuclei. 
At the same time, we are far from having a complete theo- 
retical understanding of chaos in nuclei. First, we lack an 
overall many-body theory that would permit the calcu- 
lation of nuclear spectra from the nucleon-nucleon inter- 
action within controlled approximations (except for the 
lightest nuclei). In contrast to few-degrees-of-freedom 
systems, there is also no theoretical framework such as 
the semiclassical approximation which would establish 
the connection between classical chaos and spectral fluc- 
tuation properties of the RMT type. On a more mun- 
dane level, the evidence presented above for chaos in the 
spherical shell model comes mainly from the sd-shell. Al- 
though the arguments seem generic, it would be gratify- 
ing to have similar evidence in other major shells. Would 
the sizeable correlations between spectra carrying differ- 
ent quantum numbers found in the two-body random 
ensemble persist with increasing shell size? Likewise, it 
is desirable to attain a deeper analytical understanding of 
how chaos arises within a major shell. Does the difference 
between the GOE and the two-body random ensemble 
entail other differences beyond those correlations? Also, 
the analysis of properties of the matrices C^ v in Eq. (j35|) 
within a single j -shell would be of substantial interest, 
coupled, if possible, with the proof of GOE spectral fluc- 
tuation properties for j>l, 

Chaos limits the predictability of spectral properties in 
nuclei in terms of simple models. For instance, the shell- 
model eigenfunctions become mixed ever more strongly 
as the excitation energy increases, see Section IIV.A.2I 
Such strong mixing requires the presence of many matrix 
elements of the two-body interaction and cannot be mod- 
eled in simple terms. We are not aware of any attempts 
to formulate and quantify that limitation, however. 

The way it is defined in this review, chaos is a statis- 
tical property of levels carrying identical quantum num- 
bers. Regular dynamical features in nuclei typically re- 
late several states carrying different quantum numbers. 
The coexistence of those two aspects of nuclear mo- 
tion deserves deeper analysis. Hopefully, it would also 
shed light on the dynamical properties of other fermionic 
many-body systems. 



Acknowledgments 

We owe thanks to many colleagues with whom we have 
discussed topics relating to chaos in nuclei over many 
years, and from whom we have learned much. Special 
and cordial thanks are due to those who took the trouble 
to read the first draft of this paper, and who made numer- 



51 



ous valuable and thoughtful comments and suggestions: 
S. Aberg, Y. Alhassid, O. Bohigas, T. Friedrich, G. Gar- 
vey, T. Guhr, T. Papcnbrock, A. Richter, T. H. Seligman, 
V. Zelevinsky, Y. M. Zhao. We would also like to thank 
R. Chankova, A. Chyzh, H. Friedrich, and J. F. Shriner, 
Jr., for assistance with the figures and manuscript prepa- 
ration. This work was supported in part by the U.S. 
Department of Energy grant DE-FG02-97ER41042. 



References 



[9 
[10 

[11 

[12; 

[is: 

[i4: 
[is: 

[16" 

[17; 
[is; 

[19 

[20; 

[21 

[22[ 

[23; 

[24 
[25: 

[26" 
[27 



Aberg, S., 1990, Phys. Rev. Lett. 64, 3119. 
Aberg, S., 1992, Prog. Nucl. Part. Phys. 28, 11. 
Abul-Magd A. Y., and H. A. Weidenmiiller, 1985, Phys. 
Lett. B 162, 223. 

Abul-Magd, A. Y., H. L. Harney, M. H. Simbel, and H. 
A. Weidenmiiller, 2004 Phys. Lett. B 579, 278. 
Adams, A. A., G. E. Mitchell, and J. F. Shriner, Jr., 
1998, Phys. Lett. B 422, 13. 

Agassi, D., C. M. Ko, and H. A. Weidenmiiller, 1977, 
Ann. Phys. (N.Y.) 107, 140. 

Agvaanluvsan, U., G. E. Mitchell, J. F. Shriner, Jr., and 
M. P. Pato, 2003, Nucl. Instrum. Methods Phys. Res. 
A 498, 459. 

Alhassid, Y., A. Novoselsky, and N. Whelan, 1990, 
Phys. Rev. Lett. 65, 2971. 

Alhassid, Y., and N. Whelan, 1991, Phys. Rev. Lett. 
67, 816. 

Alhassid, Y., and N. Whelan, 1993, Phys. Rev. Lett. 
70, 572. 

Alt, H., H.-D. Graf, H. L. Harney, R. Hofferbest, H. 
Lengeler, A. Richter, P. Schaardt, and H. A. Wei- 
denmiiller, 1995, Phys. Rev. Lett. 74, 62. 
H. Alt, G I. Barbosa, H.-D. Graf, T. Guhr, H. L. Har- 
ney, R. Hofferbert, H. Rehfeld, A. Richter, 1998, Phys. 
Rev. Lett. 81, 4847. 

Altland, A. and M. R. Zirnbauer, 1997, Phys. Rev. B 
55, 1142. 

Arima, A., and F. Iachello, 1975, Phys. Rev. Lett. 35, 
1069. 

Arvieu, R., F. Brut, J. Carbonell, and J. Touchard, 
1987, Phys. Rev. A 35, 2389. 
Axel, P., 1962, Phys. Rev. 126, 671. 
Bae, M. S., T. Otsuka, T. Misusaki, and N. Fukunishi, 
1992, Phys. Rev. Lett. 69, 2349. 

Balian, R., 1968, Nuovo Cim. Soc. Ital. Fis. B 57, 183. 
Balian, R., and C. Bloch, 1970, Ann. Phys. (N.Y.) 60, 
401. 

Baltes, H. P., and E. R. Hilf, 1976, Spectra of Finite 
Systems (Wissenschaftsverlag, Mannheim). 
Barbosa, C. I., T. Guhr, and H. L. Harney, 2000, Phys. 
Rev. E 62, 1936. 

Benet, L., and H. A. Weidenmiiller, 2003, J. Phys. A: 
Math. Gen. 36, 3569. 

Berry, M. V., and M. Tabor, 1977, Proc. Roy. Soc. Lon- 
don A 356, 375. 

Berry, M. V., 1981, Ann. Phys. (N.Y.) 131, 163. 
Berry, M. V., and M. Robnik, 1984, J. Phys. A: Math. 
Gen. 17, 2413. 

Berry, M. V., 1985, Proc. Roy. Soc. London A 400, 229. 
Bertulani, C. A., and V. Yu. Ponomarev, 1999, Phys. 
Rep. 321, 139. 



[2s; 
[29; 

[30 
[31 

[32; 
[33; 

[34; 
[35; 

[36; 

[37 

[38" 
[39" 

[40; 

[41 

[42[ 
[43; 
[44; 
[45; 

[4e; 

[47 

[4s; 
[49; 

[50; 

[51 

[52; 
[53; 
[54 
[55; 

[5e; 

[57 

[5s; 



Bilpuch, E. G, A. M. Lane, J. D. Moses, and G E. 
Mitchell, 1976, Phys. Rep. 28, 145. 
Blocki, J., J. Skalski, and W. J. Swiatecki, 1995, Nucl. 
Phys. A 594, 137. 

Bohigas, O., and J. Flores, 1971, Phys. Lett. B 34, 261. 
Bohigas, O., R. U. Haq, and A. Pandey, 1983, in: Nu- 
clear Data for Science and Technology, edited by K. H. 
Bockhoff (Reidel, Dordrecht), p. 809. 
Bohigas, O., M. J. Giannoni, and C. Schmit, 1984, Phys. 
Rev. Lett. 52, 1. 

Bohigas, O., and M. J. Giannoni, 1984, Lecture Notes 
in Physics 209 (Springer, Heidelberg). 
Bohigas, O., and H. A. Weidenmiiller, 1988, Ann. Rev. 
Nucl. Part. Science 38, 421. 

Bohigas, O., and P. Leboeuf, 2002, Phys. Rev. Lett. 88, 
092502. 

Bohigas, O., and M. P. Pato, 2004, Phys. Lett. B 595, 
171. 

Bohle, D., A. Richter, W. Steffen, A. E. L. Dieperink, 
N. Lo Iudice, F. Palumbo, and O. Scholten, 1984, Phys. 
Lett. B 137, 27. 

Bohr, N., 1936, Nature 137, 344. 
Bohr, A., 1951, Phys. Rev. 81, 134. 
Bohr, A., 1952, K. Dan. Vidensk. Selsk. Mat. Fys. 
Medd. 26, No. 14. 

Bohr, A., and B. R. Mottelson, 1952, K. Dan. Vidensk. 
Selsk. Mat. Fys. Medd. 27, No. 16. 
Bohr, A., and B. Mottelson, Nuclear Structure, W. A. 
Benjamin, New York, Vol. 1 (1969) and Vol. 2 (1975). 
Bollinger, L. M., and G E. Thomas, 1968, Phys. Rev. 
C 171, 1293. 

Bollinger, L. M., and G E. Thomas, 1970, Phys. Rev. 
C 2, 1951. 

Brack, M., J. Damgaard, A. S. Jensen, H. C. Pauli, V. 
M. Strutinsky, and G Y. Wong, 1972, Rev. Mod. Phys. 
44, 320. 

Brack, M., and R. K. Bhaduri, 1997, Semiclassical 
Physics (Addison Wesley, New York). 
Brink, D., 1955, D. Phil. Thesis, Oxford University, un- 
published. 

Brink, D., J. Neto, and H. A. Weidenmiiller, 1978, Phys. 
Lett. B 80, 170. 

Brody, T. A., J. Flores, J. B. French, P. A. Mello, A. 
Pandey, and S. S. M. Wong, 1981, Rev. Mod. Phys. 53, 
385. 

Brown, B. A., and G Bertsch, 1984, Phys. Lett. B 148, 
5. 

Brown, B. A., and B. H. Wildenthal, 1988, Ann. Rev. 
Nucl. Part. Sci. 38, 29. 

Bulgac, A., and D. Kusnezov, 1996, Phys. Rev. E 54, 
3468. 

Bunakov, V. E., 1999, Yadernaya Fisica 62, 5 [Phys. 
Atomic Nucl. 62, 1 (1999)]. 

Caldeira, A. O., and A. J. Leggett, 1983, Ann. Phys. 
(N.Y.) 149, 374. 

Carlson, B. V., L. F. Canto, S. Cruz-Barrios, M. S. 
Hussein, and A. F. R. de Toledo Piza, 1999, Ann. Phys. 
(N.Y.) 276, 111. 

Carlson, B. V., M. S. Hussein, A. F. R. de Toledo Piza, 
and L. F. Canto, 1999, Phys. Rev. C 60, 014604. 
Casati, G, F. Valz-Gris, and I. Guarneri, 1980, Lett. 
Nuov. Cim. 28, 279. 

Cassing, W., and W. Norenberg, 1985, Nucl. Phys. A 
433, 467. 



52 



[59] Caurier, E., G. Martinez-Pinedo, F. Nowacki, J. Reta- [87 

mosa, and A. P. Zuker, 1999, Phys. Rev. C 59, 2033. 
[60] Caurier, E., G. Martinez-Pinedo, F. Nowacki, A. Poves, 

and A. P. Zuker, 2005, Rev. Mod. Phys. 77, 427. [. 

[61] Chevalley, C, 1946, Theory of Lie Groups, (Princeton 

University Press, Princeton). [89 

[62] Conway, J., 1990, A Course in Functional Analysis, 

(Springer, New York). [90 

[63] Dembowski, C, B. Dietz, T. Friedrich, H.-D. Graf, H. 

L. Harney, A. Heine, M. Miski-Oglu, A. Richter, 2005, [91 

Phys. Rev. E 71, 046202. 
[64] De Pace, A., A. Molinari, and H. A. Weidenmiiller, [92 

2007, Ann. Phys. (N.Y.) 322, 2446. 
[65] De Shalit, A., and I. Talmi, 1963, Nuclear Shell Theory [93; 

(Academic, New York). 
[66] Dias, H., M. S. Hussein, N. A. de Oliveira, and B. H. [94 

Wildenthal, 1989, J. Phys. G: Nucl. Part. 15, L79. 
[67] Dieperink, A. E. L., O. Scholten, and F. Iachello, 1980, [95 

Phys. Rev. Lett. 44, 1747. 
[68] Dietz, B., T. Guhr, H. L. Harney, and A. Richter, 2006, [96 

Phys. Rev. Lett. 96, 254101. 
[69] D0ssing, T., B. Herskind, S. Leoni, A. Bracco, R. A. [97 

Broglia, M. Matsuo, and E. Vigezzi, 1996, Phys. Rep. 

268, 1. [' 

[70] D0ssing, T., A. P. Lopez-Martens, T. L. Khoo, T. Lau- 

ritsen, and S. Aberg, 2004, in The Labyrinth in Nuclear [99 

Structure, edited by A. Braco and C. A. Kalfas, (AIP, [100 

Melville, NY), p. 164. [101 

[71] Duflo, J., and A. P. Zuker, 1995, Phys. Rev. C 52, R23. 
[72] Dyson, F. J., and M. L. Mehta, 1963, J. Math. Phys. 4, [102 

701. 
[73] von Egidy, T., A. N. Behkami, and H. H. Schmidt, 1986, [103 

Nucl. Phys. A 454, 109; ibid. 1989, A 481, 189. 
[74] Ellegaard, C, T. Guhr, K. Lindeman, J. Nygard, and [104 

M. Oxborrow, 1996, Phys. Rev. Lett. 77, 4918 . 
[75] Enders, J., T. Guhr, N. Huxel, P. von Neuman-Cosel, 

C. Rangacharyulu, and A. Richter, 2000, Phys. Lett. B [105 

486, 273. 
[76] Enders, J., T. Guhr, A. Heine, P. von Neuman-Cosel, [106; 

V. Yu. Ponomarcv, A. Richter, and J. Wambach, 2004, 

Nucl. Phys. A 741, 3. [107 

[77] Fermi, E., E. Amaldi, O. D'Agostino, F. Rasetti, and 

E. Segre, 1934, Proc. Roy. Soc. A 146, 483; ibid. 149 [108 

522 (1935). 
[78] Flambaum, V. V., A. A. Gribakina, G. F. Gribakin, and [109 

M. G. Kozlov, 1994, Phys. Rev. A 50, 267. 
[79] Frazier, N., B. A. Brown, and V. Zelevinsky, 1996, Phys. [110 

Rev. C 54 1665. 
[80] French, J. B., and S. S. M. Wong, 1970, Phys. Lett. B [111 

33, 449. 
[81] French, J. B., V. K. B. Kota, A. Pandey, and S. Tomso- [112 

vie, 1985, Phys. Rev. Lett. 54, 2313. 
[82] Garg, J. B., J. Rainwater, J. S. Petersen, and W. W. [113 

Havens, Jr., 1964, Phys. Rev. 134, B985. 
[83] Garrett, J. D., J. Q. Robinson, A. J. Foglia, and H.-Q. [114 

Jin, 1997, Phys. Lett. B 392, 24. 
[84] Gaudin, M., 1961, Nucl. Phys. 25, 447. [115; 

[85] A. Gavron, A. Gayer, J. Boissevain, H. C. Britt, J. 

R. Nix, A. J. Sierk, P. Grang, S. Hassani, H. A. Wei- [116; 

dennriiller, J. R. Beene, B. Cheynis, D. Drain, R. L. 

Ferguson, F. E. Obenshain, F. Plasil, G. R. Young, G. 

A. Petitt and C. Butler, 1986, Phys. Lett. B 176, 312. [117 

[86] Grange, P., and H. A. Weidenmiiller, 1980, Phys. Lett. 

B 96, 26. 



Grossmann, C. A., M. A. LaBonte, G. E. Mitchell, J. D. 

Shriner, J. F. Shriner, Jr., G. A. Vavrina, P. M. Wallace, 

2000, Phys. Rev. C 62, 024323. 

Gu, J. Z., and H. A. Weidenmiiller, 1999, Nucl. Phys. 

A 660, 197. 

Gu, J. Z., and H. A. Weidenmiiller, 2001, Nucl. Phys. 

A 690, 382. 

Guhr, T., and H. A. Weidenmiiller, 1990a, Ann. Phys. 

(N.Y.) 199, 412. 

Guhr, T., and H. A. Weidenmiiller, 1990b, Chem. Phys. 

146, 21. 

Guhr, T., A. Miiller-Groeling, and H. A. Weidenmiiller, 

1998, Phys. Rep. 299, 189. 

Gurevich, I. I., and M. I. Pevsner, 1957, Nucl. Phys. 2, 
575. 

Gutzwiller, M., 1990, Chaos in Classical and Quantum 
Mechanics (Springer, Berlin). 

Haake, F., 2001, Quantum Signatures of Chaos, 2nd edi- 
tion (Springer, Berlin). 

Hackenbroich, G., and H. A. Weidenmiiller, 1995, Phys. 
Rev. Lett. 74, 4118. 

Hamoudi, A., R. G. Nazmitdinov, E. Shahaliev, and Y. 
Alhassid, 2002, Phys. Rev. C 65, 064311. 
Haq, R. U., A. Pandey, and O. Bohigas, 1982, Phys. 
Rev. Lett. 48, 1086. 

Harish-Chandra, 1957, Am. J. Math. 79, 87. 
Harney, H. L., 1984, Z. Phys. A 316, 177. 
Harney, H. L., A. Richter, and H. A. Weidenmiiller, 

1986, Rev. Mod. Phys. 58, 607. 

Haxel, O., J. H. D. Jensen, and H. E. Suess, 1949, Phys. 

Rev. 75, 1766. 

Heiss, W. D., R. G. Nazmitdinov, and S. Radu, 1994, 

Phys. Rev. Lett. 72, 2351. 

Herskind, B., B. Lauritzen, K. Schiffer, R. A. Broglia, 

F. Barranco, M. Gallardo, J. Dudek, and E. Vigezzi, 

1987, Phys. Rev. Lett. 59, 2416. 

Heusler, S., S. Midler, A. Altland, P. Braun, and F. 

Haake, 2007, Phys. Rev. Lett. 98, 044103. 

Hofmann, H., and P. J. Siemens, 1977, Nucl. Phys. A 

275, 464. 

Hofmann, H., F. A. Ivanyuk, C. Rummel, and S. Ya- 

maji, 2001, Phys. Rev. C 64, 054316. 

Honma, M., T. Otsuka, B. A. Brown, and T. Mizusaki, 

2002, Phys. Rev. C 65, 061301. 

Hussein, M. S., and M. P. Pato, 2000, Phys. Rev. Lett. 

84, 3873. 

Iachello, F., and I. Talmi, 1987, Rev. Mod. Phys. 59, 

339. 

Iachello, F., and A. Arima, 1987, The Interacting Boson 

Model (Cambridge Press, Cambridge). 

Imry, Y., 2002 Introduction to Mesoscopic Physics, 2nd 

edition, (Oxford Press, Oxford). 

Itzykson, C, and J. B. Zuber, 1980, J. Math. Phys. 21, 

411. 

Johnson, C. W., G. F. Bertsch, and D. J. Dean, 1998, 

Phys. Rev. Lett. 80, 2749. 

Johnson, C. W., G. F. Bertsch, D. J. Dean, and I. Talmi, 

1999, Phys. Rev. C 61, 014311. 

Jolie, J., R. F. Casten, P. Cejnar, S. Heinze, E. A. Mc- 

Cutchan, and N. V. Zamfir, 2004, Phys. Rev. Lett. 93, 

132501. 

Koehler, P. E., J. L. Ullmann, T. A. Bredeweg, J. M. 

ODonnell, R. Reifarth, R. S. Rundberg, D. J. Vieira, 

and J. M. Wouters, 2007, Phys. Rev. C 76, 025804. 



r>: J > 



[118] Kota, V. K. B., 2001, Phys. Rep. 347, 223. 

[119] Kota, V. V. B., 2007 J. Math. Phys. 48, 053304. 

[120] Kruecken, R., A. Dewald, P. von Brentano, and H. A. 
Weidenmiiller, 2001, Phys. Rev. C 64, 064301. 

[121] Lane, A. M., and R. G. Thomas, 1958, Rev. Mod. Phys. 
30, 257. 

[122] Leboeuf, and P., J. Roccia, 2006, Phys. Rev. Lett. 97, 
010401. 

[123] Leviandier, L., M. Lombardi, R. Jost, and J. P. Pique, 
1986, Phys. Rev. Lett. 56, 2449. 

[124] Lewenkopf, C, and V. Zelevinsky, 1994, Nucl. Phys. A 
569, 183c. 

[125] Liou, H.I., H. S. Camarda, S. Wynchank, M. Slagowitz, 
G. Hacken, F. Rahn, and J. Rainwater, 1972, Phys. Rev. 
C 5, 974. 

[126] Liou, H. I., H. S. Camarda, and F. Rahn, 1972, Phys. 
Rev. C 5, 1002. 

[127] Lombardi, M., and T. H. Seligman, 1993, Phys. Rev. A 
47, 3571 (1993). 

[128] Lombardi, M., O. Bohigas, and T. Seligman, 1994, 
Phys. Lett. B 324, 263. 

[129] J. E. Lynn, 1968, The Theory of Neutron Resonances 
(Clarendon, Oxford). 

[130] Mahaux, C, and H. A. Weidenmiiller, 1969, Shell- 
Model Approach to Nuclear Reactions (North-Holland, 
Amsterdam). 

[131] Matsuo, M., T. D0ssing, E. Vigezzi, and S. Aberg, 1997, 
Nucl. Phys. A 620, 296. 

[132] Martinez-Pinedo, G., A. P. Zuker, A. Povcs, and E. 
Caurier, 1997, Phys. Rev. C 55, 187. 

[133] Mayer, M. G., 1949, Phys. Rev. 75, 1969. 

[134] Mayer, M. G., and J. H. D. Jensen, 1955, Elementary 
Theory of Nuclear Shell Structure (Wiley, New York). 

[135] McDonald, S. W., and A. N. Kaufman, 1979, Phys. Rev. 
Lett. 42, 1189. 

[136] Mehta, M. L., 2004, Random Matrices, 3rd edition (Aca- 
demic Press, New York). 

[137] Meredith, D. C, S. E. Koonin, and M. R. Zirnbauer, 
1988, Phys. Rev. A 37, 3499. 

[138] Meredith, D. C, 1993, Phys. Rev. E 47, 2405. 

[139] Mitchell, G. E., E. G. Bilpuch, J. F. Shriner, Jr., and 
A. M. Lane, 1985, Phys. Rep. 117, 1. 

[140] Mitchell, G. E., and J. F. Shriner, Jr., 2001, Physica 
Scripta T 90, 105. 

[141] Moller, P., J. R. Nix, W. D. Myers, and W. J. Swiatecki, 
1995, At. Data Nucl. Data Tables 59, 185. 

[142] Molinari, A., and H. A. Weidenmiiller, 2006, Phys. Lett. 
B 637, 48. 

[143] Mon, K. F., and J. B. French, 1975, Ann. Phys. (N.Y.) 
95, 90. 

[144] Mottelson, B. R., 1992, Nucl. Phys. A 557, 717c. 

[145] Nature, 1936, 137, 351. 

[146] Neutron Cross Section, Sigma Center, Brookhaven Na- 
tional Laboratory BNL 325, Suppl. 2, Brookhaven, N. 
Y. (1964). 

[147] Norenberg, W., and H. A. Weidenmiiller, 1980, Intro- 
duction to the Theory of Heavy-Ion Reactions, 2nd edi- 
tion (Springer, Berlin). 

[148] Olofsson, H., Aberg, O. Bohigas, and Leboeuf, 2006, 
Phys. Rev. Lett. 96, 042502. 

[149] Ormand, W. E., and R. A. Broglia, 1992, Phys. Rev. C 
46, 1710. 

[150] Pandey, A., 1981, Ann. Phys. (N. Y) 134, 110. 

[151] Papenbrock, T., and H. A. Weidenmiiller, 2004, Phys. 



152; 

153 
154 
155 

156 

157; 

158; 

159 
160 

161 
162 

163 

164 

165 

166 

167 

168; 

169 

170 

171 

172 
173 

174 

175 

176; 

177' 
178 : 
179^ 

180 
181 
182 
183 



Rev. Lett. 93, 132503. 

Papenbrock, T., and H. A. Weidenmiiller, 2005, Nucl. 
Phys. A 757, 422. 

Papenbrock, T., Z. Pluhar, and H. A. Weidenmiiller 
2006, J. Phys. A: Math. Gen. 39, 9709. 
Papenbrock, T., and H. A. Weidenmiiller, 2006, Phys. 
Rev. C 73, 014311. 

Papenbrock, T., and H. A. Weidenmiiller, 2007, Rev. 
Mod. Phys. 79, 997. 

Percival, I. C, 1973, J. Phys. B 6, L229. 
Pieper, S. C, and R. B. Wiringa, 2001, Ann. Rev. Nucl. 
Part. Sci. 51, 53. 

Porter, C.E., 1965, Statistical Theories of Spectra: Fluc- 
tuations, Academic Press, New York. 
Ragnarsson, I., and S. Aberg, 1986, Phys. Lett. B 180, 
191. 

Raman, S., T. A. Walkiewicz, S. Kahane, E. T. Jurney, 
J. Sa, Z. Gacsi, J. L. Weil, K. Allart, G. Bonsignori, and 
J. F. Shriner, Jr., 1991, Phys. Rev. C 43, 521. 
Richter, A., 1995, Prog. Part. Nucl. Phys. 34, 213. 
Roccia, J., and P. Leboeuf, 2007, Phys. Rev. C 76, 
014301. 

Rosenzweig, N., and C. E. Porter, 1960, Phys. Rev. 120, 
1698. 

Sakhr, J., and N. D. Whelan, 2003, Phys. Rev. E 67, 
066213. 

Samyn, M., S. Goriely, M. Bender, and J. M. Pearson, 
2004, Phys. Rev. C 70, 04439. 

Sargeant, A. J., M. S. Hussein, and A. N. Wilson, 2005, 
in Nuclei and Mesoscopic Physics, East Lansing, edited 
by V. Zelevinsky, (AIP, Melville, NY), p. 46. 
Shimizu, Y. R., F. Barranco, R. A. Broglia, T. D0ssing, 
and E. Vigezzi, 1992, Phys. Lett. B 274, 253. 
Shriner, Jr., J. F., G. E. Mitchell, and E. G. Bilpuch, 
1987, Phys. Rev. Lett. 59, 435. 

Shriner, Jr., J. F., E. G. Bilpuch, and G. E. Mitchell, 
1989, Z. Phys. A 332, 45. 

Shriner, Jr., J. F., G. E. Mitchell, and T. von Egidy, 
1991, Z. Phys. A 338, 309. 

Shriner, Jr., J. F., C. A. Grossmann, and G. E. Mitchell, 
2000, Phys. Rev. C 62, 054305. 
Soloviev, V. G., 1995, Nucl. Phys. A 586, 265. 
Sommermann, M., and H. A. Weidenmiiller, 1993, Eu- 
rophys. Lett. 23, 79. 

Stafford, C. A., and B. R. Barrett, 1999, Phys. Rev. C 
60, 051305. 

Stephens, F. S., et al, 2005, Phys. Rev. Lett. 94, 
042501. 

Strutinsky, V. M., 1966 Yad. Fiz. 3, 614; Sov. Journ. 
Nucl. Phys. 3, 499 (1966). 

Strutinsky, V. M., 1967, Nucl. Phys. A 95, 420. 
Strutinsky, V. M., 1968, Nucl. Phys. A 122, 1. 
Twin, P. J., B. M. Nyako, A. H. Nelson, J. Simpson, M. 
A. Bentley, H. W. Cranmer-Gordon, P. D. Forsyth, D. 
Howe, A. R. Mokhtar, J. D. Morrison, J. F. Sharpey- 
Schafer, G. Sletten, 1986, Phys. Rev. Lett. 57, 811. 
Verbaarschot, J. J. M., and P. J. Brussard, 1979, Phys. 
Lett. B 87, 155. 

Vigezzi, E., R. A. Broglia, and T. D0ssing, 1990, Nucl. 
Phys. A A 520, 179c. 

Vigezzi, E., R. A. Broglia, and T. D0ssing, 1990, Phys. 
Lett. B 249, 163. 

Vretenar, D., and Y. Alhassid, 1992, Phys. Rev. C 46, 
1334. 



54 



[184] Watson, III, W. A., E. G. Bilpuch, and G. E. Mitchell, 

1981, Z. Phys. A 300, 89. [194 

[185] Weidenmiiller, H. A., 1980, Progress in Particle and Nu- 
clear Physics 3, 49. [195 

[186] Weidenmiiller, H. A., 1985, in: Nuclear Structure 1985, 

edited by R. Broglia, G. B. Hagemann, B. Herskind (El- [196 
sevier, Amsterdam) p. 213. 

[187] Weidenmiiller, H. A., 1993, Phys. Rev. A 48, 1819. [197 

[188] Weidenmiiller, H. A., P. von Brentano, and B. R. Bar- [198 : 
rett, 1998, Phys. Rev. Lett. 81, 3603. 

[189] Whitehead, R. R., A. Watt, D. Kelvin, and A. Conkie, [199 

1978, Phys. Lett. B 76, 149. 

[190] Wilczynski, J., 1973, Phys. Lett. B 47, 484. [200 

[191] Wilson, W. M., E. G. Bilpuch, and G. E. Mitchell, 1975, 

Nucl. Phys. A 245, 285. [201 

[192] Wildenthal, B. H., 1984, Prog. Part. Nucl. Phys. 11, 5. 

[193] Yoshinaga, N., N. Yoshida, T. Shigehara, and T. Cheon, 



1993, Phys. Rev. C 48, R509. 

Yoshinaga, N., A. Arima, and Y. M. Zhao, 2006, Phys. 

Rev. C 73, 017303. 

Zelevinsky, V., B. A. Brown, N. Frazier, and M. Horoi, 

1996, Phys. Rep. 276, 85. 

Zelevinsky, V., and A. Volya, 2004, Phys. Rep. 391, 

311. 

Zelevinsky, V., 2007, private communication. 

Zhao, Y. M., and A. Arima, 2001, Phys. Rev. C 64, 

041301. 

Zhao, Y. M., A. Arima, and N. Yoshinaga, 2002, Phys. 

Rev. C 66, 034302. 

Zhao, Y. M., A. Arima, and N. Yoshinaga, 2004, Phys. 

Rep. 400, 1. 

Zhao, Y. M., A. Arima, N. Shimizu, K. Ogawa, N. 

Yoshinaga, and O. Scholten, 2004, Phys. Rev. C 70, 

054322. 



